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We construct the generic phase diagrams encoding the topologically distinct localized and delo¬ 
calized phases of noninteracting fermionic quasiparticles for any symmetry class from the tenfold 
way in one, two, and three dimensions. To this end, we start from a massive Dirac Hamiltonian per¬ 
turbed by a generic disorder for any dimension of space and for any one of the ten symmetry classes 
from the tenfold way. The physics of Anderson localization is then encoded by a two-dimensional 
phase diagram that we deduce from the topology of the space of normalized Dirac masses. This 
approach agrees with previously known results and gives an alternative explanation for the even-odd 
effect in the one-dimensional chiral symmetry classes. We also give a qualitative explanation for the 
Gade singularity and Griffiths effects in the density of states using the first homotopy group of the 
normalized Dirac masses in two dimensions. Finally, this approach is used to analyze the stability 
of massless Dirac fermions on the surface of three-dimensional topological crystalline insulators. 


I. INTRODUCTION 

Anderson localization is the nonperturbative phe¬ 
nomenon by which the plane-wave solutions to a linear 
differential equation become exponentially localized upon 
the breaking of translation symmetry by a local ran¬ 
dom potential. 1 2 Anderson localization always prevails 
for metals in one-dimensional (ID) space and for met¬ 
als without spin-orbit coupling in tw o-dimensional (2D) 
space, however small the disorder. 3 It requires the dis¬ 
order to be lar ge e nough, of the order of the Fermi en¬ 
ergy, otherwise. 3 ^ Until the experimental discovery of 
the integer quantum Hall effect (IQHE) in 1980,^ the 
most important challenge brought about by Anderson 
localization had been to understand the metal-insulator 
transition. 

The IQHE is defined by the quantized value of the 
Hall conductivity at fixed filling fraction and by a sharp 
(quantum) transition between two consecutive quantized 
values of the Hall conductivity, the plateau transition in 
short, when the filling fraction is tuned. It is explained by 
the topological character of the Hall conductivity when 
the chemical potential is in between two consecutive Lan¬ 
dau levels of a two-dimensional electron gas subjected 
to a uniform magnetic field and by the fact that this 
topological character is retained in the regime of Ander¬ 
son locahzationPEffi] The IQHE teaches two important 
lessons. First, there can be topologically distinct insu¬ 
lating phases of electronic matter. Second, direct con¬ 
tinuous (quantum) transitions between these phases are 
possible. Several approaches have been used to study 
the plateau transition from an analytical and a computa¬ 
tional point of view. Effective models such as non-linear 
sigma models (NLSMs)jD® 3 ^quantum network models,^ 
and Dirac fermion^ were proposed. In parallel, the two- 
parameter scaling theory of Khmelnitskii and Levine et 
al. was verified through numerous large-scale numerical 
simulations!- 16 In this paper, we will generalize the ap¬ 
proach pioneered by Ludwig et al. in Ref. [IS] by which 


they showed that the minimal continuum model that cap¬ 
tures the IQHE in both the clean and disordered limits is 
a Dirac Hamiltonian with random mass and gauge fields. 

Integer quantum Hall states are examples of topolog¬ 
ical insulators and superconductors defined as follows. 
Topological insulators and superconductors are states of 
matter such that fermionic quasiparticl es are gapped in 
the bulk but gapless on the boundary. 1 ' 18 There are 
many mater ials that realize topological insulators, say 
B^SeaPAlSl Topological crystalline insulators differ from 
topological insulators in that it is a spatial symmetry as 
opposed to a “generic” symmetry such as time-reversal 
symmetry that endows the boundary states with a topo¬ 
logical protection. SnTe materials are examples of topo¬ 
logical crystalline insulators supporting an even number 
of massless Dirac cones on their mirror-symmetric bound- 
ariesP^U Weyl and Dirac semimetals are yet another 
family of noninteracting fermionic Hamiltonians display¬ 
ing gap-closing points, although in the Brillouin zone 
now, that are topologically stablePsHSU 

A defining property of (strong) topological insulators 
and superconductors in d-dimensional space is that their 
boundary states are immune to Anderson localization 
provided the disorder strength is smaller than the bulk 
gap. 27 This defining property was used by Schnyder et 
al. to establish a classification of topological insulators 
and superconductors as follows For a local disor¬ 

der that is not too strong, its effect on extended bound¬ 
ary states dispersing through the gap of the bulk states 
can be treated as a problem of Anderson localization in 
(d — 1) dimensions. Following the classification by Alt- 
land and Zirnbauer of metallic or superconducting quan¬ 
tum dots (zero-dimensional space) in terms of ten uni¬ 
versality classes,P® 3 Schnyder et al. examined the pos¬ 
sibility of adding either a Wess-Zumino-Witten (WZW) 
term or a Z 2 topological term to the NLSM describ¬ 
ing (d — l)-dimensional transport along a boundary un¬ 
der the assumption that the presence of such topologi¬ 
cal terms prevents Anderson localization of the bound- 






2 


ary states. There follows a table with the ten Altland- 
Zirnbauer (AZ) symmetry classes as rows and the di¬ 
mensionality d of space as columns (see Table [T]) . The 
entries of this table, the tenfold way of topological in¬ 
sulators and superconductors, are identified as either a 
trivial insulator/superconductor when no WZW and no 
Z 2 topological term are allowed, a Z topological insula¬ 
tor/superconductor when a WZW term is allowed, or a 
Z 2 topological insulato r/sup erconductor when a Z 2 topo¬ 
logical term is allowed!^® 0 

An alternative derivation of this table was proposed 
independently by Kitaev using K theory.^ K theory en¬ 
ables a systematic study of a homotopical structure of 
gapped Hamiltonians in the bulk (upon the imposition 
of periodic boundary conditions, say) in the clean limit. 
Appropriately normalized gapped Hamiltonians define “a 
classifying space” for each AZ symmetry class. The ten¬ 
fold way of topological insulators follows from the homo- 
topy groups of the classifying spaces, with the additional 
insight that it obeys a periodicity that originates from the 
Bott periodicity in K theory. Since the periodic table is 
also obtained by using a representation theory of massive 
Dirac Hamiltonians P2 the approach relying on K theory 
also delivers a systematic classification of Dirac masses 
entering Dirac Hamiltonians satisfying the AZ symmetry 
constraints.^ 

The goal of this paper is to reverse the logic of Schny- 
der et al. Instead of classifying topological insulators 
as an exercise in evading Anderson localization on the 
boundary, we seek to understand how topology deter¬ 
mines Anderson localization in the bulk of a disordered 
insulator. To this end, we marry the approach pioneered 
by Ludwig et al. for the IQHE with that of Kitaev. In any 
dimension d and for any AZ symmetry class, we study all 
random Dirac Hamiltonians that support a Dirac mass. 
In particular, we show that the homotopy groups of the 
classifying spaces, regarded as the space of normalized 
Dirac masses in each AZ symmetry class, are essential to 
establish the phase diagram encoding Anderson localiza¬ 
tion for these random Dirac Hamiltonians. 

Our approach gives a qualitative understanding of a 
generic global phase diagram of Anderson localization for 
any AZ symmetry class in one, two, three, etc. dimen¬ 
sions. Our phase diagram is consistent with previously 
known results and gives a simple and alternative explana¬ 
tion for the even-odd effect in the one-dimensional chiral 
classes pIHSU We can also give a qualitative explanation 
for the Gade singularity® ® and Griffiths effectS^in two 
dimensions for five out of the ten AZ symmetry classes 
in terms of the first homotopy group of the classifying 
spaces. 

Our approach can also be adapted to analyze the sta¬ 
bility of massless Dirac fermions on the surface of 3D 
topological crystalline insulators. We establish the con¬ 
ditions for the stability of the boundary states of statis¬ 
tical topological crystalline insulators (STCIs) that are 
protected by a crystalline (here a reflection) symmetry 
on average for the ten AZ symmetry classes. In this way, 


we show that the alloy Sn 1 _ ir Pb x Te has surface states 
that are protected by the local symmetry class All and 
the average reflection symmetry. 

This paper is organized as follows. In Sec.[Hj we review 
the classification of topological insulators and supercon¬ 
ductors using the K theory and present classifying spaces 
of normalized Dirac masses in Dirac Hamiltonians. In 
Sec. |III[ we relate the topology of classifying spaces to 
Anderson localization. In Secs. |IV[ |V} and m we apply 
our method to ID, 2D, and 3D systems, respectively, to 
infer global phase diagrams. In Sec. m we briefly ex¬ 
plain how one may use the classifying spaces to study the 
stability of boundary states of topological insulators. In 
Sec. |VIIl[ we study the stability of the boundary states in 
a STCI. We present a summary in Sec. |IX[ In Appendix 
A we present explicit constructions of ID massive Dirac 
Hamiltonians and discuss some topological properties of 
the space of normalized Dirac masses. Finally, we discuss 
one-loop renormalization-group flows of random masses 
and potentials for 2D Dirac Hamiltonians in Appendix 
B. 


II. CLASSIFYING SPACES OF NORMALIZED 
DIRAC MASSES 


The key object of this paper is the set of allowed nor¬ 
malized Dirac masses for each dimension d of space, rank 
of the Dirac matrices, and AZ symmetry classes. This set 
can be attached a topology, i.e., it becomes a topological 
space. As a topological space, it is given by a classifying 
space associated with the extension problem of a Clif¬ 
ford algebra, which we review in this section. Classifying 
spaces are characterized by homotopy groups that obey 
Bott periodicity, a property that carries over to the al¬ 
lowed normalized Dirac masses. The homotopy groups 
characterizing the allowed normalized Dirac masses are 
used to characterize the physics of Anderson localization 
in Secs. [IHJ|!V| and[Vj 


This section is devoted to the systematic construction 
of the classifying spaces V d that are given in the sixth 
column of Table [TJ To this end, we shall proceed in the 
following steps. First, we review the algebraic definition 
of Clifford algebras. Second, we review the relation be¬ 
tween Clifford algebras and Dirac Hamiltonians. Third, 
we review the tenfold way for Clifford algebras. Fourth, 
we review the tenfold way for classifying spaces. 


A. Definitions of Clifford algebras and their 
classifying spaces 


The complex Clifford algebra 

dq = {Al? ' ' ' : G' } 


(2.1a) 
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TABLE I. The ten Altland-Zirnbauer (AZ) symmetry classes and their topological classification in terms of the zeroth 
homotopy groups of their classifying spaces. Two complex and eight real symmetry classes are characterized by the presence 
or the absence of time-reversal symmetry ( T ), particle-hole symmetry (C), and chiral symmetry (T). Their presence is 
complemented by the sign multiplying the identity in T 2 = ±1 or C 2 = ±1, and by 1 for T. Their absence is indicated by 0. 
For each symmetry class and for the dimension d of space, the relevant extension problem with its solution in terms of the 
classifying space V d = limjv-nx> V d r . n are given in the fifth and sixth columns, respectively. The zeroth homotopy groups of 
the classifying spaces for d= 0,... ,7 are given in the last eight columns. Each entry with a nontrivial zeroth homotopy group 
of the classifying space defines a noninteracting topological (superconductor) insulator. 


Class 

T 

c 

r 

Extension 

V d 

MV d - 

-o) ’ro(Vd-i) 

7T0 {V d=2 ) 

MVd= 

=3) n 0 (V d=4 ) 

MVd- 

= 5 ) TT 0 (Vd =6 ) 

7r o(i / d = 7) 

A 

0 

0 

0 

Cl d —> Cl d+ 1 

Co+d 

z 

0 

z 

0 

z 

0 

z 

0 

AIII 

0 

0 

1 

Cl d+ 1 

“^ Cl d+ 2 

C\+ d 

0 

z 

0 

z 

0 

z 

0 

z 

AI 

+1 

0 

0 

C^0,d+2 

—• h Cl 1<d+2 

Ro-d 

z 

0 

0 

0 

z 

0 

Z 2 

Z 2 

BDI 

+1 

+1 

1 

Cl d + 1,2 

Cl d+ 13 

R-l-d 

Z 2 

z 

0 

0 

0 

z 

0 

Z 2 

D 

0 

+1 

0 

ci d , 2 

^ ci dt 3 

R-2-d 

z 2 

Z 2 

z 

0 

0 

0 

z 

0 

Dill 

-1 

+1 

1 

Cl d , 3 

= > Cl dt 4 

R?,-d 

0 

z 2 

z 2 

z 

0 

0 

0 

z 

All 

-1 

0 

0 

ci 2 , d 

-*• Cl 3:d 

-^4 — d 

z 

0 

Z 2 

Z 2 

z 

0 

0 

0 

CII 

-1 

-1 

1 

Cl d + 3,0 

ci d+3} i 

R-5-d 

0 

z 

0 

Z 2 

Z 2 

z 

0 

0 

C 

0 

-1 

0 

Cl d + 2,0 

=^ Cl d+ 2 t 1 

-^6 — d 

0 

0 

z 

0 

Z 2 

z 2 

z 

0 

Cl 

+1 

-1 

1 

Cl d + 2,1 

=* Cl d + 2,2 

Rl-d 

0 

0 

0 

z 

0 

z 2 

Z 2 

z 


is a complex vector space C 29 of dimension 2 q that is Given a representation of the complex Clifford algebra 
spanned by the basis with the basis elements Cl q , “the extension problem” denoted 




A 1 j * • • 5 n q 


n =0,1, (2.1b) 


whereby the multiplication rule 

{e L ,e L ,} = 2<5 / 


(2.1c) 


applies l. i! = 1,..., q. The vector space Cl q is closed 
under multiplication of any two of its elements owing to 
Eq. (2.1c). Hence Cl q is also an associative algebra. 

The real Clifford algebra 


dp,q — { e l> • ■ • > e p l e p+ h • ■ • > e p+q} 


(2.2a) 


is a real vector space R 2P+9 of dimension 2 p+q that is 
spanned by the basis with the basis elements 


ri , 
P+<? 


p+q 

=nw 


ni,...,n p+q = 0,!, (2.2b) 


i=l 


whereby the multiplication rule 


{ e n e t'} = 2 7?,, 


p times q times 


(2.2c) 


= diag (-1,..., -1, +1,..., +1), 


applies for i,i! = 1 ,...,p + q. The vector space Cl pq 
is closed under multiplication of any two of its elements 


owing to Eq. (2.2c). Hence Cl p is also an associative 
algebra. 


Cl q —> Cl q+1 (2.3) 

consists in identifying “the classifying space” C q that 
parametrizes the representation of the generator e q+1 
present in Cl q+1 but absent in Cl q . Similarly, given a 
representation of the real Clifford algebra Cl pq , there are 
two possible extension problems. There is the extension 
problem 


Clp^q > Clp , q+ 1 (2.4) 

that consists in identifying the classifying space R q _ p 
that parametrizes the representation of the generator 
e p+9+1 present in Cl p q+l and thus satisfying e p+q+1 = 
+1, but absent in Cl p q . There is the extension problem 

Cl p+l q (2.5) 

that consists in identifying the classifying space R p _ q+2 
that parametrizes the representation of the generator 
e p+1 present in Cl p+1 and thus satisfying e 2 +1 = —1, 
but absent in Cl p q . The latter extension problem can be 
related to the former one through the homeomorphism 

Cl p q ® Cl 0 2 — Cl q p+ 2 , (2-6) 

where Cl 0 2 is a linear algebra of real two-dimensional 
matrices and does not affect the extension problem. Clas¬ 
sifying spaces depend on the difference q—p only because 
of the homeomorphism 

Clp+ l,g-)-l — ® ^1,1 — dp,q ® ^^0,2: (2-7) 
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where dropping Cl 0 2 does not affect the extension prob¬ 
lem. 

Now, K-theory makes the remarkable prediction that 
there are two families of complex classifying spaces C 0 
and C lt while there are eight families of real classifying 
spaces R 0 ,..., R 7 . In other words, the dependence on q 
enters modulo two for the complex classifying spaces, 

C g+2 ^ C q , (2.8a) 

while it enters modulo eight for the real classifying spaces, 

Rq- 1-8 — Rq ! (2.8b) 

according to Bott’s periodicity.^ All families of classify¬ 
ing spaces V in Table [TT| are labeled by the integer number 
N entering the rank 


r = 


r • N 

min 


(2.9) 


assumed for the representation of the Clifford algebras. 


B. Definition of minimal massive Dirac 
Hamiltonians 

We assume that space is d-dimensional. The kinetic 
part of a translation-invariant Dirac Hamiltonian is 

d 

?4i„(fc) = E^’ ( 2 - 10a ) 

i=1 

where k = (k 11 ..., k d ) £ is the momentum and 
a = (aq,..., a d ) are the (Hermitian) Dirac matrices that 
obey the algebra 

{a i ,a j } = 2d tJ , i,j = l,...,d. (2.10b) 

On the one hand, we assume that the dimensionality of 
the representation of the matrices a is sufficiently large 
so that there exists at least one Hermitian matrix f3 such 
that it anticommutes with all the components of a and 
it squares to the identity matrix, 

{/3,?4 in (fc)} = 0, P 2 = l. (2.10c) 

It is then possible to write the translation-invariant mas¬ 
sive Dirac Hamiltonian 


^( fc ) = ^2 k z a i +mP, 


(2.10d) 


i=1 


where m £ R is a mass. However, the matrix /3 with its 
mass to (i.e., mass matrix in short) may not be unique. 
For example, if the Dirac matrices are chosen to be of 
rank two, then there are two linearly-independent mass 
matrices anticommuting with each other in d = 1, one 
possible mass matrix in d = 2, and none in d = 3. On 
the other hand, the translation-invariant massive Dirac 


Hamiltonian (2.lOd) becomes reducible for sufficiently 


large rank of the Dirac matrices. In any of the AZ sym¬ 
metry classes, we may start from a sufficiently large ma¬ 
trix representation of the translation-invariant massive 
Dirac Hamiltonian (2. lOd), which we then reduce until 


we reach the rank of the Dirac matrices below which we 
would lose all mass matrices. In this way, one obtains for 
each AZ symmetry class and dimension d an irreducible 
translation-invariant massive Dirac Hamiltonian of the 
form \2. 1 Od ) which is of minimum rank r min (not neces¬ 
sarily unique in that more than one distinct mass matrix 
may be possible). 

To which AZ symmetry class the Dirac Hamiltonian 


« = £< 


d 

\dx-i 


■ ft m{x) 


( 2 . 11 ) 


belongs depends on whether it is possible to construct a 
combination from the triplet of operations 


T = T K, C = C K, T, (2.12a) 


for time-reversal symmetry (TRS), particle-hole symme¬ 
try (PHS), and chiral symmetry (CHS), respectively, (K 
denotes the operation of complex conjugation and T, C, 
and T are matrices sharing the same rank as the Dirac 
matrices a) such that 


= ±1, 

[t,c} = o, r 2 = i, 

(2.12b) 

TRS: 

[T,H]= 0, 

(2.13a) 

PHS: 

{^1 = 0, 

(2.13b) 

CHS: 

{T,H} = 0. 

(2.13c) 


(We have performed a global choice of gauge for which 
[T, C\ = 0 holds.) Equations (2.13) are equivalent to 


TRS: 

[T,P\={T,cx} = 0, 

(2.14a) 

PHS: 

{C,0} = [C,a]= 0, 

(2.14b) 

CHS: 

{r,/3} = {r,a} = o. 

(2.14c) 


Observe here that the antiunitarity of T and C inter¬ 
changes the action of commutators and anticommutators 
when acting on the mass relative to the kinetic part of 
the Dirac Hamiltonian (2.11). 


C. The tenfold way for the Clifford algebras 


We are ready to combine the Clifford algebras and their 
classifying spaces from Sec. |II A| with the AZ classification 
of massive Dirac Hamiltonians from Sec. Em 

We associate with each AZ symmetry class, with each 
dimension d of space, and with any rank of the Dirac 
matrices a = ( a 1 ,...,a d ) equal to or larg er than the 
minimal rank as defined below Eq. (2. lOd), a Clifford 
algebra according to the following rules. 
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TABLE II. Complex and real classifying spaces are a list of ten topological spaces that are built out of the compact Lie groups 
U(N), O(N), and Sp(N) - the unitary, orthogonal, and symplectic matrix groups, respectively - and their quotients as is given 
in the column “Classifying space”. The symbols denoting them are given in the column “Label”. The number N is related 
to the rank r = r min N of the Dirac matrices that form the Clifford algebras, i.e., N = 1,2,... is the number of copies of the 
minimal massive Dirac Hamiltonian of rank r min . Let p = 0, 1, ... index the pth homotopy group. The complex classifying 
spaces obey the periodicity condition n p (C q ) = n p+2 (C q ) for q = 0,1. The real classifying spaces obey the periodicity condition 
n p (R q ) = 7 T p+s (R q ) for q = 0,..., 7. Hence an exhaustive list of the homotopy groups of the classifying spaces is shown in 
the columns “7r 0 (l/)”, “7r 1 (V')”,...,“7r 7 (V’)”. In each homotopy column, the three entries Z hold for N larger than an integer 
(infinity included) that depends on the order of the homotopy group and the classifying space, the two entries Z 2 hold for N 
larger than an integer that also depends on the order of the homotopy group and the classifying space. The entry 0 is a short 
hand for the group {0} made of the single element 0. 


Label 

Classifying space V 

MV) MV) MV) MV) MV) MV) MV) MV) 

C 0 

U% =0 {U(N)/[U(n) X U(N — n)]} 

z 

0 

z 

0 

z 

0 

z 

0 

c, 

U(N) 

0 

z 

0 

z 

0 

z 

0 

z 

Rg 

U» =0 {O(N)/[O(n) x 0(N - n)]} 

z 

Z 2 

Z 2 

0 

z 

0 

0 

0 

Ri 

O(N) 

Z2 

z 2 

0 

z 

0 

0 

0 

z 

r 2 

0(2N)/U(N) 

Z 2 

0 

z 

0 

0 

0 

z 

Z 2 

Rs 

U (2N)/Sp(N) 

0 

z 

0 

0 

0 

z 

Z 2 

Z 2 


U£= 0 {Sp(N)/[Sp(n) x Sp(N-n)]} 

z 

0 

0 

0 

z 

z 2 

z 2 

0 

r 5 

Sp(N) 

0 

0 

0 

z 

Z 2 

Z 2 

0 

z 

Rg 

Sp(N)/U(N) 

0 

0 

z 

Z 2 

Z 2 

0 

z 

0 

Rr 

U(N)/0(N) 

0 

z 

z 2 

Z 2 

0 

z 

0 

0 


The symmetry classes A and AIII are associated with 
the complex Clifford algebra 

A: Cl d+1 = {P,a}, (2.15a) 

AIII: Cl d+2 = {P, T,c*}, (2.15b) 

respectively. 

For the remaining eight symmetry classes, it is always 
possible to define the operation J that satisfies the rela¬ 
tions 

{T, J} = {C, J} = [F, J] = [a, J] = [/?, J] = 0, (2.16) 

and plays the role of an imaginary unit for the real Clif¬ 
ford algebras as J 2 = —1. The symmetry classes AI, 
BDI, D, Dill, All, CII, C, and Cl are associated with 
the real Clifford algebras-®- 


AI: 

Ch,d+2 

= {J/3;T,7V,a}, 

(2.17a) 

BDI: 

Cld+ 1,3 

= {Ja,TCJ; C, CJ,/3 }, 

(2.17b) 

D: 

Cld,3 = 

{Ja;C,CJ,/3}, 

(2.17c) 

Dill: 

Cl dA = 

{Ja;C,C'J,TC'J,/3}, 

(2.17d) 

AIL 

Cl 3 ,d = 

{J/3,T,7V;a}, 

(2.17e) 

CII: 

Cl d + 3,1 

= {Ja,C,CJ,TCJ;/3}, 

(2.17f) 

C: 

Cl d + 2,1 

= {Ja,C,CJ;/3}, 

(2-17g) 

Cl: 

Cld+ 2,2 

= { Ja : C , CJ]TCJ,f3 }, 

(2.17h) 


D. The tenfold way for the classifying spaces V 


The definition of the classifying space V associated 


with the Dirac Hamiltonian (2.11) depends on the AZ 
symmetry class to which it belongs, the dimensional¬ 
ity d of space, and the rank r of the Dirac matrices 
a = (ay,..., a d ), which is assumed larger than the min¬ 
imal rank r min defined in Sec. IIB The classifying space 
V encodes the fact that the mass matrix j3 in Eq. ( |2. 11 ) 
might not be unique for given d and the rank of the Dirac 
matrices a = (aq,..., a d ). The construction of the clas¬ 
sifying space V proceeds in the following steps. 


Step 1: To each AZ symmetry class, we assign the 
following pair of Clifford algebras differing by one gen¬ 
erator, holding the dimensionality d of space fixed and 
the rank of the Dirac Hamiltonian fixed. The Clifford 
algebra to the left of the arrow in the fifth column “ex¬ 
tension” from Table |T] is obtained after removing from 
the t enfold list of Clifford algebras defined by Eqs. (2.15) 
and (2.17) one generator, namely the mass matrix J (3 for 
the symmetry classes AI and All and the mass matrix (3 
otherwise. This gives the tenfold list 


A: 

AIII: 


respectively. 


Cl d = {a}, 
ci d +i = {r, a} , 


(2.18a) 

(2.18b) 
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and 


AI: 

^o,d +2 — {i T, TJ , a:}, 

(2.19a) 

BDI: 

Cl d+h2 = {Jcx,TCJ-,C,CJ}, 

(2.19b) 

D: 

Cl d ' 2 = {Jcc,C,CJ}, 

(2.19c) 

Dill: 

Cl d ,3 = {Joc;C, CJ,TCJ}, 

(2.19d) 

All: 

C/ 2id = {T,TJ;a}, 

(2.19e) 

CII: 

Cl d + 3,0 = { Jcv.,C,CJ,TCJ ;}, 

(2.19f) 

C: 

Cld+ 2,0 = {da, C, CJ; }, 

(2.19g) 

Cl: 

Cl d+2tl = {Ja,C,CJ-,TCJ}. 

(2.19h) 


Step 2 : For each AZ symmetry class, we seek all dis¬ 
tinct Hermitian matrices sharing the same rank as the 
Dirac matrices a = (aq,..., a d ) such that they can 
be added to the list of generators entering the corre¬ 
sponding Clifford algebra from the tenfold list defined 
by Eqs. (2.181 and (2.19) so as to deliver the corre¬ 


sponding Clifford a lgebra from the tenfold list defined 
by Eqs. (2.151 and (2.17). These mass matrices form a 


set V, the classifying space. Determining V is an exam¬ 
ple of the extension problem in K-theory. A character¬ 
ization of the classifying space V can be deduced from 
K-theoryJ 34 l 45 l The outcome of this exercise is listed in 
the second column of Table |n| (owing to the Bott period¬ 
icity) for the case when the rank r of the Dirac matrices 
a = (aq,..., a d ) is 


r = 




( 2 . 20 ) 


where r min is the rank of the minimal representation for 
the Dirac Hamiltonian (2.11) and TV is an integer. Ex¬ 
plicit examples for the construction of V can be found in 

Ref. [35] 

The zeroth homotopy group of the classifying spaces 
is given in the third column of Table [II] This homotopy 
group has the following physical consequences, as shown 
by Kitaev in Ref. |3H Imagine that d-dimensional space 
is divided into two halves. Both halves share a (d — 1)- 
dimensional boundary. Whenever the zeroth homotopy 
group of the classifying space in the tenfold way is non¬ 
vanishing, consider the Dirac Hamiltonian (2.11) with 
the mass term interpolating smoothly across the (d — 1)- 
dimensional boundary between two fixed elements in the 
classifying space characterized by distinct values of the 
zeroth homotopy group. By the very definition of a ho¬ 
motopy group, this is only possible if the mass term van¬ 
ishes along the (d — 1)-dimensional boundary separating 
the two halves of d-dimensional space. As was shown 
by Jackiw and Rebbi for the sym metry class BDI when 
d = lpSl the Dirac Hamiltonian ( 2.1 1| ) must then sup¬ 
port a zero mode that is extended along the (d — 1)- 
dimensional boundary but exponentially localized away 
from it. This is the defining property of a topological 
insulator (superconductor). Hence, by combining the ze¬ 
roth homotopy group of the classifying spaces given in 
the third column of Table |TT] with the Bott periodic¬ 
ity (2.8), one infers which of the AZ symmetry classes 


allows a topological (superconductor) insulator for any 
given dimensionality d of space. The periodic table for 
topological (superconductors) insulators follows 122IMI 

Another application of the zeroth homotopy groups as¬ 
sociated with the classifying spaces is relegated to the 
Appendix [Aj where we show in which sense one may de¬ 
clare that a mass matrix is “unique”. 

All higher homotopy groups of the classifying spaces 
are given in column four to ten of Table [TT] for they obey 
the periodicities 

TTpiCj) = n p+ 2 (Cj), p = 0,1,2,..., (2.21a) 

for the complex classes j = 0,1 and 

n p (Rj) = TT p+ 8 {R j ), p = 0,1,2,..., (2.21b) 

for the real classes j = 0,1,..., 7. They also obey the 
translation rules 

n P ( C j) = 7r P + 1 ( C 'i+ 1 ), P = 0,l,2,..., (2.22a) 

for the complex classes with the integers j and j + 1 
defined modulo 2 and 

^ P { R j) = n p+i( R j-i)’ P = 0,1,2,..., (2.22b) 

for the real classes with the integers j and j — 1 defined 
modulo 8. If we combine these translations rules with 
the definition of V d given in Table [T| we find the relation 

* D (V d ) = n 0 (V d _ D ) (2.23) 

for any D = 0,1, 2,..., d. 

E. Existence and uniqueness of normalized Dirac 
masses 

We summarize the properties of d-dimensional Dirac 
Hamiltonians that follow from the topology of the classi¬ 
fying spaces V d for the ten AZ symmetry classes. 

1. For each symmetry class, there exists a minimum 
rank r min (d), an even integer equal to or larger than 
the integer 2, for which the d-dimensional Dirac 
Hamiltonian supports a mass matrix and below 
which either no mass matrix or no Dirac kinetic 
contribution are allowed by symmetry. 

2. Suppose that the rank of the d-dimensional mas¬ 
sive Dirac Hamiltonian is r = r min (d) TV. When 
there exists a mass matrix squaring to unity (up 
to a sign) that commutes with all other symmetry- 
allowed mass matrices, we call it the unique mass 
matrix. There are the following three cases depend¬ 
ing on the entries in the 7r 0 (V d ) column of Table]!] 

(a) 7r 0 (V d ) = Z: there is always a unique mass 
matrix for TV > 1. 
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(a) *o(V d , r ) = Z 



FIG. 1. Connectedness of the compact topological space 
V dr parametrized by the normalized Dirac masses. The ze¬ 
roth homotopy groups i^o{V dr ) index the disconnected parts 
of the compact topological space V d r . The zeroth homotopy 
group n 0 (V d r N ) is either (a) Z (in the limit N — > oo), 
(b) Z 2 , or (c) {0} = 0 according to Table [IT] (a) When 
7r 0 (limjv-t.oo V d r N ) = Z, V dr N is disconnected and 
given by the union of path-connected and compact topologi¬ 
cal subspaces indexed by the half-integers or integers defined 
in Eq. (3.8c l. For odd N, V d r N is the union of path- 
connected and compact topological subspaces labeled by neg¬ 
ative and positive half integers. The total “volume” of the 
union of all subspaces with negative half-integer labels equals 
that of the union of all subspaces with positive half-integer 
labels when the mean values for the Dirac masses vanish. For 
even N, V d r N is the union of path-connected and compact 
topological subspaces labeled by the integers (3.8c I. The total 
“volume” of the union of all subspaces with strictly negative 
integer labels equals that of the union of all subspaces with 
strictly positive integer labels when the mean values for the 
Dirac masses vanish. The subspace labeled by the integer 0 
has a distinct (larger) “volume” when the mean values for 
the Dirac masses vanish, (b) When TTo{V d r ) = Z 2 , V d<r is the 
union of two path-connected and compact topological sub¬ 
spaces that are indexed by the integers ( 3. 15d[ > and are of 
equal “volume” when the mean values for the Dirac masses 
vanish, (c) When i^o[y d r ) = {0} = 0, V d r is a path-connected 
and compact topological space. 


(b) 7r 0 (H d ) = Z 2 : there is a unique mass matrix 
/3 min only when N = 1. When N is an even 
integer, for any given mass matrix, there ex¬ 
ists another matrix that anticommutes with 
it. When N is an odd integer larger than 1, 
the matrix /3 min (g> plays a role similar to 
that of the unique mass matrix in that the two 
matrices ±/3 min <8> Ijy belong to different con¬ 
nected components of V (j . even though there 
exists a normalized mass matrix that neither 



FIG. 2. For each realization of a random potential perturb¬ 
ing the massless Dirac Hamiltonian defined in d-dimensional 
Euclidean space R d , we may decompose R“ into open sets 
(domains) of linear size £ dis . In each of these domains, the 
normalized Dirac masses correspond to a unique value of 
their zeroth homotopy group. At the boundary between do¬ 
mains differing by the values taken by their zeroth homotopy 
group, the Dirac masses must vanish. Such boundaries sup¬ 
port quasi-zero-energy boundary states. When it is possible 
to classify the elements of the zeroth homotopy group by the 
pair of indices A and B, we may assign the letters A or B to 
any one of these domains as is illustrated. When the typi¬ 
cal “volume” of a domain of type A equals that of type B, 
quasi-zero modes undergo quantum percolation through the 
sample and thus establish either a critical or a metallic phase 
of quantum matter in d-dimensional space. 


commutes nor anticommutes with /3 min ® Ijy 

(c) 7T 0 (V d ) = 0: for any given mass matrix, there 
exists another mass matrix that anticommutes 
with it. That is, there is no unique mass ma¬ 
trix for N > 1. 

When (i) N = 1, (ii) the Dirac Hamiltonian has a 
unique mass matrix that realizes topologically dis¬ 
tinct ground states for different signs of its mass, 
and (iii) the mass is varied smoothly in space, then 
domain boundaries along which the mass vanishes 
are accompanied by massless Dirac fermions, whose 
low-energy Hamiltonian is of rank r min (d)/2. It fol¬ 
lows that r min (d - 1) = r min (d). 

We illustrate in Appendix [A] these properties for ID 
Dirac Hamiltonians in the ten AZ symmetry classes. 


F. Relationship to higher homotopy groups 


The zeroth homotopy group of a topological space in¬ 
dicates if it is path connected and, if not, how to in¬ 
dex all its disti nct subspaces that are path connected. 
Equation (2.23) relates the zeroth to the higher homo¬ 
topy groups of the classifying spaces V d defined in Table 










[IJ Equation (2.23) can be given the following interpreta¬ 
tion. 

We recall that the homotopy group n n (X) is the set 
of homotopy classes of maps / : S n -A X between the 
unit sphere S n in (n + l)-dimensional Euclidean space 
and the topological space X. Following Teo and Kane in 
Ref. EH we identify D in Eq. ( |2.23[ ) as the dimensionality 
of the sphere S D that surrounds a defect in d-dimensional 
space. For point, line, and surface defects, D = d — 1, 
D = d — 2, and D = d — 3, respectively. In other words, 
D +1 is the codimension of such defects in d-dimensional 
space. A homotopy group TT D (V d ) with more than one 
element signals that defects of codimension D + 1 in the 
normalized Dirac masses can be indexed by a topological 
number. For given D and d , Eq. (2.231 [both d and d — D 
on the left- and right-hand sides of n D (V d ) = n 0 (V d _ D ), 
respectively, are defined either modulo 2 or modulo 8 de¬ 
pending on the AZ symmetry class] dictates through the 
homotopy reduction D ^ D—D = 0 and the dimensional 
reduction d -A d — D which five of the ten AZ symmetry 
classes support topological defects in their normalized 
Dirac masses. In particular, point defects inherit the 
topological numbers from the zeroth homotopy group of 
V 1 . Any Dirac Hamiltonian with a defective normalized 
Dirac mass of topological character supports eigenstates 
with a vanishing energy eigenvalue (zero modes) that are 
bound to the defect in the directions transverse to it. 
These zero modes are robust to any local perturbation 
as long as they respect the AZ symmetry class and they 
are not too strong. 


III. ANDERSON LOCALIZATION AND THE 
ZEROTH HOMOTOPY GROUP OF THE 
CLASSIFYING SPACES 


We assume space to be d-dimensional. Furthermore, 
we assume that a microscopic lattice model with lattice 
spacing a describing noninteracting fermions propagating 
in a static random environment respecting one of the ten 
symmetry constraints from the AZ symmetry classes is 
captured by the Dirac Hamiltonian 

d r\ 

= - bV(®) + ..., (3.1a) 

»=l 1 * 

in the low-energy and long-wavelength limit. The rank 
of the Dirac matrices is 


N, 


(3.1b) 


where r min is the smallest rank that admits a Dirac mass 
matrix and N = 1,2,.... Hence, for any x £ R d , there 
exists a Dirac mass matrix V(x) of rank r that competes 
with the kinetic contribution parametrized by the Dirac 
matrices a = (aq,...a d ) [i.e., V(x) anticommutes with 
a]. Th e mat rix elements of the Dirac mass matrix V ( x ) 
in Eq. (3.1a) are random functions of x £ R d . The dots 


represent all other static random vector and scalar po¬ 
tentials allowed by the AZ symmetry class. We fix the 
chemical potential /i to be /i = 0. 

The matrix elements of the Dirac mass matrix V(x) are 
assumed to be random functions that change smoothly 
in space on the length scale of £ dis (» a). Their corre¬ 
lations are assumed local in that these matrix elements 
that are not related by the AZ symmetries are uncor¬ 
related up to an exponential precision beyond the finite 
length scale £ dis a, e.g., (disorder averaging is denoted 
by an overline) 


V(x) =: m /3 0 , (3.1c) 


and 


~tr[V (x) — m /3 0 ] [V(y) — m /3 0 ] =: g 2 e |x 

(3. Id) 


with all higher cumulants vanishing. The choice of the 
normalized Dirac mass matrix /? 0 from the classifying 
space V d r will be done in Secs. 


Ill A 


IIIB 


and IIIC 


such a way that the parameter space (m, g) £ R. x [0, oo[ 
captures the phase diagram representing the competi¬ 
tion between delocalized and all topologically distinct 
localized phases of noninteracting fermions in a given d- 
dimensional AZ symmetry class. The former phase is fa¬ 
vored by the Dirac kinetic contribution. The latter phase 
is favored by the Dirac masses. In other words, localized 
(insulating) phases are favored by large |m|, whereas in¬ 
creasing g generates more density of states in the band 
(mass) gap and, in doing so, helps delocalization. 

The main result of this section is captured by Fig. [T] 
a consequence of Table [Tj and Fig. [2j Let V d r be the 
classifying space associated with any typical realization 
of the random Dirac Hamiltonian & We then have 
the following phases (at the band center e = 0 of the 
quasiparticle dispersion): 


(a) If limjv->-oo 7To (V d r N ) = Z, there are N + 1 topo¬ 
logically distinct insulating phases that are sepa¬ 
rated pairwise by either a critical point or a metal¬ 
lic phase. 

(b) If T^o{V dr N ) = Z 2 , there are two topologically 
distinct insulating phases that are separated by ei¬ 
ther a critical point or a metallic phase. 

(c) If 7 T Q (V dr N ) = 0, there is only a topologically 

trivial insulating phase. 


When d > 2 (where the equality holds for class All, CII, 
D, and Dill), the ground state is metallic for sufficiently 
large g. 

The key intuition to support items (a) and (b) is the 
following. Since the spatial variation of V{x) is as¬ 
sumed to be smooth for any realization of the disorder, 
d-dimensional space can be decomposed into open sets 
(domains) with the characteristic size £ dis such that 
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(1) for each domain the values taken in it by V(x) 
can be assigned an index from the zeroth homo- 
topy group 7 T 0 (V d r ), 

(2) and det[V(a;)] = 0 along the boundary of each do¬ 
main. 

Figure [2] is an illustration of this decomposition of d- 
dimensional space when either 7r 0 (V^ r ) = Z or n 0 {V d r ) = 
Z 2 . There are then gapless modes bound to the bound¬ 
aries of these domains. We use a semiclassical picture in 
analogy with the Chalker-Coddington network model of 
the IQHE. When no connected boundaries defined by the 
condition det[H(a:)] = 0 percolate across d-dimensional 
space, we expect an insulating phase. However, if a con¬ 
nected boundary along which det[V(a:)] = 0 percolates 
across d-dimensional space, we expect departure from 
an insulating phase. Quantum mechanics modifies this 
percolating picture by turning it into that of a Chalker- 
Coddington-like (quantum) network model in dimension 
d, whereby the scattering matrix at each node of the net¬ 
work is fixed by the AZ symmetry class and r = r min N. 


A. Case of the zeroth homotopy group Z 


In each dimension d of space, there are three AZ sym¬ 
metry classes whose classifying spaces V dr N are the 
unions of one of the three Grassmannian manifolds, i.e., 
the classifying spaces are any one of 

U {u(N)/[U(n)xU(N-n)}}, (3.2a) 

n= 0, " ,N 

U {0(N)/[0(n)x0(N-n)}} : (3.2b) 

n— 0,-" ,N 

U {SpWPp(n) X Sp(N - n)]}. (3.2c) 

n= 0,-" ,N 


These unions of Grassmannian manifolds are realized by 
the topological space of normalized Dirac masses when¬ 
ever there exists a unique (up to a sign) normalized mass 
matrix that commutes with all other allowed normalized 
mass matrices. In the limit N —► oo, the zeroth homo¬ 
topy group of either one of these Grassmannian manifolds 
is Z. 

When the rank (3.1b I of the Dirac Hamiltonian ( |3.1a ) 
belonging to any one of these AZ symmetry classes is the 
minimal one, iV = 1 and there is a unique (up to a sign) 
normalized Dirac mass matrix /3 0 of rank r min such that 


V ( x ) = m{x) /3 0 


(3.3) 


in Eq. (3.la I. The Dirac Hamiltonian (3.1a I with the 


uniform mass m > 0 is topologically distinct from the 
Dirac Hamiltonian (3.1a) with the uniform mass m < 0. 


Correspondingly, the classifying spaces (3.2) reduce to 


U U(l)/[U(n)xU(l-n)]*{-l,+l}, 

n= 0,1 

U 0(l)/[0(n)x 0(1-n)]-{-!,+!}, 


n=0,l 


n=0,l 


'I; +1}) 

(3.4a) 

‘I; +1}) 

(3.4b) 

{ — 1j +1}' 

(3.4c) 


If the Dirac mass matrix V(x) of minimal rank is ran¬ 
dom with the statistical correlation (3.1c) and (3.Id), we 


may decompose d-dimensional space in disjoint domains 
as depicted in Fig. [2] A typical domain has the linear size 
£ dis . In it m{x) ^ 0 with a given sign, along its bound¬ 
ary m(x) = 0, and a sign change is only permissible 
across this boundary into another domain with opposite 
and constant sign of the mass m(x) ^ 0. Any boundary 
separating two domains with opposite signs of m(x) ^ 0 
binds gapless boundary states. Whenever two boundaries 
approach each other within a distance much smaller than 
£dis> boundary states undergo an elastic quantum scat¬ 
tering process dictated by the AZ symmetry class. In 
other words, Fig. [2] defines a quantum network model, 
whereby incoming plane waves along the boundaries de¬ 
fined by the condition m(x) = 0 scatter off each other 
elastically at the nodes of this network of boundaries. 
The mean value m of the random Dirac mass m(x) dic¬ 
tates the relative volume occupied by the domains with 
sgn[m(a:)] = +1 relative to the volume occupied by the 
domains with sgn[77i(a:)] = —1. When m = 0, both vol¬ 
ume are typically equal, in which case the domain bound¬ 
aries percolate across the system and the boundary states 
are delocalized and signal either a critical or a metallic 
phase. 

The situation is different when N = 2. Indeed, the 


unions of Grassmannian manifolds (3.2) are now com¬ 
prised of the pair of Grassmannians 


U(2)/[U(2) x 17(0)], 
0(2)/[0(2) x 0(0)], 
Sp(2)/[Sp(2) x Sp( 0)], 


U(2)/[U(0) x 17(2)], (3.5a) 
O(2)/[O(0) x 0(2)], (3.5b) 
Sp(2)/[Sp(0) x Sp( 2)], (3.5c) 


with the dimensions 0, 0, and 0, respectively, and the 
Grassmannians 


U(2)/[U(1) x 17(1)], 

(3.6a) 

0(2)/[0(l) x 0(1)], 

(3.6b) 

Sp{2)/[Spil) x Sp(l)], 

(3.6c) 


with the dimensions 2, 1, and 4, respectively. Accord¬ 
ingly, there are three topologically distinct insulating 
phases when N = 2. To derive the dimensions of these 
Grassmannian manifolds, we used the fact that U(n), 
0(n), and Spin) have the dimensions n 2 , n[n — l)/2, 
and n{2n + 1), respectively. 

Imagine that d-dimensional space is randomly decom¬ 
posed into three types of domains according to the rule 
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that the random Dirac mass matrix V(x) is associated 
with only one of the three path-connected Grassmannian 
manifolds making up the classifying space (3.6) in each 
domain. The domain boundaries bind gapless states, per¬ 
colation of which leads to delocalization or criticality be¬ 
tween localized phases as in the case of N = 1. The 
transitions can be induced by changing the parameter 
m. Unlike the N = 1 case, however, the ground state 
at m = 0 and for any non-vanishing g not too strong is 
most likely a localized phase. Indeed, localization is most 
likely to occur because d-dimensional space is randomly 
partitioned into domains such that most of the domains 
are characterized by a random Dirac mass matrix V ( x ) 
associated with the path-connected Grassmannian man¬ 
ifolds of largest dimension in Eq. (3.6). 

This difference between the N = 1 and N = 2 cases 
is not accidental. The same difference holds between the 
cases of odd N and even N integers, namely that the tun¬ 
ing m = 0 delivers typically a critical or metallic phase of 
quantum matter in d-dimensional space when N is odd, 
while it delivers typically a localized phase of quantum 
matter in d-dimensional space when N is even, as we now 
explain. 

When the classifying space V d r N is any one of the 
three unions of Grassmannian manifolds (3.2), we can 
always choose to represent the Dirac Hamiltonian (3.1a) 
with 


ol = a min ® 1 j\r (3.7a) 

for the Dirac kinetic contribution and 

V(x) = fi min ®M(x) (3.7b) 

for the Dirac mass contribution, whereby ct min and /3 min 
represent the Clifford algebra with rank r min , while is 
a unit NxN matrix and M(x) is a random NxN Hermi- 
tian matrix which is a smooth continuous function of x. 
For a given realization of the random matrix M(x), we 
may partition d-dimensional space into domains whose 
boundaries are defined by det [M(a:)] = 0. Each domain 
can be assigned a topological index as follows. We may 
index each element of 7r 0 (Vb N ) with the values of v 
defined by 


1 {-1/2, +1/2}, if A 7 = 1, 

v : = - tr{sgn[M(a;)]} G < {-1,0,+1}, if N = 2, 

I and so on, if N > 2, 

(3.8a) 

where the Hermitian matrix M(x) is diagonalized by the 
unitary matrix U(x), 


M =: t/ 1 ’ diag (A 1; ..., X N ) U, (3.8b) 

and 

sgn[M(*)] := U t diag ..., |pp j U. (3.8c) 


Each domain with V(x) a smooth function of x has 
thereby been assigned the topological index v. 

We still need to choose th e para mete r spac e (m,g) G 
ffi. x [0, oo[ announced in Eqs. (3.1c) and (3.Id), when the 
zeroth homotopy group of the classifying space is Z. We 
must distinguish two cases. 

Case when g > 0 is a relevant perturbation to 
the clean critical point: We select the normalized 
Dirac mass matrix 


Po : = Pmin ® ( 3 ' 9a ) 

that anticommutes with all the components of a and 
commutes with all Dirac mass matrices allowed for given 
r and symmetry constraints (see also Appendix A 3 a I. 
The matrix (3 0 is the unique mass matrix introduced in 
Sec. HE We define the parameter space (m,g) G M x 
[0, oo [ through the probability distribution of the Dirac 
mass matrix V(x) given by 


V(x) =: m/3 0 , 


(3.9b) 


and 


—■ o - 2 P -\ x ~y\/£dit 


-tv{[V(x) - m (3 0 \[V(y) - m/? 0 ]} =: g“ e 

(3.9c) 

with all higher cumulants vanishing. With the definitions 
(3.9b) and (3.9c) for the probability distribution of V(x), 
the point m = g = 0 is a massless Dirac critical point. 
The critical point m = g = 0 separates two insulating 
phases with v = ±N/2 along the horizontal axis g = 0 in 
parameter space. By assumption, g > 0 is relevant in the 
vicinity of the clean critical point at m = g = 0. This is 
the rule when d = 1, in which case there appears for g > 0 
N — 1 additional localized phases, with v = —(N/2) + 
1, —(N/2) + 2,..., (N/2) — 1 that are pairwise separated 
by lines of critical points, all of which emerge from the 
critical point m = g = 0, as was demonstrated in Ref. 1551 
(It is the dimerization denoted by / in Ref. [55] that plays 
the role of m.) For those symmetry classes in d = 2 
for which g > 0 is (marginally) relevant in the vicinity 
of the clean critical point at m = g = 0, we conjecture 
that N — 1 additional localized phases that we may again 
label with i/ = —(N/2) + l,-(N/2) + 2,..., (N/2) - 1 
are stabilized when g > 0. Consecutive localized phases 
are either separated by a line of critical points or by a 
metallic phase. The relevant symmetry classes are AIII, 
BDI, and CII in d = 1 and A and C in d = 2. 

Case when g > 0 is an irrelevant perturbation to 
a clean critical point: We define the parameter space 
(m,g) G R x [0, oo[ through the probability distribution 
of the Dirac mass matrix V(x) given by 


V(x) =: m/3 0 + V 0 , (3.10a) 

where /3 0 again commutes with all other mass matrices 
permitted by the symmetry class, and 

“ tr |[U(x) -V(x)][V(y) - V(y)]j =: g 2 e ~ lx ~ vl/UB . 

(3.10b) 
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Here, V 0 is any mass matrix permitted by the symmetry 
that satisfies the condition 


^0 ~~ /^min ® M 0 , 


(3.10c) 


where the N x N Hermitian matrix M 0 has N nondegen¬ 
erate eigenvalues. The existence of the r x r Hermitian 
matrix V 0 is required to obtain N -f 1 distinct localized 
phases in the phase diagram by changing the parameter 
m in the clean limit g = 0. The prescription (3.101 ap¬ 


plies to the symmetry class D in d = 2 and all symmetry 
classes with 7r 0 (V) = Z in d > 3. 

More discussions on the relationship between the RG 
flow of g and the probability distribution of the Dirac 
mass matrix V(x) are presented in Sec. VIB for 3D dis¬ 
ordered systems and in Appendix |B 3 c| for the 2D systems 
of the symmetry class D. 

The zeroth homotopy group of the topological space 


V 


d,r 


N encodes the connectedness of V d r 


n- 


The 


“volume” of each path-connected component 

U{N)/[U{n) x U{N-n)], 
0(N)/[0(ri) x O(N-n)], 
Sp(N)/[Sp(n) x Sp(N - n )], 

of V d n is measured by the dimension 

2 n{N -n) = N 2 - n 2 - (TV - n) 2 , 

N(N - 1) - n(n - 1) 


i(N — n) = 


(TV — n)(N — n — 1) 


(3.11a) 

(3.11b) 

(3.11c) 


(3.12a) 


(3.12b) 


4 n{N -n) = N(2N + 1) - (TV - n)[2(N - n) + 1] 

— n{2n + 1), (3.12c) 

respectively, as is depicted schematically in Fig. ER 
When TV is odd, one can always write 


^ d,r min N — A U B, 


(3.13a) 


where 


A := (J {U{N)/[U(N-n)xU(n)]} (3.13b) 


n= 0,...,- 


and 


B := |J {U(N)/[U(n)xU(N-n)]}, (3.13c) 

n=^,...,N 


and similarly with the substitutions U —> O, Sp. The ex¬ 
istence of a critical or metallic phase of quantum matter 


when m = 0 in the parametrization (3.91 follows from 


repeating the argumentation for the N = 1 case that is 
captured by Fig. [2} 

When N is even, one can always write 


V r 


d ’ r min jv - d F 0 U V + , 


(3.14a) 


where 

V_ := (J {U(N)/[U{N-n) x U{n)}}, (3.14b) 

V 0 := U(N)/[U(N/ 2) x U(N/ 2)], (3.14c) 

and 

V + := |J {U(N)/[U(n) x U(N - n)]}, (3.14d) 

n=f +1,...,JV 


and similarly with the substitutions U —> O, Sp. Since 
the index v assigned to localized phases is an odd function 
of m, the ground state at m = 0 and g non-vanishing but 
not too large in the parametrization (3.91 is expected 
to be in the localized phase, for d-dimensional space is 
partitioned into domains of Dirac mass matrices which 
are predominantly drawn from V 0 . 


B. Case of the zeroth homotopy group Z 2 


In each dimension d of space, there are two AZ symme¬ 
try classes whose classifying spaces V dr N are home- 
omorphic to either the orthogonal group 0(N ) or the 
quotient space 0(2N)/U(N) and thus have the zeroth 
homotopy group Z 2 . The former case is called the first 
descendant Z 2 . The second case is called the second de¬ 
scendant Z 2 . If so, we can always choose to represent the 
Dirac Hamiltonian (3.1a) with 


a = a min ( 8 ) Ijv 


(3.15a) 


for the Dirac kinetic contribution and 


V{x) = p min ®M(x) 


(3.15b) 


for the Dirac mass contribution in each domain where 
the Dirac mass matrix V(x) is continuous and invertible. 
Here, /9 min is a r min /2 x r min /2 matrix such that, when 
tensored with the antisymmetric Pauli matrix <r 2 , o; min 
and p min ® <j 2 deliver a representation of the Clifford al¬ 
gebra of rank r min . Finally, Ijy is a unit N x N matrix 
and M(x) is a 2 N x 2 N Hermitian matrix which is also 
antisymmetric, i.e., 

M{x) = M\x), M(x) = —M t ( x). (3.15c) 

Each domain in the partition of d-dimensional 
space into domains defined by the boundaries where 
det[M(a;)] = 0 can be assigned a topological index as fol¬ 
lows. We may index each element of 7r 0 (V dr N ) with 
the values of v = 0,1 defined by 


(-i r 


Pf [i M(x)\ 
•\/det [i M(x)] 


(3.15d) 


This topological number v is well-defined because i M(x) 
is a real-valued antisymmetric matrix and det [iM(x)\ is 
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positive. Each domain has thereby been assigned the 
Z 2 -valued topological index v. 

When N is odd, we choose the parameter space 
(m,g)elx [0, oo [ by selecting 

A) : = Anin ® ® liv (3.16a) 

in the probability distribution of the Dirac mass matrix 
V(x) given by 


V(x)=:m /? 0 , (3.16b) 


and 

“ tr {[y (x) - m f3 0 ][V(y) - m/? 0 ]} =: g 2 e ~ lx ~ vl/idi % 

(3.16c) 


with all higher cumulants vanishing. In the clean limit 
g = 0, the point m = 0 is a massless Dirac critical point 
separating the two insulating phases with v = 0 (m > 0) 
and v = 1 (m < 0). Given g > 0, the d-dimensional space 
is decomposed into domains, and we may identify the do¬ 
mains labeled by A in Fig. [2] with v = 0 and the domains 
labeled by B in Fig. [2] with v = 1 for any realization of 
the random Dirac mass matrix V(x). Hence tuning the 
mean value m to some critical value (e.g., m = 0) real¬ 
izes a situation where the domains A and B appear with 
equal probability and the domain boundaries percolate. 
This tuning stabilizes either a critical point or a metallic 
phase in d-dimensional space. 

When N is e ven, we do not adopt the probability dis¬ 
tribution (3.16), for it leads to the massless Dirac point 


m = g = 0 separating two insulating phases belonging to 
the same topological phase with v = 0 in the clean limit 
g = 0. For example, when we consider N — 2 and the 
first descendant Z 2 (see Appendix |A 3^ .), we may choose 

M(x) =m 2 0 (x) a 2 (8> t 0 + m 2 A 33 ) ct 2 <8> 

+ m 23 (x) a 2 <g> t 3 + m 12 (x) v 1 <g> r 2 , (3.17) 


where the quadruplets a^ and r„, /r = 0,1, 2,3, are made 
of the unit 2x2 matrix and the three Pauli matrices, 
respectively. In this case 


(- 1 )" = sgn(m 2 0 + mj 2 - m 21 - m 23 ). (3.18) 

If we define the controlling parameter m as in Eq. ( |3T6l ), 
we have m 2 0 = m and 2 = m 2 x = m 2 3 = 0, for 
which the localized phase with v = 0 always appear for 
non-vanishing m with sufficiently small g. Instead, we 
choose a probability distribution such that the massless 
Dirac point m = g = 0 in the parameter space (m,g) £ 
R x [0, oo[ separates two insulating phases with v = 0 and 
v = 1 in the clean limit g = 0. For the above example of 
N = 2, such a choice is given by 


V{x) = 0, (3.19a) 


m 2,o + TO 1,2 - m 2,i - m 2,3 = : m > (3.19b) 


and 

1 tr [V(x) V{y)} =: g 2 e -l*»-wl/Cai.. 


(3.19c) 


For general values of N, we may adopt m 
instead of Eq. (|3.16h . 


Pf [i M(x)\ 


C. Case of the zeroth homotopy group {0} 


In each dimension d of space, there are five AZ symme¬ 
try classes whose classifying spaces V d r are compact and 
path-connected topological spaces and have thus a van¬ 
ishing zeroth homotopy groups. We can always choose 
to represent the Dirac Hamiltonian (3.1a) with 


1 


N 


(3.20) 


for the Dirac kinetic contribution, and then we choose 
/3 min arbitrarily from the allowed r min x r min normalized 
Dirac mass matrices which anticommutes with all the 
components of a m in- 

For any of these five AZ symmetry classes, we choose 
the parameter space (m, g) £ R x [0, oo[ by selecting 

A) := Anin ® IjV (3.21a) 

in the probability distribution of the Dirac mass matrix 
V(x) given by 


V(x) =: m/3 0 , (3.21b) 

“ tr {[V’(aj) -m/3 0 ][V(y) - m^ 0 ]} =: g 2 

(3.21c) 


with all higher cumulants vanishing. 

In these five AZ symmetry classes, the phase diagram 
has only a single localized phase that is adiabatically con¬ 
nected to a topologically trivial band insulator with re¬ 
ducing disorder. 

Nevertheless, there can be anomalies of the conductiv¬ 
ity and of the density of states at e = 0 when higher than 
the zeroth homotopy group of V d r have more than one 
elements, as we shall illustrate in Sec. [V] 


IV. APPLICATION TO ID SPACE 

The effects of static and local disorder are always 
strong in ID space. The ballistic transport for a finite 
number N of ID channels is unstable to disorder; An¬ 
derson localization rules. However, there are exceptions 
to the rule, i.e., the localization length and the density 
of states are divergent at the boundaries of insulating 
phases in the phase diagrams of the symmetry classes 
AIII, BDI, CII, D, and DHL The main result of this 
section is that, according to Sec. m these anomalies are 
caused by the nature of the disconnectedness of the topo¬ 
logical space parametrized by the normalized Dirac mass 
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TABLE III. Interplay between the topology of the normalized Dirac masses in a given AZ symmetry class and Anderson 
localization when space is one-dimensional. For each AZ symmetry class from the first column, (i) the second column gives the 
minimal rank r min of the Dirac matrices for which a Dirac mass is allowed by the symmetries, (ii) the third column gives the 
corresponding topological space for the normalized Dirac masses V d=1 r=r N , (iii) the fourth column gives the corresponding 
zeroth homotopy group n 0 {V d=l r=r N ), and (iv) the fifth column gives the corresponding panel from Fig. [3| Each panel 
from Fig. 131 is a phase diagram for quasiparticles at the Fermi energy, which is fixed to e = 0 in the cases of tne three chiral 
classes (AIII, BDI, and CII) and the four BdG classes (C, Cl, D, and Dill). Parameter space is two-dimensional with the mean 
value m £ R of the characteristic Dirac masses as horizontal axis and the characteristic disorder strength g > 0 as vertical 
axis. The last column characterizes transport through the dependence of the mean conductance as a function of the length L 
of the disordered region when the random masses are identically distributed with a vanishing mean value m = 0. Insulating 
and critical phases are characterized by a mean conductance that decays exponentially and algebraically fast with L/(NC), 
respectively, where £ is the mean free path. The entries “even-odd” indicate that the conductance decays algebraically for odd 
N and exponentially fast for even N. In each of the zeroth homotopy column, the three entries Z hold in the limit N —» oo, 
while the entry 0 is a short hand for the group {0} made of the single element 0. 


AZ symmetry class 
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Zj 
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critical 
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z 
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even-odd 
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4 
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0 
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insulating 

Cl 

4 

Ro 

0 

(c) 

insulating 


as encoded by the zeroth homotopy group of the relevant 
classifying space in the third column of Table m 

We can apply the lessons from Sec. m to deduce the 
qualitative phase diagram for disordered ID wires in any 
symmetry class. The phase diagram is spanned by the 
control parameters m and g defined in Sec. m as shown 
in Fig. [3) 

We note that, regarding transport of quasiparticles of 
any non-vanishing energy e, there are only three relevant 
symmetry classes, the standard symmetry classes A, AI, 
and AIL The chiral symmetry classes AIII, BDI, and CII 
and the BdG symmetry classes C, Cl, D, and Dill show 
their characteristic transport properties at the band cen¬ 
ter e = 0 of their quasiparticle spectra. In other words, 
a non-vanishing e induces a crossover from one chiral or 
BdG class to one standard symmetry classP^ 

Before presenting our results, we review the relevant 
literature. 


A. A brief review 

For one-dimensional lattice models with local random 
potentials, analytical results for the conductance and 
density of states have been obtained using the Fokker- 
Planck (FP) equation obeyed by the Lyapunov exponents 
of the ID transfer matrix EI) t3 9 l 48 t {H] 

These results establish that the metallic phase [in 


the so-called quasi-one-dimensional (ID) limit defined by 
taking the number N of one-dimensional channels to in¬ 
finity in a suitable wajP^-j is always unstable to disorder 
away from the band center. Upon increasing the length 
of the ID geometry, disorder drives a crossover to an (An¬ 
derson) insulating phase.^ESl However, the band center 

e = 0 is anomalous for seven of the ten AZ symmetry 
c l a g CCG j36 H39l48l5m 55l 

On the one hand, the chiral symmetry classes AIII, 
BDI, and CII display an even-odd effect in the parity of 
the number N of ID transverse channels by which the 
disorder-averaged conductance at the band center alter¬ 
nates between exponential decay for even N and alge¬ 
braic decay for odd iVP^On the other hand, the disorder- 
averaged thermal conductance at the band center for any 
one of the superconducting symmetry classes D and Dill 
shows a mean conductance with algebraic decay for any 
^yj5il These results hold for Gaussian disorder with van¬ 
ishing mean. 

The density of states in the neighborhood of the band 
center also signals that the band center is a critical en¬ 
ergy for the symmetry classes AIII, BDI, CII, D and Dill 
separating two ID Anderson insulating phases, as it fails 
to follow the dependence predicted by random matrix 
theory (i.e., a suppression of the density of states re¬ 
sulting from the enhanced level repulsion caused by the 
spectral symmetry about the band center) j :j 8 l 52 l 54 l I n con¬ 
trast, the disorder-averaged conductance is exponentially 
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(a) AIM, BDI, CM 

N odd 


(b) D, Dill 

N even 



(c) A, Al, All, C, Cl 

g A 


Insulator 


> 

m 


FIG. 3. Qualitative quantum phase diagrams for ID disordered wires with the quasiparticle energy fixed at e = 0. The 
horizontal axis is parametrized by the characteristic value m G R of the disorder-averaged Dirac masses allowed by the symmetry 
class. The vertical axis is parametrized by the characteristic value g > 0 taken by the strength of the generic static and local 
random mass disorder allowed by the symmetry class. Arrows on the phase boundaries and on the horizontal axis indicate 
flows under renormalization group transformations. 


suppressed with the length of the wire in the symmetry 
classes C, Cl, A, AI, and All for all energies while the 
density of states does not deviate from its expected be¬ 
havior at the diffusive (unstable) fixed point. 

The Fokker-Planck equation (the DMPK equation) 
obeyed by the Lyapunov exponents describing the trans¬ 
fer matri x in ID contains universal data of geometric ori- 
gin.EHIS 7 in each of the ten AZ symmetry classes, the ID 
transfer matrix defines a noncompact symmetric space, 
of whi ch th e Lyapunov exponents are the radial coordi¬ 
nates.^® 9 An infinitesimal increase in the length of the 
disordered region for one of the ten symmetry classes in¬ 
duces an infinitesimal Brownian motion of the Lyapunov 
exponents that is solely controlled by the multiplicities of 
the ordinary, long, and short roots of the corresponding 
classical semi-simple Lie algebra under suitable assump¬ 
tions on the disorder (locality, weakness, and isotropy be¬ 
tween all channels). The long and short roots have a very 
special meaning with regard to Anderson localization. 
When the ID transfer matrix describes the stability of 
the metallic phase of noninteracting fermions perturbed 
by one-body local random potentials in the bulk of a ID 
lattice model, the multiplicity of the short root entering 
the Brownian motion of the Lyapunov exponents always 
vanish EMMHsHIS Moreover, the multiplicities of the long 
roots also vanish for t he five ID symmetry classes AIII, 
BDI, D, Dill, and CIl! 36 l 51 l This vanishing is the signa¬ 
ture of the existence of critical points separating topo¬ 
logical insulating phases in ID wires. However, when the 
ID transfer matrix describes the ID boundary of a 2D 
topological band insulator moderately perturbed by local 
random potentials, the multiplicities of the short roots is 
non-vanishing in the Brownian motions of the Lyapunov 
exponents in the five AZ symmetry classes A, AIIpB D, 
Dill, and C. Correspondingly, the conductance is of order 
one along the infinitely long boundary, i.e., the insulating 
bulk supports extended edge states. These extended edge 
states can be thought of as realizing a ID ballistic phase 
robust to disorder. The interpretation of a DMPK equa¬ 
tion with non-vanishing multiplicity of the short root is 


that it describes transport along the boundary of a 2D 
topological insulator. 

The same conclusions for the stability of a bulk ID 
metallic phase}® 1 '®! and for the stability of a ID bal¬ 
listic phase along the boundary of a 2D band insulator 
perturbed by moderate disordei^®l also follow from de¬ 
scribing these disordered noninteracting fermionic mod¬ 
els with NLSMs. The role played by the multiplicities of 
the short roots in order to evade Anderson localization 
within the ID transfer matrix approach is played by the 
presence of a topological term in the NLSM with the tar¬ 
get spaces corresponding to the symmetry classes A, All, 
D, Dill, and 


B. 7T 0 (P d= i iT .) =: 


According to the third column of Table m the topo¬ 
logical space V d=1 r=r N that parametrizes the space 
of Dirac masses is the union of the Grassmannian man¬ 
ifolds indexed by the number of ID channels N for the 
symmetry classes AIII, BDI, and CII. 

Figure [3](a) is the schematic phase diagram for the chi¬ 
ral symmetry classes AIII, BDI, and CII at the band 
center e = 0. A ID disordered wire with N channels is 
characterized by A^ + l distinct insulating phases that are 
indexed by the set of half-integers and integers with N 
odd and even, respectively. The flow of g along the phase 
boundary separating two neighboring insulating phases 
in Fig. p i a) is to the infinitely strong disorder fix points 
g —> The topological argument for this prop¬ 

erty is the following. The topological space V d=1 r=r N 
parametrized by the Dirac masses allowed in any one of 
the three chiral symmetry classes is disconnected. One 
Dirac mass is singled out, the Dirac mass whose Dirac 
matrix commutes with all other Dirac masses present 
when N > 1: /3 0 defined in Eq. (3.9a). The value m taken 
by averaging this mass over its probability distribution is 
the control parameter m that triggers a quantum phase 
transition between two distinct neighboring insulating 
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phases in Fig. [3^a) . At m = 0, a disordered wire is in the 
v = 0 localized phase when TV is even and on the phase 
boundary between the v = ± | insulating phases when TV 
is odd. This qualitative even-odd effect in TV that origi¬ 
nates with the parity of the number of disconnected com¬ 
ponents TV +1 of V d=1 r=r N is captured quantitatively 
by the Fokker-Planck approach IH H 39 I 48 IIIH111 the NLSM 
approach ®^®® 7 or the application of nonconnnutative 
geometrjfsM 71 ! and quantum-entanglement techniques.^ 
One-parameter scaling is lost at the phase boundaries 
in the vicinity of the band centerPSHSS] Recently, ana¬ 
lytic results encoding a two-parameter scaling similar to 
that characterizing the integer quantum Hall effect were 
obtained for the symmetry classes AIII and BDI from 
the supersymmetric NLSM approach.^ Moreover, the 
fan-like phase diagram depicted in Fig. |4|a), in which 
TV +1 localized phases emerge from m = g = 0, has been 
confirmed for the symmetry class BDI from a numerical 
study of the Lyapunov exponents of transfer matrices.^ 
Note that a nonlinearity in the band dispersion makes 
the fan-like diagram asymmetric about m = 0 in tight- 
binding models. 

The anomalous behavior at the phase boundaries can 
be deduced from the results obtained previously by 
the Fokker-Planck approach.^SHSUHUHtMl ^t any pi lase 
boundary the dimensionless Landauer conductance g 
shows anomalous dependence on the wire length L ( 3 > 
Nt), 


In <7 oc — 



In q oc-In | - 

y 2 {N£ 


(4.1a) 


and the density of states near e = 0 has the strongest 
divergence, the Dyson singularity, 


p(e) oc 


1 


£r[ln(|e|r)] c 


(4.1b) 


where r = TV 2 £, £ is the mean free path, and the Fermi 
velocity has been set to unity. The universal critical be¬ 
havior (4.1) is a manifestation of an (unstable) infinite- 


disorder quantum critical point controlled by Griffiths 
effects E2H23 

In an insulating phase the dimensionless conductance 
has exponential dependence on L , 


In oc — —, (4.2) 

where the localization length £ is of order TV £ and di¬ 
verges at phase boundaries. The exponent in the power- 
law e dependence of p(e) increases from —1 as m is 
changed from its value at a phase boundary. Half way be¬ 
tween two consecutive transitions, the density of states is 
the closest to the prediction from random matrix theory, 

p(s) oc (£r ) m ° _1 | ln(£r)|, (4.3) 


where m 0 = 1, 2,4 are the multiplicities of the ordinary 
roots for the chiral symmetry classes BDI, AIII, and CII, 


respectively.® However, the density of states (4.3) is not 
quite the one expected from random matrix theory, for it 
acquires a multiplicative logarithmic correction, as shown 
in Ref. [351 The density of states interpolates between 
these two limiting functions as a function of the control 
parameter m .' 8 


C. n 0 (V d=hr ) = : 


According to the third column of Table m the topo¬ 
logical space V d=1 r=r N that parametrizes the normal¬ 
ized random Dirac mass in ID has two disconnected com¬ 
ponents for the BdG symmetry classes D and DHL Ac¬ 
cordingly, there should be two topologically distinct in¬ 
sulating phases separated by a phase-boundary line, at 
which the localization length and the density of states 
dive rge. Tuning the control parameter m def ined in 
Sec. |III b| to zero selects the phase-boundary line , 51 52 Ml 
around which one-parameter scaling is broken 152155 77 
Figure |3](b) is the phase diagram for the BdG symme¬ 
try classes D and Dill at the Fermi level £ = 0. A ID 
disordered wire is characterized by two insulating phases 
that are separated by a phase boundary along which the 
disorder strength g flows to the infinitely strong disorder 
fix point, the Dyson fix point. This phase boundary is lo¬ 
cated at m = 0. The topological argument for this prop¬ 
erty is the following. The topological space V d _ 1 r=r N 
parametrized by the Dirac masses allowed in any one of 
these two BdG symmetry classes is made of two discon¬ 
nected components. For any realization of the random 
potential these two components are indexed by the Z 2 
index defined in Eq. (3.8). Hence the control parame¬ 


ter m that drives the quantum phase transition between 
the two distinct insulating phases can be chosen to be 
the disorder-averaged value over this PfafRan. On the 
other hand, had we started from the symmetry class BDI 
and weakly broken the time-reversal and chiral symme¬ 
tries by the addition of weak perturbations that bring the 
Dirac Hamiltonian to one in the symmetry class D, we 
may then keep the control parameter m of the symmetry 
class BDI along the horizontal axis in parameter space. 
There follows the same phase diagram as in Fig. [4](a) with 
the caveat that the TV + 1 unperturbed localized phases 
with the Z topological numbers v become TV +1 localized 
phases with alternating Z 2 topological numbers v after 
perturbation.® 

The anomalous behavior at the phase boundary can 
be deduced from the results obtained previously by 
the Fokker-Planck approach.®® The anomalous depen¬ 
dence on the wire length L of the dimensionless Landauer 
conductance g is the same as Eq. (4.la I, and the density 


of states p(e) exhibits the Dyson singularity (4.1b). In 
the insulating phases the conductance has the exponen¬ 
tial dependence on L , Eq. (4.2). The density of states 


p{e) has the power-law dependence on the excitation en¬ 
ergy £ with the exponent continuously varying from —1 
of the Dyson singularity to the exponent of the random 
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matrix theorj® (0 and 1 for class D and Dill, respec¬ 
tively), as |m | is increased from the critical point m = 0. 
The fact that the symmetry classes D and Dill share with 
the chiral symmetry classes the same Dyson singularity 
is known as “superuniversality” )5£L 7S 


This result is also understood nonperturbatively as a con¬ 
sequence of the fact that the 2D NLSM that encodes the 
physics of Anderson localization in the symmetry class 
AI has a coupling, the inverse of the value of the longitu¬ 
dinal conductivity in the diffusive regime, that flows to 
strong coupling ! 82 * 83 * 


D. 7T 0 (V d=ltr ) = {0} 

According to the third column of Table m the topo¬ 
logical space V d=1 r _ r N that parametrizes the normal¬ 
ized random Dirac masses is path-connected for the sym¬ 
metry classes A, AI, All, C, and Cl. Figure |3jc) is the 
schematic phase diagram for the symmetry classes A, AI, 
All, C, and Cl at e = 0 and to all symmetry classes 
away from e = 0. A ID disordered wire is always lo¬ 
calized in these symmetry classes. The topological argu¬ 
ment for this property is the following. The topological 
space V d=1 r parametrized by the Dirac masses allowed in 
any one of these five symmetry classes is path connected. 
Consequently, the sign of m has no topological mean¬ 
ing. The allowed Dirac masses select a unique insulating 
phase. 

The density of states in the vicinity of the quasiparticle 
energy e = 0 is affected by the enhanced level repulsion 
in the BdG symmetry classes C and Cl, but this effect 
is captured by random matrix theon® s i nC e the local¬ 
ization length is never divergent) 51 * 52 * p c (e) oc e 2 and 

Pci(e) oc |e|. 


V. APPLICATION TO 2D SPACE 

This section is dedicated to Anderson localization in 
the ten 2D AZ symmetry classes when a local disorder 
is present. First, we review known facts about these 
2D disordered systems. Unlike the Fokker-Planck ap¬ 
proach for which we only know how to extract controlled 
analytical results in ID* 79 ! the NLSM approach delivers 
perturbative renormalization-group (RG) results in any 
dimensions close to 2DP§1 We then show how the topol¬ 
ogy of classifying spaces combined with perturbative RG 
flows allow to deduce the phase diagram that delimits 
the topologically distinct localized and metallic phases 
for each 2D AZ symmetry class. This analysis delivers 
a unifying explanation for many past analytical and nu¬ 
merical results regarding Anderson localization in the 2D 
symmetry classes when the disorder is local. 


A. A brief review 

In the symmetry classes AI, all the single-particle 
states are exponentially localized upon perturbation by 
a random potential. This result is understood perturba- 
tively as the phenomenon of weak localization by which 
the first non-vanishing quantum correction to the longi¬ 
tudinal conductivity caused by a disorder is negative! 81 ! 


In the symmetry class A, the phase diagram of single¬ 
particle states consists of multiple of localized phases 
that are distinguished by their quantized Hall conduc¬ 
tivity and are separated by ph ase bo und aries at which 
the localization length diverges.® ItHGSHI The existence 
of critical points separating neighboring localized phases 
can be understood as a result of the topological term 
which can be added to the 2D NLSM for class AP^l These 
critical states levitate with increasin g diso rder strength 
until they “annihilate” upon merging! 84 * 85 * 

In the 2D symmetry class All, only the single-particle 
states close to the edges of the Bloch bands are exponen¬ 
tially localized in the presence of a weak disorder. 82 83 
As the disorder strength is increased, the lower and up¬ 
per mobility edges separating in energy the Lifshitz tails 
from the extended single-particle states approach each 
other until they merge and all single-particle states are 
exponentially localized for sufficiently strong disorder. 

Consequently, the metallic phase for noninteracting 
electrons propagating in 2D space is unstable to the pres¬ 
ence of any local disorder in the symmetry classes A and 
AI, while it is stable for a sufficiently weak disorder in 
the symmetry class All, as long as the chemical potential 
is sufficiently far from the unperturbed band edges. 

These three outcomes for the competition between the 
kinetic energy and the local random potential in 2D space 
are not exhaustive in the presence of a particle-hole or a 
chiral symmetry, as it does in the symmetry classes AIII, 
BDI, D, Dill, CII, C, and CI. Among the BdG symme¬ 
try classes, symmetry class CI is insulating very much 
in the same way as symmetry class AI isF® Symmetry 
class C can display a spin IQHE very much in the same 
way as symmetry class A can display the charge IQHEP^ 
The metallic phases in the symmetry classes D and Dill 
are stable to weak disorder EZHSU Symmetry class D can 
display a thermal IQHE very much in the same way as 
symmetry class A can display the charge IQHE.* 88 9 ™ 92 J 
Symmetry class Dill can realize a thermal Z 2 topological 
insulato r ver y much in the same way as symmetry class 
All can.* 28 * 93 Finally, the chiral classes AIII, BDI, and CII 
are anomalous at the band center with diverging local¬ 
ization length and density of states, while they crossover 
to the symmetry classes A, AI, and All, 
any non-vanishing value of the chemical 


respectively, for 
potential!® 41 


All of these features that are of topological origin can 
be understood qualitatively from the zeroth and first ho- 
motopy groups from Table IV, as we now explain below. 



















17 


TABLE IV. Interplay between the topology of the normalized Dirac masses in a given AZ symmetry class and Anderson 
localization when space is two-dimensional. For each AZ symmetry class from the first column, (i) the second column gives the 
minimal rank r min of the Dirac matrices for which a Dirac mass is allowed by the symmetries, (ii) the third column gives the 
corresponding topological space for the normalized Dirac masses V d - 2 r=r jv, (iii) the fourth column gives the corresponding 
zeroth homotopy group 7 t 0 (Vj_o r =r /v)> (iv) the fifth column gives the corresponding first homotopy group ‘K 1 (V d=2 r=r N ), 
and (v) the sixth column gives the corresponding panel from Fig. HI Each panel from Fig. Elis a phase diagram for quasiparticles 
at the Fermi energy, which is fixed to e = 0. Parameter space is two-dimensional with the mean value m £ R for the characteristic 
Dirac masses as horizontal axis and the characteristic disorder strength g > 0 as vertical axis. The last column characterizes 
transport through the dependence of the mean conductivity as a function of the linear length L of the disordered region when 
m = 0. Insulating phases are characterized by a mean conductivity that decays exponentially with L. A metallic phase has a 
mean conductivity that grows with L. A critical phase has a mean conductivity of the order of e 2 /ft. The entries “even-odd” 
indicate a critical or metallic phase for odd N and an insulating phase for even N. The entries “Gade singularity” indicate a 
diverging density of state associated with n 1 (V d=2 r ) = Z. In each homotopy column, the three entries Z hold for N larger than 
an integer (infinity included) that depends on the order of the homotopy group and the classifying space, the two entries Z 2 
hold for N larger than an integer that also depends on the order of the homotopy group and the classifying space. The entry 
0 is a short hand for the group { 0 } made of the single element 0 . 
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B. Implications of the topology of the classifying 
spaces 


Equation (3.1l with d = 2 and with the rank r = r min N 
is our starting point. It is assumed that the local ran¬ 
dom disorder entering the Dirac Hamiltonian is the most 
general r x r Hermitian matrix with random identically 
independently distributed (iid) matrix elements up to the 
constraints imposed by the selected AZ symmetry classes. 
We fix the chemical potential /i = 0 and discuss trans¬ 
port properties of quasiparticles of energy e = 0 for all 
ten symmetry classes. 

For each AZ symmetry class in the first column of Ta¬ 
ble |IV[ we give the minimum rank r min for which the 
Dirac Hamiltonian admits a Dirac mass, the classifying 
space V d=2 r , its zeroth homotopy group ^o(Vd=2 r jv )> 
its first homotopy group ^i{V d=2r N ), and its phase 
diagram for TV = 1 and N = 2 . All other entries hold 
for any N. Each panel from Fig. [3] is a phase diagram 
in a two-dimensional parameter space. The horizontal 
and vertical axes are the characteristic mean value m of 
the Dirac masses and the characteristic disorder strength 
g, respectively, introduced in Sec. HI The last column 
characterizes transport through the dependence of the 
mean conductivity as a function of the linear length L of 
the disordered region when m = 0 and £ = 0. Insulating 


phases are characterized by a mean conductivity that de¬ 
cays exponentially with L. A metallic phase has a mean 
conductivity that grows with L. Critical phases have a 
mean conductivity that is non-vanishing but finite for 
L —► oo. The entry “even-odd” indicates that the parity 
of N selects either an insulating or a critical phase, de¬ 
pending on whether N is even or odd, respectively. The 
new entry “Gade singularity” compared to the ID Table 
IH] signals a conductivity that is independent of the dis¬ 
order strength and a divergence of the density of states 
upon approaching the band center. This new entry is a 
signature of 7 r 1 (V ’ d=2 r ) = Z, as we are going to explain. 

The main result of this section is that, according to 
Sec. |III[ there are deviations in Table |IV| away from the 
insulating behavior found in the 2D symmetry classes AI 
and Cl. These deviations are attributed to the zeroth and 
first homotopy groups of the normalized Dirac masses 
being either Z 2 or Z. 


C. TV 0 {V d= 2 t r) = Z 

According to Table |IV[ the 2D AZ symmetry classes 
A, D, and C are associated with Grassmannian mani¬ 
folds. Hence the zeroth homotopy group of the topolog¬ 
ical spaces associated with the normalized Dirac masses 
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FIG. 4. Qualitative quantum phase diagrams for 2D random Dirac Hamiltonians with the quasiparticle energy fixed at e = 0. 
The horizontal axis is parametrized by the characteristic value m G R of the disorder-averaged Dirac masses allowed by the 
symmetry class. We choose the probability distribution (3.101 for the case of TV = 2 in the symmetry class D. The vertical 
axis is parametrized by the characteristic value g > 0 taken by the strength of the generic local disorder allowed by the 2D 
AZ symmetry class. The zeroth homotopy group associated with the space of normalized Dirac mass fixes the indexing of the 
topologically distinct insulating phases. The first homotopy group of the normalized Dirac masses determines if the density of 
states is regular or singular at the band center. The existence of a metallic phase, the flows along phase boundaries, and the 
dimensionality of the phase boundaries all follow from perturbative renormalization-group calculations. 


in these AZ symmetry classes is Z. According to Sec. m 
their phase diagrams must host N + 1 topologically dis¬ 
tinct insulating phases and N transition lines separating 
them. 


The arguments of Sec. m follow from the zeroth ho¬ 
motopy group of the Grassnrannian manifolds. The 2D 
symmetry classes A and C share the same trivial first ho¬ 
motopy group. In view of the fact that the 2D metallic 
phase is perturbatively unstable (due to weak localiza¬ 
tion) in both symmetry classes (a fact not encoded in 
the homotopy groups of Table IV), we expect 2D phase 
diagrams in the 2D symmetry classes A and C that are 
similar to the one in the ID chiral symmetry classes, the 


one shown in Figs. §a) and[4])b) for TV = 1 and N = 2, 
respectively. 

Contrary to the 2D symmetry classes A and C, the 
metallic phase of the 2D symmetry class D is stable (due 
to weak antilocalization) to weak disorder. Hence the 
2D symmetry class D is expected to host a (thermal) 
metallic phase when the characteristic disorder strength 
g is sufficiently large, in addition to N + 1 topologically 
distinct insulating phases. When the metallic phase is 
robust to weak disorder, topologically distinct localized 
phases can be separated by a metallic phase rather than 
by a mere phase boundary line. It is the RG flow of 
the random Dirac mass in the Dirac Hamiltonian that 
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determines which of the metallic phase or the critical 
boundary separates two topologically distinct insulating 
phases. Now, the characteristic coupling g for the dis¬ 
order is marginally irrelevant in the vicinity of the crit¬ 
ical point m = 0, g = 0 for the symmetry class D with 
N = 1, as is reviewed in Appendix[B] Consequently, there 
is a RG flow along the phase boundary between topo¬ 
logically distinct insulating phases into the critical point 
m = 0, g = 0 when N = 1. More generally, phase bound¬ 
aries separate topologically distinct insulating phases at 
small g. Thus, starting from any of the N + 1 insulating 
phase separated by phase boundaries, we have a transi¬ 
tion driven by increasing g into the metallic phase of the 
2D symmetry class D, as shown in Figs. gc) (N = 1) 
an dgd) (N = 2). With regard to the case of N > 1, say 
N = 2 as depicted i n Fig. gd), we are choosing the prob¬ 
ability distribution ( 3.10[ ) to parametrize the horizontal 
axis of the phase diagram instead of the probability dis¬ 
tribution ( |3.9[ ) used to parametrize the horizontal axis of 
the phase diagramgc) when N = 1. (See Appendix|B 3 c 


for a more detailed discussion of the relationship between 
the irrelevant RG flow of g and the expectation values 
of mass terms.) Moreover, the Grassmannian manifold 
for the 2D symmetry class D also differs from that for 
the symmetry classes A and C in that its first homotopy 
group is nontrivial and given by Z 2 . In the N = 2 case, it 
is the path-connected component 0(2)/0(l) x 0(1) from 
the space of Dirac masses (the v = 0 insulating phase) 
that has a nontrivial first homotopy group and can thus 
host Z 2 vortices, each supporting a Majorana zero mode. 
Such zero modes would lead to a Griffiths-like singularity 
at e = 0 in the density of states. However, Z 2 vortices are 
not expected to exist in any insulating phase whose space 
of Dirac masses has a trivial first homotopy group. For 
example, in the N = 2 case, the point-like components 
0(2)/0(2) x 0(0) or O(2)/O(0) x 0(2) from the space 
of Dirac masses (the v = ±1 insulating phases) have a 
trivial first homotopy group. Hence, if it is the v = ±1 
insulating phases of Dirac fermions in the symmetry class 
D with N = 2 that are realized in the vicinity of e = 0, 
then no Griffiths-like singularity is expected at e = 0 in 
the density of states. In other words, we expect a density 
of states that deviates from the one predicted by random 
matrix theory in the symmetry class D for any insulat¬ 
ing phase supporting Majorana zero modes bound to Z 2 
Dirac mass vorticesPS 


These predictions are consistent with the following nu¬ 
merical studies of 2D quantum network models. The 
cases N = 1 and N = 2 in the symmetry class A are 
the best known applications of quantum network models. 
Chalker and Coddington studied numerically in Ref. ITU 
the 2D quantum network model co rresp onding to the case 
N = 1 in the symmetry class A. 1595 ! Lee and Chalker 
studied numerically in Ref. (5B] the 2D quantum network 
model corresponding to the case N = 2 in the symme¬ 
try class A with their investigation of spin-degenerate 
Landau levels and found two quantum phase transitions 
separating three insulating phases of 2D quantum mat¬ 


ter. Two-dimensional quantum network models for the 
symmetry class C have also been studied analytically and 
numerically. 97 ®®! They deliver the same phase diagrams 
as for the 2D symmetry class A. The case of N = 1 in 
the symmetry class D can also be regularized by a 2D 
quantum network model with nodal scattering matrices 
of appropriate symmetry and rank, 97 ! or by two-band lat¬ 
tice models for Majorana fermions!®!! The former model 
shows a critical boundary with the attractive fixed point 
corresponding to the Dirac Hamiltonian at m = g = 0 
that separates two insulating phases of 2D quantum mat¬ 
ter for not too strong g, as shown in Fig. |4][c) when 
N = 1. The two insulating phases are unstable to the 
metallic phase upon increasing g, as shown in Fig. |4|c). 
There have been numerical studies on the critical prop¬ 
erties around the multicritical point where the tw o in- 
sulating and one metallic phases meet at g > 0. 92 93 
For N = 2, there are three insulating phases on top of 
which sits a metallic phase, as shown in Fig. |4](d). More¬ 
over, quasi-zero modes bound to Z 2 vortices in the Dirac 
masses are present in any insulating phase for which the 
first homotopy group is Z 2 , the v = 0 phase in Fig. |l|d) 
when N = 2. 


D. 7To(Fd=2,r) — Z 2 


According to Table |IV[ the 2D AZ symmetry classes 
Dill and All are associated with disconnected topologi¬ 
cal spaces with two path-connected components, i.e., the 
zeroth homotopy group of the topological spaces associ¬ 
ated with the normalized Dirac masses in these AZ sym¬ 


metry classes is Z 2 . According to Sec. Ill their phase 
diagrams must host two topologically distinct insulating 
phases and one critical boundary or one metallic phase 
separating them at m = 0. 

From the point of view of topology, the normalized 
Dirac masses in the 2D symmetry class Dill differ from 
those in the 2D symmetry class All in that the former has 
the first homotopy group Z 2 when N > 2, while the latter 
has a trivial first homotopy group for all N = 1,2,.... 
This difference manifests itself through a density of states 
that deviates from the one predicted from random matrix 
theory in the 2D insulating phases of the symmetry class 
Dill when N = 2,3,... [Fig. |4jf)]. The metallic phase 
is stable (due to weak antilocalization) to the presence 
of weak disorder in both 2D symmetry classes. Now, the 
characteristic coupling g for the disorder is marginally 
relevant in the vicinity of the critical point m = g = 
0 for symmetry classes Dill and All, as is reviewed in 
Appendix [B] Consequently, the metallic phase separates 
two topologically distinct insulating phases for any non¬ 
vanishing g > 0 all the way down to the critical point 
m = g = 0 and sits on top of the insulating phases as 
depicted in Figs. 0e)’B f )’ and 0g)- 

These predictions are consistent with the following 
pu blished works. Studies of random Dirac Hamiltoni- 
a np 2 HTol anc j a ne t wor k model 3 ®! in the 2D symmetry 
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class All have confirmed the robustness of the metallic 
phase all the way down to the critical point m = g = 0 in 
Fig. |4^g). The numerical study in Ref. [53] of a 2D quan¬ 
tum network model for the symmetry class Dill is also 
consistent with Fig. [4je). 


E- MV d=2ir ) = {0} 


According to Table |IV[ the 2D AZ symmetry classes 
AIII, AI, BDI, CII, and Cl are associated with path- 
connected topological spaces, i.e., the zeroth homotopy 
group of the topological spaces associated with the nor¬ 
malized Dirac masses in these AZ symmetry classes is 
the trivial group {0}. According to Sec. Ill their phase 
diagrams must host no more than one insulating phases. 

The symmetry classes AI and Cl do not support a 2D 
metallic phase, as it is unstable due to weak localization. 
For this reason, their phase diagrams in Fig. |4|h) consist 
of a single insulating phase. An example of a localized 
phase in class Cl is a dirty d 2 _„2 superconductor with 
TRS and spin SU( 2) rotation symmetry.^ 8 d 

The chiral symmetry classes AIII, BDI, and CII stand 
out in 2D space because (i) the quantum corrections to 
their longitudinal Drude conductivities vanish to any or¬ 
der in perturbation theory precisely at the band center, 
and (ii) their density of states diverges at the band center 
e = 0, as was shown by Gade and Wegner for a sublat¬ 
tice model perturbed by disorder with chiral symmetry 
in Refs. EH and ED Items (i) and (ii) are also true for 
random Dirac Hamiltonians in the 2D sym met ry classes 
AIII, BDI, and CII at the band center B 2 | 43 | i06 | 

On the one hand, the topology of the classifying spaces 
of the 2D chiral symmetry classes AIII, BDI, and CII is 
not predictive regarding item (i). On the other hand, 
the fact that the first homotopy group of the normal¬ 
ized Dirac masses in the 2D symmetry classes AIII, BDI, 
and CII is Z provides an explanation for the diverg¬ 
ing density of states upon approaching the band cen¬ 
ter alternative to the one based on the RG calculations 
from Refs. HHM31 and 1IUB1 or to the one based on Grif¬ 
fith effects from Ref. EH The physical interpretation of 
7r l( V dt™r=4) = ■Xl(yJSl,r=4) = = Z is that 

one finds the (noncontractible) unit circle S 1 in the topo¬ 
logical space of normalized Dirac masses from the 2D 
symmetry classes AIII, BDI, and CII. This subspace is 
generated by two anticommuting mass matrices as was 
pointed out in the studies of charge fractionalization in 
graphene.^ 07113 Such a unit circle, corresponding as it 
is to a pair of normalized masses, supports point-like de¬ 
fects, 17(1) vortices, which can bind an integer number 
of zero modes if isolated! 114 A finite but dilute density 
of vortices binds midgap states resulting in the singular 
density of states at the band center. We show the phase 
diagrams of the 2D chiral symmetry classes at e = 0 in 
Fig. j4j(i). In the insulating phase the density of states 
exhibits a Griffiths singularity p(e) oc \e\ a at the band 
center e = 0 with the exponent a varying continuously 


with |m | from a = — latm = 0 (Gade singularity) to the 
value expected from random matrix theory a = 1,0,3 
for class AIII, BDI, and CII, respectively. We note that 
the 2D symmetry class CII differs from the 2D symme¬ 
try classes AIII and BDI when £ / 0 and the disorder 
is weak. Indeed, any deviation of the chemical poten¬ 
tial away from the band center is a relevant symmetry¬ 
breaking perturbation that drives a crossover between the 
symmetry class CII to the symmetry class AIL As spin- 
orbit coupling favors antiweak localization, the metallic 
phase is stable to weak disorder in the symmetry class 
AIL 


VI. APPLICATION TO 3D SPACE 

In this section, we are going to employ the machinery 
of classifying spaces in order to deduce the qualitative 
phase diagrams capturing the physics of Anderson local¬ 
ization when the effects of local disorder are not always 
strong, or marginally strong as they were in ID and 2D, 
respectively. The lowest dimension in which such disor¬ 
der can be treated perturbatively in all ten AZ symme¬ 
try classes is d = 3. We shall show that the machinery 
of the classifying spaces applied to the normalized Dirac 
masses is predictive regarding the phase diagrams when 
combined with the fact that the characteristic coupling 
g for the disorder is irrelevant at the critical Dirac point 
for any d > 3. Before doing so however, we briefly review 
salient past results on Anderson localization in 3D. 

The metallic phase is stable in 3D space for the stan¬ 
dard symmetry classes A, AI, and AIL A sufficiently 
strong disorder is required to stabilize the insulating 
phase jnaHni The same is true for all chiral and sym¬ 
metry classes at the band center. 

The insulating phases in 3D space for the standard 
symmetry classes A and AI are topologically trivial. 
There are two topologically distinct insulating phases in 
the 3D symmetry class AIlP"^HSD At the band center, 
four more nonstandard symmetry classes support topo¬ 
logically distinct insulating phasesP^o] 


A. Higher homotopy groups and topological defects 


More interestingly from the point of view of this pa¬ 
per, we predict that there are 3D symmetry phases with 
topologically trivial insulating phases that support a sin¬ 
gular density of states at the band center due to the fact 
that their normalized Dirac masses support topological 
defects, as we are going to explain now. 

To investig ate t his possibility, we start from the Dirac 
Hamiltonian (3.11 with d = 3 and the rank r = r ■ N. 


It is assumed that the local disorder entering the Dirac 
Hamiltonian is the most general r x r Hermitian matrix 
with random iid matrix elements up to the constraints 
imposed by the selected AZ symmetry classes. We fix the 
chemical potential p = 0 and discuss transport properties 
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TABLE V. Interplay between the topology of the normalized Dirac masses in a given AZ symmetry class and Anderson 
localization when space is three-dimensional. For each AZ symmetry class from the first column, (i) the second column gives 
the minimal rank r min of the Dirac matrices for which a Dirac mass is allowed by the symmetries, (ii) the third column gives the 
corresponding topological space for the normalized Dirac masses V d - 3 r=r N , (iii) the fourth column gives the corresponding 
zeroth homotopy group n 0 (V d=3 r=r N ), (iv) the fifth column gives the corresponding first homotopy group ■n 3 (V d=3 r=r N ), 
(v) the sixth column gives the corresponding second homotopy group n 2 (V d=3 r=r N ), (vi) the seventh column gives the 
corresponding panel from Fig. [ 5 ] and (vii) the last column signals if the density of states per unit energy and per unit volume 
is singular at the Fermi energy, which is fixed to e = 0. A singular density of states is associated with nontrivial point defects 
characterized by i^ 2 (V d - 3 r=r N ). In each homotopy column, the three entries Z hold for N larger than an integer (infinity 
included) that depends on the order of the homotopy group and the classifying space, the two entries Z 2 hold for N larger than 
an integer that also depends on the order of the homotopy group and the classifying space. The entry 0 is a short hand for the 
group {0} made of the single element 0. 
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FIG. 5. Qualitative quantum phase diagrams for 3D random Dirac Hamiltonians with the quasiparticle energy fixed at e = 0. 
The horizontal axis is parametrized by the characteristic value m G R of the disorder-averaged Dirac masses allowed by the 
symmetry class. We choose the probability distribution (3.101 for the case of N = 2 in the symmetry classes AIII, Dill, and Cl. 
The vertical axis is parametrized by the characteristic value g > 0 taken by the strength of the generic local disorder allowed 
by the 3D AZ symmetry class. The zeroth homotopy group of the topological space associated with the normalized Dirac 
mass fixes the indexing of the topologically distinct insulating phases. The existence of a metallic phase, the flows along phase 
boundaries, and the dimensionality of the phase boundaries all follow from perturbative renormalization-group calculations. 


of quasiparticles of energy e = 0 for all ten symmetry 
classes. 


According to Sec. IIF when 7T d _ 1 (V (l ) has a nontriv¬ 


ial entry, point defects in the normalized Dirac masses 
may bound zero modes. If so, they contribute to a 
singular density of states. From the transition rule in 
Eq. (2.23), in any dimension d of space, ^ d _i{V d ) = Z 
for the chiral symmetry classes AIII, BDI, and CII, while 
7 T d _ 1 (V d ) = Z 2 in the BdG symmetry classes D and DHL 
Hence, in any dimension d of space, a random Dirac 
Hamiltonian of the form (3.1) displays a singular density 


of states at the band center in the chiral symmetry classes 
AIII, BDI, and CII and in the BdG symmetry classes D 
and Dill due to the proliferation of midgap states bound 
to Z and Z 2 point defects, respectively. Topological de¬ 
fects of higher dimensions d def = 1 ,..., d — 1 also bind 
midgap states that contribute to the density of states 
around e = 0. However, they are not expected to cause 
a singularity of the density of states at e = 0. Indeed, ei¬ 
ther the density of extended defects scales with the linear 
size L of space as L ~ dd and is thus subextensive, or the 
midgap states bound to compact extended defects with 
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the characteristic linear length L de( are a level spacing of 
order 1/L def away from e = 0. Above the dimensional¬ 
ity of space d = 3, all AZ symmetry classes support at 
least one nontrivial homotopy group. The form of the 
singularity of the density of states at the band center 
depends parametrically on the statistical distribution of 
the Dirac masses. It is generically nonuniversal, except 
at the critical points that govern the transition between 
distinct phases. 


B. ir 0 (V d=3r ) — Z 


The normalized Dirac masses realize the Grassman- 
nian manifolds in the 3D symmetry classes AIII, Dill, 
and Cl according to Table |V} We deduce that these 3D 
symmetry classes support N + 1 topologically distinct 
insulating phases separated by N boundaries by com¬ 
bining the arguments of Sec. m with the fact that the 
characteristic coupling g for the disorder is irrelevant at 
the critical Dirac point for any d > 3, see Fig. §a). 
Since the characteristic coupling g is irrelevant, we must 
distinguish the case when N = 1 from the case when 
N > 1. When N = 1, we choose the probability distri¬ 
bution (3.9). Wh en N > 1, we choose the probability 
distribution (3.10). Indeed, we recall that if we label 
the horizontal axis (the clean limit) of the phase dia¬ 
gram with m defined by Eq. ( |3.9b ), then the horizontal 
axis supports no more than two insulating phases with 
v = ±y, respectively. Given that g is irrelevant in the 
vicinity of m = g = 0, the choice (3.9b) delivers two lo¬ 
calized phases below the metallic dome. In other words, 
the phase diagram when N > 1 with the choice (3.9) is 
identical to that of Fig. [5ja) with the caveat that the 
index v = ±4 is to be replaced with the index v = 
for the two localized phases. To display N + 1 > 2 lo¬ 
calized phases in the phase diagram, it is necessary to 
have N + 1 > 2 insulating phases along the horizon¬ 
tal axis (the clean limit) of the phase diagram. This 
is achieved with the choice (3.10) that insures that the 
eigenvalues of the Dirac mass matrix j3 0 in Eq. (3.10a) 
competes with N nondegenerate eigenvalues of the ma¬ 
trix M 0 in Eq. (3.10c) The second homotopy group of the 
normalized Dirac masses is Z in the 3D symmetry class 
AIII according to Table [V] The first and second homo¬ 
topy groups of the normalized Dirac masses are Z 2 in the 
3D symmetry class Dill according to Table [VJ Hence the 
density of states (the number of states per unit energy 
and per unit volume) is singular at the band center due 
to the proliferations of midgap states that are bound to 
point-defects of the normalized Dirac masses for the 3D 
symmetry classes AIII and DHL The 3D symmetry class 
Cl has trivial first and second homotopy groups so that 
the density of states is nonsingular at the band center. 


C. n 0 (V d=3r ) — Z 2 

The normalized Dirac masses are the union of two 
path-connected compact topological spaces in the 3D 
symmetry classes All and CII according to Table [V] 
We deduce that these 3D symmetry classes support two 
topologically distinct insulating phases separated by one 
boundary by combining the arguments of Sec. |III| with 
the fact that the characteristic coupling g for the disor¬ 
der is irrelevant at the critical Dirac point for any d > 3. 
A metallic phase sits on top of these two topologically 
distinct insulating phases. Thus, we deduce the phase 
diagram for classes All and CII shown in Fig. [5](b). This 
phase diagra m is c onsistent with the self-consistent Born 
approximation 121 and non pert urbative numerical studies 
in the symmetry class AIL 11 ' On the other hand, had we 
started from the symmetry class Dill and weakly broken 
the particle-hole and chiral symmetries by the addition 
of weak perturbations that bring the Dirac Hamiltonian 
to one in the symmetry class All, we may then keep the 
control parameter m of the symmetry class Dill along the 
horizontal axis in parameter space. This would give the 
same phase diagram as in Fig. [5](a) for N > 1 with the 
caveat that the N + 1 unperturbed localized phases with 
the Z topological numbers become N +1 localized phases 
with alternating Z 2 topological numbers after perturba¬ 
tion. 

The first and second homotopy groups of the normal¬ 
ized Dirac masses from the 3D symmetry classes All and 
CII differ. The second homotopy group is trivial in the 
3D symmetry classes All according to Table [V] Hence 
the density of states is not singular in the 3D symme¬ 
try classes AIL The second homotopy group is Z in the 
3D symmetry classes CII according to Table [V] Hence 
the density of states is singular in the 3D symmetry 
classes CII. The first homotopy groups of the 3D symme¬ 
try classes All and CII are Z 2 and trivial, respectively. 
In the former case, the density of states receives a con¬ 
tribution from midgap states bound to the line defects. 
This contribution need not be singular at e = 0. 


D. Tr 0 (V d=3r ) — {0} 

The normalized Dirac masses are path-connected com¬ 
pact topological spaces in the 3D symmetry classes A, AI, 
BDI, D, and C according to Table [V] These symmetry 
classes support a single insulating phase, see Fig. [5])c) . 
Of these, only the 3D symmetry class AI has the triv¬ 
ial group for its three homotopy groups. The normalized 
Dirac masses in the 3D symmetry classes BDI and D have 
the second homotopy groups Z and Z 2 for N sufficiently 
large, respectively. Thus, the 3D symmetry classes BDI 
and D have densities of states that are singular at the 
band center because of the proliferation of midgap states 
bound to Z and Z 2 point-like defects in the normalized 
Dirac masses, respectively. Line defects of the normal¬ 
ized Dirac masses are also allowed in the 3D symmetry 
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classes A, D, and C. They also contribute to the density 
of states, but they need not give rise to a singularity at 
e = 0. We note that the density of states computed for 
a network model from the symmetry class C was found 
to be nonsingular in Ref. [1001 


VII. BOUNDARIES OF TOPOLOGICAL 
INSULATORS 

So far, we have been concerned with the interplay be¬ 
tween topology, local symmetries, and Anderson localiza¬ 
tion for d-dimensional random massive Dirac Hamiltoni¬ 
ans that capture the effects of local, and smooth disorder 
of d-dimensional lattice models such as quantum network 
models at low energies and long wavelengths. We have 
explained in Sec. |IID| how Table[l]can be used to identify 
in each dimension d of space the five AZ symmetry classes 
that support topologically distinct insulating phases of 
d-dimensional quantum matter. There is an alternative 
road to the classification of TIs and TSCs that was ex¬ 
plained in Sec. |II D[ which we now review for the sake of 
completeness. 

Any (d — l)-dimensional boundary of a d-dimensional 
topological insulators or superconductors is immune to 
Anderson localization. 28 This immunity is captured by 
the presence of a topological term in the NLSM that 
captures the physics of Anderson localization on the 
(d — l)-dimensional boundary subject to the symme¬ 
try constraint imposed by one of the ten AZ symmetry 
classes.^ Alternatively, this property holds for those AZ 
symmetry constraints that prohibit the existence of any 
Dirac mass matrix of rank r = r min < r min entering the 
random Dirac Hamiltonian that encodes single-particle 
transport at low energies and long wavelength on any 
(d — l)-dimensional boundary. In Appendix |A 4| we for¬ 
mulate the conditions for the presence (or absence) of 
Dirac mass matrices. We show that the conditions for the 
absence of normalized Dirac masses in (d— 1)-dimensions 
are equivalent to the conditions for a nontrivial zeroth ho- 
motopy group of the topological spaces associated with 
the normalized Dirac masses in d-dimensions. This al¬ 
ternative derivation of the tenfold way for TI/TSCs is 
captured by Table |Vl| for 2D and 3D spaces. 


VIII. SURFACE STABILITY OF TOPOLOGICAL 
INSULATORS WITH REFLECTION SYMMETRY 

A. General discussions 

This section is devoted to the stability of topologi¬ 
cal crystalline insulators (TCIs) with a reflection sym¬ 
metry when perturbed by a local random potential that 
retains the reflection symmetry on average. The disor¬ 
der is assumed weak in the bulk in that its characteristic 
strength is much smaller than the band gap from the d- 
dimensional bulk. In this limit, the (d — 1)-dimensional 


TABLE VI. Immunity to Anderson localization in five out 
of the ten AZ symmetry classes along any boundary of a 
topological insulator (TI) or a topological superconductors 
(TSC). This immunity is a consequence of the absence of any 
Dirac mass matrix with the rank r = r min < r min (see Ap¬ 
pendix |A 4[ ) for a random Dirac Hamiltonian capturing the 
low-energy and long-wavelength effects of local disorder on 
the boundary in the corresponding AZ symmetry class. 


class 

edge of 2D TI/TSC 

surface of 3D TI/TSC 

A 

ballistic 

- 

AIII 

- 

metallic 

AI 

- 

- 

BDI 

- 

- 

D 

ballistic 

- 

Dill 

ballistic 

metallic 

All 

ballistic 

metallic 

CII 

- 

metallic 

C 

ballistic 

- 

Cl 

- 

metallic 


boundary may be treated in isolation and represented 
in terms of massless Dirac fermions perturbed by local 
random potentials obeying appropriate symmetries. The 
material SnTe is an example of a TCI and the alloying 
with Pb in Sn 1 _ a ,Pb a ,Te can be modeled by a weak local 
random potential. 

By definition, a d-dimensional TCI is a band insula¬ 
tor that supports on any one of its (d — 1 (-dimensional 
boundaries that is invariant under the crystalline sym¬ 
metry, here a reflection symmetry, boundary states that 
disperse across the bulk gap with a linear (Dirac) dis¬ 
persion relation in the close vicinity to their band cross- 
ingdU j n other words, the reflection symmetry protects 
(d — 1)-dimensional boundary massless Dirac cones. 

A classification of d-dimensional TCI with reflec¬ 
tion symmetry has been obtained in Ref. [35] (see also 
Ref. Il23|) from studying the zeroth homotopy groups 
of the topological spaces associated with the normal¬ 
ized Dirac masses entering d-dimensional massive Dirac 
Hamiltonians obeying nonlocal reflection symmetries in 
addition to the AZ local symmetries. For each of the 
AZ symmetry classes ordered as in the first columns of 
Tables VII and VIII[ the results from Ref. S3 are sum¬ 
marized in columns three to six when space is two- and 
three-dimensional, respectively. For completeness, the 
second columns of Tables IVIII and IVIIII are the same as 
in Table [I] when d = 2 and d = 3, respectively. Hence, 
for the AZ symmetry classes BDI, Dill, CII, and Cl, 
the zeroth homotopy group of the topological space as¬ 
sociated with the normalized Dirac masses when reflec¬ 
tion symmetry is not imposed (column two) is supple¬ 
mented by the zeroth homotopy group of the topological 
space associated with the normalized Dirac masses when 
the reflection symmetry BPtMc holds (columns three to 
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TABLE VII. Topological classification of crystalline insulators (TCIs) in 2D with reflection symmetry and their edge stability 
when the reflection symmetry is retained on average, a situation for which we use the acronym STCI (statistical topological 
crystalline insulator). For each AZ symmetry class, we give in the second column the topological classification of insulators 
(TIs). Columns three to six are the TCIs for the symmetry classes BDI, Dill, CII, and Cl with the reflection symmetry R r,T ’ ,lc , 
where the superscripts ?? T = ± and r] c = ± specify if the operation for reflection commutes (+ sign) or anticommutes (— sign) 
with reversal of time T and exchange of particles and holes C, respectively. Here, we choose the convention (fJ^T’^c ) 2 — -pi. 
The third and fourth columns for the symmetry class AIII/(AI, AII)/(D, C) are the topological classification of the TCIs with 
the reflection symmetry I?’ 1 ’’ 7 , where the superscript 77 specifies if the operation for reflection commutes (+) or anticommutes 
(—) with chiral/time-reversal/particle-hole transformation, respectively. The quotation marks for the entry “Z 2 ” corresponding 
to a TCI with the R~'~ symmetry in class CII indicates instability to interavalley scattering . 35 125 The last column gives the 
STCIs. The superscript f over the trivial zeroth homotopy groups in the STCI column signals when the edge states of a TCI 
are unstable to the presence of disorder that respects the reflection symmetry on average. 


AZ symmetry class 

2D TI 

R+’ + 

2D TCI 

R > R + -~ 

R-'+ 

2D STCI 

A 

Z 

0 

- 

- 

- 

0 

AIII 

0 

z 

0 

- 

- 

z 

AI 

0 

0 

0 

- 

- 

0 

BDI 

0 

z 

0 

0 

0 

z 

D 

z 

z 2 

0 

- 

- 

Z 2 

Dill 

Z 2 

Z 2 

z 

0 

Z 2 

Z 2 

All 

Z 2 

0 

z 2 

- 

- 

o f 

CII 

0 

z 

“Z 2 ” 

0 

0 

z 

C 

z 

0 

0 

- 

- 

0 

Cl 

0 

0 

z 

0 

0 

o f 


six). We choose the convention (/pr^c ) 2 = -pi. Here, 
any one of the superscripts rj T = zt and r) c = ± takes 
the value + (—) if the operation of reflection commutes 
(anticommutes) with the operation T for the reversal of 
time and the operation C for the exchange of particles 
and holes, respectively. For the AZ symmetry classes 
AIII/(AI, AII)/(D,C), the third and fourth columns indi¬ 
cate the topological classification with the reflection sym¬ 
metry R v,ri , where the superscript 77 takes the value + 
(—) if the reflection operator commutes (anticommutes) 
with the operation for chiral/time-reversal/particle-hole 
transformation, respectively. For the symmetry class A, 
the third column shows the classification in the presence 
of a reflection symmetry. Any entry indicates the 
absence of the symmetries under reversal of time T or 
under exchange of particles and holes C. 

The question we want to address is what is the fate of 
those boundary massless Dirac cones protected by a com¬ 
bination of AZ and reflection symmetries if a boundary 
is perturbed by a local random potential belonging to a 
statistical ensemble such that (i) there exist realizations 
of the random potential that break the reflection sym¬ 
metry, (ii) even though the reflection symmetry holds on 
average. The acronym STCI for statistical topological 
crystalline insulator is used when items (i) and (ii) are 
met. 


The answer is found in the last columns of Tables IVIII 
and |VIII| in 2D and 3D space, respectively. Each en¬ 
try of the column STCI from the Tables |VII| and VIII 


gives the zeroth homotopy group of the topological space 
associated with the normalized Dirac masses in (d — 1)- 
dimensional space that are allowed by imposing no other 
symmetries than the ones from the AZ symmetry class to 
which the d-dimensional STCI belongs. The intuition for 
this result is that a d-dimensional STCI supports either a 
metallic or critical phase on any of its (d— l)-dimensional 
boundary that is invariant under the crystalline opera¬ 
tion if Anderson localization in the corresponding (d— 1)- 
dimensional AZ symmetry class is preempted for topolog¬ 
ical reasons. 

The interpretation of the non-vanishing zeroth homo¬ 
topy groups in the columns three to six from Tables m 
and |VIII| is the following. Whenever the representation 
of the d-dimensional TCI in terms of a Dirac Hamilto¬ 
nian with Dirac matrices of rank r is the minimal one, 


r = r min , ( 8 . 1 a) 

i.e., no mass matrix is compatible with the symmetries 
for any smaller rank of the Dirac matrices (the actual 
value of r min depends on the AZ class and the type of 
the R symmetry), then there exists a normalized mass 
matrix /3 with the following two properties. First, /3 is 
unique up to a sign. Second, (3 is invariant under the 
reflection represented by R , 

Rf3R~ 1 =p. (8.1b) 

Two Dirac Hamiltonians that differ through the sign ±(3 
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TABLE VIII. Topological classification of crystalline insulators (TCIs) in 3D with reflection symmetry and their surface 
stability when the reflection symmetry is retained on the average, a situation for which we use the acronym STCI (statistical 
topological crystalline insulator). For each AZ symmetry class, we give in the second column the topological classification 
of insulators (TIs). Columns three to six are the TCIs for the symmetry classes BDI, Dill, CII, and Cl with the reflection 
symmetry R' 1t,v c j where the superscripts rj T = ± and r] c = ± specify if the operation for reflection commutes (+ sign) or 
anticommutes (— sign) with reversal of time T and exchange of particles and holes C , respectively. Here, we choose the 
convention (J^'t^c) 2 = -pi. The third column for symmetry class A indicates the Z classification of TCIs with reflection 
symmetry. The third and fourth columns for the symmetry class AIII/(AI, AII)/(D, C) are the topological classification of 
the TCIs with the reflection symmetry R n,v , where the superscript r\ specifies if the operation for reflection commutes (+) or 
anticommutes (—) with the chiral/time-reversal/particle-hole transformation, respectively. The quotation mar ks for t he entry 
“Z 2 ” corresponding to a TCI with the RT'~ symmetry in class C indicates instability to intervalley scattering? 35 123 The last 
column gives the STCIs. The superscript j over the trivial zeroth homotopy groups in the STCI column signals when the 
surface states of a TCI are unstable to the presence of disorder that respects the reflection symmetry on the average. 


AZ symmetry class 

3D TI 

R + ’+ 

3D TCI 

R R - 

R~’+ 

3D STCI 

A 

0 

z 

- 

- 

- 

z 

Am 

Z 

0 

Z 

- 

- 

o f 

AI 

0 

0 

z 

- 

- 

o f 

BDI 

0 

0 

0 

0 

z 

o f 

D 

0 

z 

0 

- 

- 

z 

Dill 

z 

Z 2 

0 

z 

z 

Z 2 

All 

Z 2 

Z 2 

z 

- 

- 

Z 2 

CII 

Z 2 

0 

z 2 

Z 2 

z 

o f 

C 

0 

z 

“Z 2 ” 

- 

- 

z 

Cl 

z 

0 

0 

z 

z 

o f 


with which (3 enters select two topologically distinct crys¬ 
talline insulating phases in d-dimensional space. The 
topologically nontrivial TCI phase is the one that sup¬ 
ports a massless Dirac cone on any one of the (d — 1)- 
dimensional boundaries that are invariant under the re¬ 
flection symmetry. Such a massless Dirac cone is derived 
from the Dirac Hamiltonian a ■ p with the Dirac matri¬ 
ces a = (a l ,..., of rank r min /2 and the boundary 

momentum p £ 1 . 

Assume that a topologically nontrivial phase of a d- 
dimensional TCI is selected by the sign with which the 
mass matrix /3 enters a Dirac Hamiltonian of minimal 
rank. Furthermore, assume that, without the protec¬ 
tion arising from the reflection symmetry, this TCI is 
a topologically trivial band insulator. This implies that 
a normalized mass matrix (3' must exist such that f3' an- 
ticommutes with (3, 

(£,/3'} = 0, (8.2a) 


symmetry by m' /?' is small, 

|m'| \m\. (8.3) 

If the sign of m selects a topologically nontrivial crys¬ 
talline insulating phase in the limit m' = 0, then any 
(d — l)-dimensional boundary left invariant under the re¬ 
flection symmetry must support a massive Dirac cone 
with the gap proportional to \m'\. Such a massive Dirac 
cone is derived from the Dirac Hamiltonian a- p + fh' /3' 
with the Dirac matrices a = (a lt ..., a d _ 1 ) and (3' of 
rank r min /2 and the boundary momentum p £ K d_1 . 
The Dirac mass matrix /3' at the boundary originates 
from the Dirac mass matrix f3' in the bulk. 

Assume that (3' is multiplied by a random function 

m! : R. d —> R, x h>■ m!(x ) (8.4a) 

that satisfies the bound 


and j3' is odd under the transformation 

RI3'R- 1 = -I3' (8.2b) 


by the operation R. 

Consider a d-dimensional Dirac Hamiltonian with the 
pair of mass terms m /3 and nn! f3' that satisfies Eqs. (8.1) 
and (8.21. Assume that the breaking of the reflection 


sup |to , (:e)| < \m\, (8.4b) 

cceR d 

where m selects a topologically nontrivial crystalline in¬ 
sulating phase in the limit for which the reflection sym¬ 
metry is recovered. The fate of the boundary states that 
are protected by the reflection symmetry in the presence 
of the weak reflection-symmetry-breaking perturbation 
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(8.4) is the same as the fate of massless Dirac fermions 
in (d — 1) dimensions perturbed by a generic local ran¬ 
dom potential with vanishing mean [by item (ii)] that 
obeys no other symmetries than the ones from the AZ 
symmetry class to which the d-dimensional TCI belongs. 
The latter stability analysis can be done with the meth¬ 
ods from Secs. mm When the classifying space of the 
normalized Dirac masses in (d — 1 [(-dimensions has non¬ 
trivial zeroth homotopy group (listed in the last columns 
of Tables VII and VIII| ), it has two path-connected com¬ 
pact subspaces indexed by the ambiguity to assign a sign 
to the unique normalized Dirac mass matrix ft of the 
minimal rank (half the minimal rank of the bulk TCI, 
i.e., r min /2). The anticommutation relation (8.2b) im¬ 
plies that this pair of path-connected compact subspaces 
are homeomorphic to each other by the reflection trans¬ 
formation R. 

Assume now that we are given a statistical ensemble 
of random Dirac Hamiltonians of rank r = r min , each 
realization of which would be the same d-dimensional 
TCI, perturbed by reflection-symmetry-breaking ran¬ 
dom potentials, with the reflection symmetry R obey¬ 
ing Eq. (8.1b). We may then decompose any (d — 1)- 
dimensional boundary that is invariant under the re¬ 
flection symmetry into patches labeled by the sign of 
the Dirac mass m!(x) multiplying the /3' matrix that 
smoothly varies along this boundary, say A if the sign 
is + and B if the sign is —, as was done in Fig. [2] The 
relative size of all the regions labeled A compared to all 
the regions labeled B in Fig. [2] is fixed by the mean value 
of the mass term with the mass matrix /3'. This differ¬ 
ence in sizes vanishes if the reflection symmetry holds on 
average. The existence of a (d— 2)-dimensional subman¬ 
ifold between regions A and regions B that (i) percolates 
across the (d— 1 (-dimensional boundary and (ii) binds a 
midgap state signals a metallic or critical phase of (d— 1)- 
dimensional boundary. 

The argument of the previous paragraph breaks down 
if we relax the assumption r = r min to r = r min N with 
N = 2, 3,.... Indeed, when the rank of the Dirac ma¬ 
trices is the integer multiple N = 2,3,... of the min¬ 
imum rank r min , the operation for reflection does not 
need to interchange regions A and B from Fig. 0 any¬ 
more. Imposing the reflection symmetry on average is 
not sufficient anymore to ensure that the difference in 
the sizes of all regions labeled A and all regions labeled 
B in Fig. [2] vanishes. Stronger assumptions obeyed by 
the statistical ensemble of random masses must be ful¬ 
filled to guarantee the existence of a (d — 2(-dimensional 
path between all regions A and B that percolates across 
the (d — 1)-dimensional boundary depicted in Fig. [2] and 
binds a midgap state. To appreciate this point, we con¬ 
sider the following two cases. 

Consider first the case of the zeroth homotopy group 
Z 2 and N = 2. We can then make the choice (3.15) 
with M a 4 x 4 matrix. The reflection tran sforma tion 
R still anticommutes with the mass matrix (3.15b), as 
assumed in Eq. (8.2b). However, it leaves the Pfaffian 


(3.15d) used to index the topological sector of the mass 
term unchanged. Hence R cannot be used to construct a 
homeomorphism between the subspaces belonging to dis¬ 
tinct topological sectors when N = 2. Imposing that the 
randomness preserves the reflection symmetry on aver¬ 
age is not sufficient to enforce a metallic or critical phase 
of (d — l)-dimensional boundary. 

We can generalize this argument to arbitrary N. For 
odd V, preserving the reflection symmetry on average 
leads to a vanishing val ue for the disorder average of the 
Z 2 topological number ( |3.15d ). Either a metallic or criti¬ 
cal phase follows at the boundary of a STCI. For even N , 
preserving the reflection symmetry on average does not 
constrain the disorder average of the topological number 
(3.15d). We thus expect an even-odd effect (a sensitivity 
to the parity of N ) with regard to the stability of surface 
states of a STCI when the zeroth homotopy group is Z 2 . 

In the case of Z for the zeroth homotopy group of the 
topological space associated with the normalized Dirac 
masses in (d — 1) dimensions, there exists a unique nor¬ 
malized Dirac mass matrix, up to a sign, that commutes 
with all other permissible mass matrices. This Dirac ma¬ 
trix corresponds to /3 0 in Eq. (3.9a). The disorder aver¬ 
age value of this mass term vanishes (i.e., m = 0) due 
to statistical (average) reflection symmetry of the STCI. 
According to Tables m and m we have the following 
even-odd effect. The (d — 1 (-dimensional boundary of a 
STCI is metallic or critical for odd N, while it is localized 
for even N, as depicted in Figs. 3(a) and 4(a)-4(d). 

Next, we consider two examples of a three-dimensional 
TCI. 


B. Example in the 3D symmetry class All 

First, we consider a 3D TI in the symmetry class AIL 
We then impose a crystalline symmetry by demanding 
invariance under the reflection 


x = (x, y, z) {-x, y, z) (8.5a) 

s = (s x , S y , s z ) (s x , -Sy, -s z ), (8.5b) 

for the space and spin degrees of freedom, respectively. 

To this end, we choose the 4x4 massive Dirac Hamil¬ 
tonian 

'H(k) := ( k x s y — k y s x ) g> a x + k z s 0 ® a y +ms 0 ®d r 

( 8 - 6 ) 

The 2 x 2 unit matrix s 0 and the Pauli matrices s act on 
the spin-1/2 degrees of freedom. They enter the kinetic 
energy through a spin-orbit coupling. The 2x2 unit 
matrix a 0 and the Pauli matrices er act on the orbital 
degrees of freedom. Reversal of time is represented by 
conjugation of R{k) with 


T := is,, 


>cr 0 K, 


(8.7) 


where complex conjugation is denoted by K. The massive 
Dirac Hamiltonian (8.6) realizes a 3D insulator in the 
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symmetry class All with the Dirac matrices 
a — ( s y ® ~ s x 0 a xi s o ® a y) > /3 = s 0 0 cr z , (8.8) 

as 

TT-L(-k)T~ l =H{k). (8.9) 

The rank r = 4 of the Dirac matrices is the minimal one 
for 3D topological insulators from the symmetry class All 
represented by single-particle massive Dirac Hamiltoni¬ 
ans, 


r • =4. 

' min 


( 8 . 10 ) 


where r 0 is a 2 x 2 unit matrix and r are Pauli matrices, 
admits the normalized mass term 

j3' = s z ®a x ®T y (8.13b) 

that anticommutes with 


/3 = s 0 


(8.13c) 


and is invariant under the transformation T = \s y 0 cr 0 0 
r 0 K that implements reversal of time. 

We now impose the crystalline symmetry originating 
from the operation (8.5b) on the 3D massive Dirac Hamil¬ 
tonian (8.13aI. To this end, we define 


According to Ref. OS} the 3D Dirac Hamiltonian (8.6) 
of rank r = r min = 4 captures the low-energy and long- 
wavelength limit in the vicinity of one out of four inequiv¬ 
alent L points from the face-centered-cubic (fee) Brillouin 
zone of the band insulator SnTe for one sign of m or the 
band insulator PbTe for the opposite sign of m. The 
rank r = r min = 4 arises because (1) both SnTe and 
PbTe have a simple rocksalt structure (two atoms in the 
unit cell), (2) both display a spin-orbit coupling, and (3) 
because the splitting in energy between the conduction 
and valence bands is minimum at the center of the eight 
hexagonal faces of the fee Brillouin zone, the so-called L 
points. SnTe and PbTe belong to the same trivial insu¬ 
lating phase in the entry All from the second column of 
Table |VHTJ because this phase supports an even number 
of gapless surface states. However, the plane 


R x~ s x® a o> { R xY = s 0 ®cr 0 , { R x,T} = 0, 

(8.14a) 

and demand that 


R x R(~K, k y > k z) R x = R{ k x, k y, K)~ 


(8.14b) 


One observes that the mass matrix (3 is compatible with 
the crystalline symmetry originating from the operation 
(8.5b). However, the crystalline symmetry originating 
from the operation (8.5b) is not compatible with the mass 
matrix /?' defined by Eq. (8.13b) if 


R x '■= s x 0 (7 0 0 r 0 , T := is y 0 cr 0 0 r 0 K, (8.15a) 
is combined with the symmetry condition 
Rf H{- k x, k y, k z)®T 0 R- = U(k x ,k y ,k z )®T 0 . (8.15b) 


(k x , k y , k z ) = ( 0 , k y , k z ) ( 8 . 11 ) 


-X^^K) = (° , k y, k z) 

from the face centered cubic Brillouin zone is invariant 


under the reflection (8.5b). On this so-called mirror 


plane, the massive 3D Dirac Hamiltonian (8.6) simplifies 


to the direct sum of two massive 2D Dirac Hamiltoni¬ 
ans of rank r = 2 from the 2D symmetry class D whose 
relative Chern numbers differ by one. 

The topological space of the normalized Dirac masses 
in the 3D symmetry class All follows here from the ex¬ 
tension problem 


(8.12a) 


^-^ 2,3 — { l ' RJ, OC x , CXy , Oi z j" 

~^ ^ 3,3 = {JP,T,TJ~ a x , cx, yjQx} 

with the classifying space 

R ± = O(N) (8.12b) 

[see Eq. j [2.5[ ) and Table [TI], whose zeroth homotopy 
group i: ° 


This conclusion can be generalized in the following way. 
If we impose the symmetry under the reflection (8.5b) 


that anticommutes with reversal of time to any single¬ 
particle 3D Dirac Hamiltonian belonging to the symme¬ 
try class All, we mu st sup plement the Clifford algebra for 
class All [recall Eq. (2.17)] with a new generator J a x R~ 
that squares to +1 times the identity operator. The topo¬ 
logical space of the normalized Dirac masses in the 3D 
symmetry class All obeying this reflection symmetry fol¬ 
lows here from the extension problem 

= R J ) ^ y ’ ^ z > R^xRx } 

^ ^-^ 3,4 = { J ft i R> RJ^ ex x , Oiy , oc z , J ex x R x } 

(8.16a) 

with the classifying space 

R 0 = jv lim o U() r =0 {0(Af)/[0(n) x O(N-n)]} (8.16b) 


[see Eq. (2.5) and Table [TT] , whose zeroth homotopy 
group is 


7r 0 (i?i)=Z 2 , (8.12c) 


7T 0 (i? 0 ) = Z, (8.16c) 


as the solution. We recall that the topological number is 
Z 2 rather than Z from the fact that the tensoring 

R(k) —>• H(k) 0 r 0 , (8.13a) 


as the solution. 

The zeroth homotopy group Z for the normalized Dirac 
masses in the 3D symmetry class All obeying the crys¬ 
talline symmetry under the operation (8.5b) is the reason 
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why an even number of surface massless Dirac cones rep¬ 
resented by rank two Dirac matrices is protected, when 
no such protection is operative if we relax this crystalline 
symmetry and deal with the 3D symmetry class AIL In 
other words, an insulator from the 3D symmetry class 
All that would support an even number of surface Dirac 
cones in the clean limit is topologically trivial because a 
generic local random potential from the symmetry class 
All gaps out any even number of surface Dirac cones. 
Since the clean limit of SnTe can be approximated at 
low energies by the massive Dirac Hamiltonian (8.6) ten- 
sored with the N x N unit matrix where TV = 4 is the 
even number of Dirac cones centered around the N = 4 
inequivalent L points in the fee Brillouin zonewe con¬ 
clude that SnTe is a topologically trivial insulator under a 
generic single-particle local perturbation belonging to the 
symmetry class All, while SnTe is a topologically non¬ 
trivial insulator whose topologically distinct phases are 
indexed in Z, if the crystalline symmetry under the oper¬ 
ation (8.5b) is imposed. Accordingly, the corresponding 
entries in Table IVIIII are Z 9 and Z in the second and 


It remains to explain the 


fourth columns, respectively, 
entry Z 2 in the last column. 

We now assume that the 8x8 Dirac Hamiltonian 
(8.13a) is perturbed by a local random potential that 
(l) belongs to the symmetry class All, (2) may break 
the crystalline symmetry (8.14) for any given realization, 
(3) but preserves the crystalline symmetry (8.14) on the 
average. Hence the mass matrix /?' defined in Eq. (8.13b) 
enters the random potential multiplied by a random func¬ 
tion m! : R 3 —> R, x H>• m'(x) that varies smoothly in 
space but averages to zero everywhere in 3D space. This 
mass matrix anticommutes with the mass matrix (8.13c). 
Any surface of the TI that is invariant under the opera¬ 
tion (8.5b) can be divided, for any realization of the ran¬ 
dom function m!, into the regions A in which sgnm'(x) 
is + and the regions B in which sgnm'(x) is —, as is 
done in Fig. [2j Along the one-dimensional boundary be¬ 
tween two path-connected regions A and B, m'(x) = 0. 
Such an edge binds a pair of midgap helical edge states. 
The typical difference between the area of all the regions 
A and the area of the regions B is fixed by the mean 
value of the random function m! . This typical difference 
is zero if the mean value of the random function m! van¬ 
ishes everywhere on the surface. If so, a pair of helical 
midgap states are bound to a one-dinrensional path that 
percolates across the surface. Hence they are extended. 

The parameter x in the alloy Sn 1 _ a ,Pb a .Te is associated 
with a local random potential as the selection of the sites 
on which Sn is to be substituted by Pb is essentially a 
random Poisson process to a first approximation. In the 
clean l imit x —» 0, the projection of the massive 3D Dirac 
cones (8.6) at the four inequivalent L points onto one of 
the mirror symmetric surface Brillouin zone delivers an 
even number of 2D massless Dirac cones. Specifically, 
four Dirac cones appear on the [001] surface (N = 2) LL^AI 
while six Dirac cones appear on the [111] surface (N = 3). 
The disorder in the alloy Sn 1 _ x Pb x Te belongs to the sym¬ 


metry class AIL It is also believed to preserve, on aver¬ 
age, the reflection symmetry present in the clean limit 
(8.6). Thus, the alloy Sn 1 _ a ,Pb x Te is an example of the 
STCIsfMlMI According to the even-odd effect in the case 
of the zeroth homotopy group Z 2 discussed in Sec. VIII A[ 


the gapless surface states of Sn 1 _ a .Pb z ,Te are stable (ex¬ 
tended) on the [111] surface, but they are unstable (lo¬ 
calized) on the [001] surface for generic disorder that pre¬ 
serves the reflection symmetry on average. Nevertheless, 
the four Dirac cones at the [001] surface can also be sta¬ 
ble if the C 4 rotation symmetry interchanging four Dirac 
cones is preserved on average. Alternatively, if one of 
the two reflection symmetries, which are present at the 
[001] surface, is broken by a lattice distortion, only two 
out of the four Dirac cones remain massless. This pair 
of surviving massless Dirac cones is immune to Anderson 
localization in the presence of weak disorder that pre¬ 
serves one of the two reflection symmetries at the [001] 
surface on average. 


C. Example in the 3D symmetry class AI 

Second, we consider a 3D TI in the symmetry class AI 
upon which we impose a crystalline symmetry under 

x = (x,y,z) (~x,y,z) (8.17) 


that anticommutes with reversal of time. 

To this end, we choose the 8x8 massive Dirac Hamil¬ 
tonian (r min = 8 in the 3D symmetry class AI) 

'H(k) := k x A 233 + k y X 023 + k z X 002 + E m j Xj . 

i 

(8.18a) 

Here, 

x ^up ’■= s p ® o v <g> T P , n,v,p = 0,1, 2, 3, (8.18b) 

where each quartet (s 0 , s 1: s 2 , s 3 ), (u 0 , oq, <r 2 , cr 3 ), and 
(r 0 , t 1; r 2 , t 3 ) of 2 x 2 Hermitian matrices is built out 
of the 2x2 unit matrix together with the three Pauli 
matrices, and I denotes the collective index taking the 
values 


1 = 001,013,133,333. (8.18c) 

Represent reversal of time T by complex conjugation 

K, 


T := X 000 K. 


(8.19) 


Represent the reflection (8.17) that anticommutes with 
T by 


R~ := X 020 . (8.20) 

One verifies that 

TTL(-k) T = H(k). (8.21) 
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Furthermore, one verifies that imposing 


R x'U-{-k x ,k y ,k z ) R x =H(k x ,k yl k z ) 


x ’ y z) 


( 8 . 22 ) 


demands of the fou r mass es entering the massive 8x8 
Dirac Hamiltonian (8.18a) that three of them vanish, 


m 0 i 3 = m! 3-, = m 


133 


'333 


= 0 . 


(8.23) 


Since only one normalized mass matrix is compatible 
with the AZ symmetry class AI and the reflection sym¬ 
metry (8.22) in 3D space, it follows that the topological 


space of normalized Dirac masses in the minimal repre¬ 
sentation has two compact and path-connected compo¬ 
nents {±A 001 }. Consequently, a bou ndary that is left 
invariant by the reflection symmetry (8.17) must host 


boundary states that are described by a massless 4x4 
Dirac equation in the low-energy and long-wavelength 
limit. Accordingly, we must find a non-vanishing entry 
in the corresponding entry of the R column of Ta¬ 
ble |VIII| In order to conclude that this entry is Z, we 
must double the rank of the Dirac Hamiltonian (8.18a) 


and verify that one mass commutes with all other allowed 
masses. Alternatively, we may reach the same conclu¬ 
sion from the fact that, without the reflection symmetry, 
the extension problem Cl 0 5 —► Cl 15 has the classify¬ 
ing s pace i? 5 with n 0 (R 5 ) = 0 for solution according to 
Eq. ( |2.5| ), while, with the reflection symmetry, the ex¬ 
tension problem Cl 0 6 —> Cl 16 has the classifying spac e 
R 4 with 7r 0 (i? 4 ) = Z for solution accord ing to Eq. ( |2.5| ). 
The fact that the reflection symmetry (8.22) in the 3D 


symmetry class AI preempts more than one mass matrix 
is the reason for which a 3D STCI in the symmetry class 
AI has no protected surface states when the reflection 
symmetry holds only on average. 


D. Relation to surface states of weak topological 
insulators 

Surface states of weak topological insulators (WTIs) 
are robust to local disord er that pre serves the transla¬ 
tion symmetry on average ESHtlSi 127 This is understood 
in the same manner as with STCIs. In a STCI, the aver¬ 
age reflection symmetry tunes the surface states to be on 
the critical boundary or in the metallic phase at m = 0 
for odd N. In a WTI, the average translation symmetry 
tunes the surface states to be on the critical boundary or 
in the metallic phase at m = 0 for odd N, as was shown 
for Z 2 WTI of IV = 1 in Ref. tl041 When N is an even 
integer, the average translation symmetry is not suffi¬ 
cient for the surface states of a WTI to evade Anderson 
localization, as in STCIs. 


IX. SUMMARY 

The problem of Anderson localization for noninteract¬ 
ing fermions realizing random Dirac Hamiltonians was 


revisited from the perspective of the topologies of the 
spaces associated with Dirac mass matrices allowed by 
symmetries. 

It was first shown that the topological space V of 
normalized Dirac masses can be characterized in a sys¬ 
tematic fashion with the help of the homotopy groups 
7r 0 (V), 7r 1 (V'), ..., 7r d _ 1 (U) by imposing on V the Altland- 
Zirnbauer (AZ) symmetries obeyed by all realizations of 
the local disorder. (The AZ symmetries are local.) For 
any spatial dimension d and for any n = 0,... , d — 1, 
one finds three d-dimensional AZ symmetry classes for 
which 7 T n (V) = Z and two d-dimensional AZ symmetry 
classes for which Tr n (V) = Z 2 . In particular, the zeroth 
homotopy group is nontrivial whenever the topological 
space V is disconnected. The zeroth homotopy group 
7r 0 (U) indexes the path-connected subspaces of V and 
enumerates the topologically distinct insulating phases 
in d dimensions. By considering a mapping of conduct¬ 
ing channels, formed along boundaries between topo¬ 
logically distinct insulating regions, to a d-dimensional 
random quantum network model, we showed (i) which 
pairs of topologically distinct insulating phases are sep¬ 
arated by either a metallic or a critical phase and (ii) 
that the control parameter across this intervening metal¬ 
lic phase or critical boundary involves either the disorder- 
averaged value taken by the Dirac mass that commutes 
with all other Dirac masses when n 0 (V) = Z, or the 
disorder-averaged value over the Pfaffian of the Dirac 
mass when 7r 0 (V) = Z 2 . Qualitative phase diagrams in a 
two-dimensional parameter space spanned by this control 
parameter and the characteristic disorder strength were 
deduced for each d-dimensional AZ symmetry when the 
chemical potential is at the band center by combining 
7r 0 (V) with a one-loop RG calculation. 

We have also explained under what conditions the self¬ 
averaging density of states per unit energy and per unit 
volume is singular in the vicinity of the band center and 
when this singularity becomes universal. A sufficient con¬ 
dition for a singular density of states is that the homotopy 
group 7r d _ 1 (V r ) is nontrivial and point defects bounding 
zero-energy states may appear. However, topology alone 
is not sufficient to predict the nature of the singularity. 

Finally, we have extended the study of the interplay be¬ 
tween the AZ symmetry classes, topology, and Anderson 
localization by allowing crystalline (reflection) symme¬ 
tries to hold on average, a situation for which we reserve 
the acronym STCI (statistical topological crystalline in¬ 
sulator). Thereto, we find protected boundary states 
in five out of ten AZ symmetry classes by demanding 
that it also obeys a nonlocal reflection symmetry on av¬ 
erage. We have shown that the stability against disorder 
of (d — l)-dimensional boundary states of d-dimensional 
TCIs is understood from the phase diagram of (d — 1)- 
dimensional Dirac Hamiltonians with random mass. We 
have also pointed out the existence of an even-odd depen¬ 
dence of the stability on the parameter N that specifies 
the matrix dimension of Dirac Hamiltonians. Finally, we 
have shown that the alloy Sn-^^Pb^Te is an example of 
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a STCI in the 3D symmetry class All with an average 
reflection symmetry such that its surface states are pro¬ 
tected by the zeroth homotopy group tt 0 (V) = Z 2 for 
surface normalized Dirac masses. 

Before closing this section, we comment on two limita¬ 
tions of our approach. First, the random Dirac Hamilto¬ 
nians that we have studied are low-energy effective mod¬ 
els that should be derived from tight-binding models on 
regular lattices. In the context of topological insulators, 
Dirac Hamiltonians are only good effective models at en¬ 
ergy scales below the band gap. Our phase diagrams for 
Dirac Hamiltonians with random masses indicate that 
the metallic phases are stable even in the limit of strong 
disorder, whereas strong disorder should generically re¬ 
sult in localization in lattice models. There is no contra¬ 
diction here, since there is an implicit assumption in our 
analysis of Dirac Hamiltonians that the disorder strength 
is always smaller than the energy scale below which Dirac 
Hamiltonians are valid approximations to the underlying 
lattice models. We note that several attempts were pub¬ 
lished to study within a NLSM approach the localization 
transitions in the strong-disorder regime (of lat tice mod¬ 
els) in the chiral and All symmetry classes EUHliSSI 

Second, our approach was based on K-theory. K- 
theory aims at identifying stable homotopy groups of 
topological spaces. From the perspective of band theory, 
stability is to be interpreted as a topological attribute 
that is robust against any addition of trivial bands. Thus, 
there are cases to which our theory cannot be applied. An 
example is the localization problem of a band insulator 
characterized by a nontrivial Hopf map. A band insu¬ 
lator characterized by a Hopf map- 28 131 132 [a nontrivial 
element of 7r 3 (S' 2 )] is unstable to the addition of trivial 
bands. Hence such a Hopf map is not expected to be 
captured by the zeroth homotopy groups of the topolog¬ 
ical spaces of normalized Dirac masses 7r 0 (V). Indeed, 
one cannot write down a Dirac Hamiltonian with a con¬ 
stant mass term supporting a non-trivial Hopf map, i.e., 
no Dirac Hamiltonian with a constant mass term corre¬ 
sponds to a band insulator with a non-trivial Hopf map. 
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Appendix A: Explicit constructions of the compact 
topological spaces V in ID 


The properties of Sec. [TT] are central to this paper. 
They are applications of deep results in algebraic topol¬ 
ogy. They can be verified through explicit constructions 
of massive Dirac Hamiltonians for any given dimension d 
of space, as we are going to do when d = 1 in Secs. HH 
and |A 2| We shall then motivate why the properties of 
Sec. |II| hold in any d = 1,2,... in Sec. |A 3| We close 


Appendix [A] with a discussion of massless Dirac Hamil¬ 
tonians in Sec. ro 


1. Examples of Dirac Hamiltonian of rank 2 


For simplicity, we represent a generic Dirac Hamilto¬ 
nian of rank r = 2 that obeys periodic boundary (ring 
geometry) conditions by 

'H(k) := r 3 k + r 3 A x + r 2 M 2 + r 3 + t 0 A 0 . (Al) 


The Fermi velocity and the Planck constants have been 
set to unity. The unit 2x2 matrix is denoted by r 0 
and t = (ti,t 2 ,t 3 ) are three Pauli matrices. The real 
number A x couples as a vector potential induced by a 
magnetic flux at the center of a ring would do. The 
real number A 0 couples as a scalar potential would do, 
i.e., as the chemical potential. The real parameters 
and M 2 couple to the two available Pauli matrices that 
anticommute with the kinetic energy r 3 fc. The four real 
parameters A x , M 1 , M 2 , and A 0 furnish an exhaustive 
parametrization of a 2 x 2 Hermitian matrix. 


When A x = A 0 = 0, the Dirac Hamiltonian (Al) has 
the eigenvalues 

e(fc) = ± \Jk 2 + M 2 , M 2 := M\ + M|. (A2) 


The real parameters M ± and M 2 thus parametrize a Dirac 
mass. More precisely, we may write 


t 2 M 2 + Ti M x =Mf3(0) (A3a) 

where 

M := \/M? + M?, 9 := arctan 

V 1 2 Mi (A3b) 

/ 3(9 ) := r 3 cos# + r 2 sin0, 

in which case (a = r 3 ) 

(3 2 (9) = 1, {/3(0),a} = O. (A4) 

Thus, the normalized Dirac mass f3{9) is parametrized by 
the angle 0 < 9 < 2ir. As a topological set, the normal¬ 
ized Dirac masses {/?(0)|O < 9 < 27t} are homeomorphic 
to the circle S 1 . They are also homeomorphic to E/(l) 
through the map 

/? 12 :[O,27rKH(l),0^[/3(0)] 12 . (A5) 
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Symmetry class A: The Dirac Hamiltonian (Al) 


with A x , M-|. M 2 , A 0 smooth and non-vanishing func¬ 
tions of the position x along the ring is said to belong to 
the AZ symmetry class A. As we have seen, the set 


V£i. r=2 :={/3(0)|O<0<2 tt} =:S 1 


(A6) 


a circle, is the topological space of normalized Dirac 
masses associated with the Dirac Hamiltonian by adding 
to the kinetic contribution with the Dirac matrix a the 
normalized mass matrix (3(0) for 0 < 0 < 2ir. We note 
that V^ =1 r=2 and U( 1) are homeomorphic as topologi¬ 
cal spaces. Thus, they share the same honrotopy groups. 
On the other hand, V^L X r=2 is not a group under matrix 
multiplication while [7(1) is. 

In addition to the conservation of the fermion num¬ 


ber, we may impose TRS on the Dirac Hamiltonian (Al) 
There are two possibilities to do so. 

Symmetry class All: If charge conservation holds 
and TRS is imposed through 


H(k) = +t 2 H*(—k) r 2 , 


(A7a) 


The topological spaces V^ 1 r=2 and XJ( 1) are homeo¬ 
morphic. (This homeomorphism is not a group homo¬ 
morphism, for y d i lr=2 is not a group while [7(1) is.) 
Consequently, they share the same homotopy groups. 

The standard symmetry classes A, All, and Al can be 
further constrained by imposing the CHS. This gives the 
following three possibilities. 

Symmetry class AIII: If charge conservation holds 


together with the CHS 


H(k) = — Tj H(k) t x , 

(A9a) 

then 


kL(k) =T 3 k + T 3 A 1 +r 2 M 2 . 

(A9b) 


There is a unique mass matrix. The topological space of 
normalized Dirac masses obtained by adding to the Dirac 
kinetic contribution a mass matrix squaring to unity and 
obeying the CHS is 

2 = {± r 2 }- (A9 c) 


then 


n(k) =r 3 k + t 0 A 0 


(A7b) 


Symmetry class CII: It is not possible to write down 
a 2 x 2 Dirac equation in the symmetry class CII. For 
example, imposing 


No mass matrix is permissible if TRS squares to minus 
the identity. The topological space of normalized Dirac 
masses in the symmetry class All is the empty set 


V AU - 

v d=l,r=2 — 


(A7c) 


Because of the fermion-doubling problem 11331 the only 
way to realize (A7b I as the low-energy and long wave¬ 


length limit of a lattice model is on the boundary of a 
two-dimensional topological insulator in the symmetry 
class AIL 

Symmetry class Al: If charge conservation holds and 
TRS is imposed through 


n(k) = +7i n*(-k) iq, (A8a) 

then 

/ H(k) = r 3 k + t 2 M 2 + t x M x + t 0 A 0 . (A8b) 


The same mass matrix as in the symmetry class A is 
permissible if TRS squares to the identity. The honre- 
omorphy between the allowed masses in the symmetry 
classes A and Al is accidental. It does not hold for larger 
representations of the Dirac matrix as we shall verify ex- 
plicity when considering a rank 4 Dirac matrix below. 
The topological space of normalized Dirac masses ob¬ 
tained by augmenting the Dirac kinetic contribution by 
adding a mass matrix squaring to unity and obeying this 
TRS is 


%(k) = —t x 'H(k) t x , %(k) = +r 9 H*(—k) r 2 , 

(A10) 

enforces the symmetry class Dill, for composing the CHS 
with the TRS delivers a PHS that squares to the unity 
and not minus the unity. In order to implement the sym¬ 
metry constraints of class CII, we need to consider a 4 x 4 
Dirac equation (see Sec. A 2). 


Symmetry class BDI: If charge conservation holds 
together with 


%(k) = —t x H(k) t x , 


H(k) = +t x H*(—k) r lt 

(Alla) 


then 


T~L(k) = r 3 k + t 2 M 2 . (Allb) 

There is a unique mass matrix. The topological space 
of normalized Dirac masses obtained by adding to the 
Dirac kinetic contribution a mass matrix squaring to 
unity while preserving TRS and PHS (a product of TRS 
and CHS), both of which square to unity, is 

V®-=a = i ±r 2i- ( Allc ) 

Now, we move on to the four Bogoliubov-de-Gennes 
symmetry classes with PHS. 

Symmetry class D: If we impose PHS through 

H(k) = - H*(-k ), (Al2a) 


then 


V^i, r = 2 “ im\0 <0<2n} =: S 1 . (A8c) 


'H(k) = t 3 k + t 2 M 2 . (Al2b) 






32 


There is a unique mass matrix. The topological space of 
normalized Dirac masses obtained by adding to the Dirac 
kinetic contribution a mass matrix squaring to unity and 
preserving the PHS squaring to unity is 

V£i,r =2 = {irj. (A12c) 

Symmetry class Dill: If we impose PHS and TRS 


through 



n(k) = -n*(-k), 

U{k) = +T 2 H*{-k)r 2 , 

(A13a) 

respectively, then 

'H(k) = t 3 k. 

(Al3b) 


No mass matrix is permissible if TRS squares to minus 
the identity. The topological space of normalized Dirac 
masses in the symmetry class Dill is the empty set 


worn 

l.r—2 


(A13c) 


Because of the fermion-doubling problemp^J the only 
way to realize (A13bl as the low-energy and long wave¬ 
length limit of a lattice model is on the boundary of a 
two-dimensional topological superconductor in the sym¬ 
metry class DHL 

Symmetry class C: If we impose PHS through 


H(k) = -r 2 H*{-k)T 2 , (A14a) 


then 


'H(k) = r 3 A 1 + t 2 M 2 + t x M x . (Al4b) 


PHS squaring to minus unity prohibits a kinetic energy 
in any Dirac Hamiltonian of rank 2 in the symmetry class 
C. 

Symmetry class Cl: If we impose PHS and TRS 
through 

TL(k) = — t 2 T-L*{—k ) r 2 , H(k) = +t x H*(—k) t x , 

(A15a) 

respectively, then 


'H{k) = t 2 M 2 + t x M x . (Al5b) 


PHS squaring to minus unity prohibits a kinetic energy 
in the symmetry class CL 


2. Examples of Dirac Hamiltonian of rank 4 


Eq. ( |Al| ). A second unit 2x2 matrix is denoted by er 0 
and cr = (oq, er 2 , cr 3 ) are another set of three Pauli matri¬ 
ces. The summation convention over the repeated index 
v = 0,1,2,3 is implied. There are four real-valued pa¬ 
rameters for the components A x v with v = 0,1, 2, 3 of an 
17(2) vector potential, eight for the components M x v and 
M 2 v with v = 0,1, 2,3 of two independent 17(2) masses, 
and four for the components A 0 v with v = 0,1, 2,3 of 
an 17(2) scalar potential. As it should be there are 16 
real-valued free parameters (functions if we opt to break 
translation invariance). 


Symmetry class A: The Dirac Hamiltonian (A16) 


with A l in M-| M 2v) Aq v smooth and non-vanishing 
functions of the position x along the ring is said to belong 
to the AZ symmetry class A. All eight mass matrices of 
rank r = 4 in the ID symmetry class A can be arranged 
into the four pairs (M 1 M 2 „) with v = 0,1, 2,3 of anti¬ 
commuting masses. The topological space of normalized 
Dirac masses obtained by adding to the kinetic energy a 
mass matrix squaring to the unit matrix is® 


V d t 1>r =4 ~\P = 



U £ 17(2) } . (Al7a) 


As a topological space, it is thus homeomorphic to 
1/(2) ~ 17(1) x 517(2) ~ S 1 x 5 3 , an interpretation ren¬ 
dered plausible by the parametrization 

V d t hr=4 = {M ■ X+N ■ Y\M 2 = cos 2 6, N = tan 6M) 



=: S 1 x 5 3 , 


M 

:= {M 2 Q,M x l ,M 12 ,M X 3 j , 


N 

:= (— M x M 21 , M 2 2 , M 2 3 j , 


X 

:= (r 2 ® ct 0 , Tj <g> oq, -rq <g> cr 2 , t x 

(g)cr 3 ) , 

Y : 

= (-Tj <g> CT 0 , t 2 ® <7 1; t 2 ® a 2 , t 2 <g) cr 3 ) , 


(Al7b) 


We note that V^L X r=4 is not closed under matrix multi¬ 
plication so that it does not carry the group structure of 

U( 2). 

In addition to the conservation of the fermion number, 
we may impose TRS on the Dirac Hamiltonian (A16) 
There are two possibilities to do so. 

Symmetry class All: If charge conservation holds 
together with TRS through 


H(k) = +t x (g> a 2 %*{— k) t x 0 cr 2 , (Al8a) 


then 


Consider the generic Dirac Hamiltonian of rank r = 4 
that obeys periodic boundary conditions (ring geometry) 


n(k) :=t 3 <8> cr 0 k + t 3 <g> cq, A X v +r 2 ®a v M 2 v 
+ T x ®a u M l v + r 0 ® t t v A 0 v . 


(A16) 


The Fermi velocity and the Planck constants have been 
set to unity. The matrices r 0 and t were defined in 


U(k) =T 3 ®a 0 k+ Y2 T 3® a v A i,» + t 2 ®(t 0 M 2 0 

i/=1,2,3 

+ T x 0a o M 10 + r 0 O ct 0 A 0 0 . 

(A18b) 

Observe that by doubling the Dirac Hamiltonian ( A7b[ ) , 
we went from no mass matrix to two anticommuting 
mass matrices. The topological space of normalized Dirac 









33 


masses obtained by adding to the Dirac kinetic contribu¬ 
tion a mass matrix squaring to unity and obeying this 
TRSis obtained by imposing the TRS (A18a) on Eq. 
(|Al7a|) to be 


K An 


d—l,r=4 


o u 

0, 


17 = +a 2 U t <j 2 £ 17 (2) 


(A18c) 

As a topological space, V£p[ r=4 can be shown to be 
homeomorphic to U (2) / Sp{\) ~ U( 1) x SU(2)/SU(2) ~ 
17(1), an interpretation rendered plausible by the 
parametrization 


V d A 2lr =4 = {M-X\M 2 = l}=:S\ 
M := (M 2 0 ,M 10 ) , 

X := (r 2 g a 0 ,T x (g) cr 0 ). 


(Al8d) 


Symmetry class AI: If charge conservation holds to¬ 
gether with TRS through 


H{k) = +t x (gi <T 0 H*(—k) t x (g cr 0 , 


(A19a) 


then 


T~L{k) = r 3 g <r 0 fc + r 3 g <x 2 di j2 + ^ (r 2 g er„ M 2 i , 

i/=0,l,3 

+ g cr y M 1>t , + T 0 (g A 0> „). 

(A19b) 


Symmetry class AIII: If charge conservation holds 
together with the CHS 


H(k) = —t 4 (g (T 0 'H(k) t x (g CT 0 , 


(A20a) 


then 


1 ° v M 2 ,v) 


H(k) = r 3 ®a 0 k+ ^2 (r 3 <g a v A l v + r 2 

y=0,l,2,3 

(A20b) 

The Dirac mass matrix r 2 g cr 0 M 2 0 that descends from 
Eq. (A9b) commutes with the triplet of anticommuting 


mass matrices r 2 gcr 1 M 2 4 , r 2 gcj 2 M 2 2 , and t 2 ®<j 3 M 2 3 . 
The topological space of normalized Dirac masses ob¬ 
tained by adding to the Dirac kinetic contribution a mass 
matrix squaring to the unit matrix and obeying the CHS 


v- 


AIII 
d—l,r=4 ‘ 


:= <P = r 2 


A : = UI n U\ 


m,n = 0,1,2, m + n = 2, U £ 17(2), 
m times n times N 


I m,n ■= diag(—1,..., —1, +1,..., +1)j. 

(A20c) 

As a topological space, it is thus homeomorphic to 
17(2)/[t7(2) x 17(0)]U17(2)/[17(1) x 1/(1)] U 17(2)/[17(0) x 
17(2)], as is also apparent from the parametrization 


There are six mass matrices of rank r = 4 in the ID sym¬ 
metry class AI that can be arranged into the three pairs 
(M-, v , M 2 v ) with v = 0,1,3 of anticommuting masses. 
The topological space of normalized Dirac masses ob¬ 
tained by adding to the Dirac kinetic contribution a mass 
matrix squaring to unity and obeying this TRS is 


^=l,r=4 “ < P = 



U = +U T £ 17(2) 


(A19c) 


As a topological space, V ( ^ 1 r=4 can be shown to 
be homeomorphic to U(2)/0(2) ~ 17(1)/{±1} x 

SU(2)/U(1) ~ S 1 x S 2 , an interpretation rendered plau¬ 
sible by the parametrization 

17 d AI i,r=4 = {M • X + N ■ Y\M 2 = cos 2 9,N = tan 0M) 
=: S 1 x S 2 , 

M := (^M 2 0 ,M 11 ,M 13 ) , 

N := (— M lf0 , M 21 , M 2 3 ) , 

X := (r 2 g cr 0 , t x g ct 4 , t 4 g ct 3 ) , 

Y := (—r 4 g fi 0 , t 2 g O’!, r 2 g cr 3 ) . 

(A19d) 


The standard symmetry classes A, All, and AI can be 
further constrained by imposing the CHS. This gives the 
following three possibilities. 


V^lr =4 = {±r 2 g cr 0 } U {M • X\M 2 = 1} , 

M := (M 2j i, M 2>2 , M 23 ) , (A20d) 

X := (r 2 g cr 4 , r 2 g cr 2 , r 2 g a 3 ) , 

(recall that S' 2 ~ SU(2)/U(\) so that t/(2)/[17(l) x 
17(1)] ~ S 2 ). 

Symmetry class CII: If charge conservation holds 
together with CHS and TRS 


’H(k) = — r 4 g cr 0 'H(k) t x g ct 0 , 
?!(£;) = +t 4 g cr 2 T~L*{—k) r 4 g cr 2 , 


(A21a) 

(A21b) 


respectively, then 


^(*0 = ^3 


>cr 0 k 


1 -^2,0- 


r 3 ® + r 2 

^=1,2,3 

(A21c) 

There is a unique mass matrix, as was the case in 
Eqs. (A9b) and (Alibi. The topological space of nor¬ 


malized Dirac masses obtained by adding to the Dirac 
kinetic contribution a mass matrix squaring to the unit 
matrix and obeying CHS and TRS squaring to minus 
unity is 


yen 
v d.=\,r=i 


■= [P S 17 d t™ r=4 P = (r 4 g cr 2 ) P* (r 4 g ct 2 )}. 


(A21d) 


As a topological space, Hjr 1 ] r=4 can be shown to be 
homeomorphic to Sp(l) / Sp(lj x Sp(0) U Sp(l)/Sp(0) x 
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Sp(l), an interpretation rendered plausible by the 
parametrization 

V£L l h r=4 = { ±t 2 ® cto}- (A21e) 

Symmetry class BDI: If charge conservation holds 
together together with CHS and TRS 

H(k) = — t x 0 a 0 H(k) t x 0 cr 0 , (A22a) 

H(k) — +r t 0 a 0 H*(—k) r 4 0 cr 0 , (A22b) 

respectively, then 


As a topological space, V^ =1 r=4 is homeomorphic to 
0(2), as Vf =lr=A -^US 1 ~ 1/(1) x Z 2 ~ 0(2). 

Observe that V'jl 1 r=4 is not closed under matrix mul¬ 
tiplication, i.e., it does not carry the group structure of 
0 ( 2 ). 

Symmetry class Dill: If we impose PHS and TRS 
through 

H(k) = —H*(—k), H(k) = +r 2 0 cr 0 H*(—k) t 2 0 <r 0 , 

(A25a) 

then 


7~L(k) = T 3 0cr o /c + T 3 0 cr 2 A 12 + 55 r 2 0 cr„ M 2 ; ,. 

i/=0,l,3 

(A22c) 

The Dirac mass matrix t 2 0 a 0 M 2 0 that descends from 
Eq. (Allb) commutes with the pair of anticommuting 
mass matrices t 2 ®o 1 M 2 and t 2 0ct 3 M 2 3 . The topolog¬ 
ical space of normalized Dirac masses obtained by adding 
to the Dirac kinetic contribution a mass matrix squaring 
to the unit matrix and obeying CHS and TRS squaring 
to unity is 


T/BDI 

v d= l,r=4 ■- 


{P G V ( tl\] r=i 


P = ( T 1 0 a o) P* ( T 1 ® cr 0 )}- 


(A22d) 


As a topological space, Vf™ r=A can be shown to be 
homeomorphic to 0(2)/[0(2) x 0(0)] U 0(2)/[0(l) x 
0(1)] U 0(2)/[0(0) x 0(2)], an interpretation rendered 
plausible from the parametrization 


K?= D iU = {±^2 « cr 0 } U {M • X\M 2 = 1} , 

M := (Af 2jl , M 2 3 ) , (A22e) 

X := (t 2 0(7 1 ,t 2 0ct 3 ), 


(recall that S 1 ~ 0(2)/[0(l) x 0(1)]). 

Now, we move on to the four BdG symmetry classes. 
Symmetry class D: If we impose PHS through 


n(k) = - n*{-k ), (A23a) 

then 

H(k) = r 3 0 a 0 k + t 3 0 a 2 A 12 + ^ t 2 0 a v M 2 v 

w=0,l,3 

+ H 0 cr 2 M 12 + r 0 0 cr 2 A 0i2 . 

(A23b) 


There are four Dirac mass matrices. None commutes 
with all remaining ones. However, each of them is anti¬ 
symmetric and so is their sum. The topological space of 
normalized Dirac masses obtained by adding to the Dirac 
kinetic contribution a mass matrix squaring to the unit 
matrix and obeying PHS squaring to unity is 

V[L hr =i = {M ■ X\ M 2 = 1} U {N ■ Y\N 2 = 1} , 

M := (M 21 ,M 2 S ) , N := {M 2 0 ,M 12 ) > 

X := (r 2 0cr 1 ,r 2 0fT 3 ), Y := (r 2 0 <t 0 ,t 4 0 a 2 ). 

(A24) 


T~L(k) = t 3 0 <Tq k + r 3 0 cr 2 A ± 2 + t x 0 cr 2 M x 2 . (A25b) 


Observe that there is only one Dirac mass matrix [there 
was none in Eq. (A13b)]. Moreover, this Dirac mass ma¬ 
trix is Hermitian and antisymmetric. The topological 
space of normalized Dirac masses obtained by adding to 
the Dirac kinetic contribution a mass matrix squaring to 
the unit matrix and obeying this PHS and this TRS is 


Vd= l,r=4 = {±7i ®<r 2 }. 


(A25c) 


As a topological space, V£Jl I r=4 is homeomorphic to 
0(2)/J7(l). 

Symmetry class C: If we impose PHS through 


■H(fc) =—r 0 0 cr 2 ■?!*(—fc) r 0 0 (t 2 , (A26a) 


then 

'H(k) = t 3 0 cr 0 k + ^ r 3 ® + T 2 ®a 0 M 2 0 

i/=1,2,3 

+ 51 T 1®<Tv M 1,v+ 55 T V®< J v A V,v 

i/=l,2,3 i/= 1,2,3 

(A26b) 

There are four mass matrices that anticommute pairwise. 
The topological space of normalized Dirac masses ob¬ 
tained by adding to the Dirac kinetic contribution a mass 
matrix squaring to the unit matrix and obeying this PHS 
is 


H £ 1 , r=4 = {M-X|M 2 = l} =: S 3 , 

M := (M 2 q, M 4 4 , M 12 , M 13 ) , (A26c) 

X := (t 2 0 cr 0 , t 4 0 a 11 r 4 0 <7 2 , r 4 0 cr 3 ). 

As a topological space, V^ =[ r=4 is homeomorphic to 
Sp( 1) since we have Sp( 1) ~ SU( 2) ~ S' 3 . Observe that 
KEl,r=4 is not closed under matrix multiplication, i.e., 
it does not carry the group structure of Sp(l). 

Symmetry class Cl: If we impose PHS and TRS 
through 


T~L(k) = —t 0 0 a 2 %*(—k) r 0 0 a 2 , 
'H(k) = +r 4 0 erg T~L*{—k) r 4 0 u 0 , 


(A27a) 

(A27b) 





35 


respectively, then 

'H(k) = t 3 g cr 0 fc + r 3 g cr 2 A 12 + r 2 g a 0 M 20 

+ Y ( T i ® M l,v + t 0 g o v A 0 R . ( A27c ) 

i/=l,3 

There are three mass matrices that anticommute pair¬ 
wise. The topological space of normalized Dirac masses 
obtained by adding to the Dirac kinetic contribution a 
mass matrix squaring to the unit matrix and obeying 
this PHS and this TRS is 

V d c Yr =4 = {M, ■ x t \ Mf = 1} =: S 2 , 

M := (^M 2 0 ,M 11 ,M 1 ^ , (A28) 

X := (t 2 g) ct 0 , Tj_ g a 1 , t x g <r 3 ). 

As a topological space, V d l x r=4 is homeomorphic 
to Sp{l)/U{l) since we have the homeomorphism 
Sp(l)/U(l) ~ SU(2)/U{1) ~ S 2 . 

3. Interpretation of 7r 0 (H) = Z,Z 2 ,0 

Assume that we have chosen the dimensionality d of 
space and the AZ symmetry class. We may then increase 
the rank r of the Dirac matrices until we reach the small¬ 
est rank r min for which we may write the massive Dirac 
Hamiltonian 


y = r\ - -L 777 R . -L 

' *'min min ^min 1 * * * : 


(A29) 


where the (d + 1) matrices a min v ..., a min d , /3 min an¬ 
ticommute pairwise and square to the unit r min x r min 
matrix, to is a real-valued (mass) parameter, and ..” 
accounts for all scalar and vector gauge contributions as 
well as for the possibility of additional mass terms. Let 
V dr be the compact topological space associated with 
the normalized Dirac masses that enter on the right-hand 
side of Eq. (A29). There are then three possibilities when 


^oi^d.r . ) A 0. 


a. Case n 0 (V) A 0 


Assume that no other normalized r n 
trix than /3 min enters the right-hand side of Eq. ( |A29| . It 
then follows that 


matrices that have become available in the given AZ sym¬ 
metry class upon increasing r from from r min to r min N 
with N = 2,3,.... We may then write a generic Dirac 
mass term as 


/3 m in g M= /3 min 



(A31) 


where the N x N Hermitian matrix M squares to the 
unit N x N matrix U^, the N x N Hermitian matrix 
A is traceless, both M and A are restricted by the AZ 
symmetry class, and v is the topological index defined in 
Eq. (3.8). The mass matrix /3 min g 11 N is unique up to a 


sign. The matrix M must obey the polar decomposition 


M = U 


U, 


t 


N—n,n ± N—n,n N—n,n ’ 
IN—n,n n't -^n)> 


(A32) 


where the matrices U N _ form either U(N),0(N), or 
Sp(N) depending on the AZ symmetry class. (For the 
symplectic case, a matrix in Sp(N) has qu aternions as 
matrix elements.) The representation (A32) is unique up 
to multiplication of U N _ n n from the right by the block 
diagonal matrix diag(/7 JV _ n , U n ), where the nxn matrix 
U n is taken from either U(n),0(n ), or Sp(n) d epen ding 
on the AZ symmetry class. The representations ( |3T| and 
(|3.8|) follow. 


Case of first descendant 7T 0 (V d = R ± ) = Z 2 : There 
is a pair of two anticommuting normalized Dirac mass 
matrices that span V rl 2 upon increasing r from r ■ 
to 2r min . We can choose a representation of the Dirac 
mass matrix as follows. If AT = 2, 3,... and if we define 
the Hermitian r min /2 x r min /2 matrix p min by 


Anil) = : Pn 


> t 2, 


we may then write 

Pmin ® M = p min g (r 2 A 2 + r x A x ) 


(A33a) 


(A33b) 


for the generic normalized Dirac mass. Here, the 2Nx2N 
matrix M is Hermitian and antisymmetric because the 
N x N matrix A 2 is Hermitian and symmetric (A 2 com¬ 
mutes with the operation K for complex conjugation), 
while the N x N matrix A 1 is Hermitian and antisym¬ 
metric (A x anticommutes with the operation K for com¬ 
plex conjugation). The matrix M is a solution of the 
extension problem 


CL 2 = {ir 3 ; K, iK} -> C7, 3 = {ir 3 ; K, iK, M}. (A34) 


V d , r . = {±/3 min }. 


(A30) 


Hence the ground state of T~L min for to > 0 cannot be 
smoothly deformed into the ground state of H min for to < 
0 without closing the gap proportional to \m\. There are 
then three possible outcomes upon increasing the rank in 

Eq. ( A29| ) from r min to r min N with N = 2, 3,_ 

Case 7r 0 (V^) = Z: There is one normalized Dirac mass 
matrix which commutes with all normalized Dirac mass 


Indeed, it is the matrix ir 3 that enters in the kinetic 
contribution to the Dirac Hamiltonian in the exam¬ 


ple (A23b) for a one-dimensional Dirac Hamiltonian in 


the symmetry class D. The explicit representation of 
Eq. (3.151 for this example thus follows from solving with 
Eq. (A33) the extension problem (A34). 


It is instructive to compare E qs. (A33a ) and (A33b ) 
with the explicit representations (A12b| and (A23b) cor¬ 


responding to the one-dimensional symmetry class D 






























36 


with r = r min and r = 2r min , respectively. This com¬ 
parison suggests that the choice A 2 = ± 11 N and A x = 0 
is special for any odd N. Indeed, one may verify that 
it is impossible to construct an antisymmetric matrix A x 
of odd rank N that connects smoothly, without closing a 
gap, between the choices A 2 = +1^ and A x = 0 on the 
one hand, and A 2 = — and A x = 0 on the other hand. 
Hence the choices A 2 = ±t N and A 1 = 0 represent two 
topologically distinct phases when TV is odd. 

Finally, we observe that there is no unique Dirac mass 
matrix that commutes with all other Dirac mass matrices 
for IV > 1. For example, the Dirac mass matrix with 
A 2 = and A 1 = 0 neither commutes nor anticommutes 
with the Dirac mass matrix for which 


/ COS0 

0 

°1 

1 ° _i 

°\ 

0 

II 

CN 

COS0 

0 

, A 3 = sind +i 0 

° 

\ 0 

0 


O 

O 

0/ 


(A35) 

Case of second descendant (V d ~ A) — ^2 : 
There is a pair of three anticommuting normalized Dirac 
mass matrices that span V d 2r upon increasing r from 
r min to 2r min . We can choose a representation of the 
Dirac mass matrix as follows. If N = 2, 3,... and if we 
define the Hermitian r min /2 x r min /2 matrix p min by 

Anin = : Pmin® ^2, (A36a) 


we may then write 

Pmin ® M = Pmin ® ( r 2 A + T 3 A 3 + T 1 A 1 + T 0 A 0 ) 

(A36b) 

for the generic Dirac mass matrix. Here, the 2 N x 2 N 
matrix M is Hermitian and antisymmetric because the 
N x N matrix A 2 is Hermitian and symmetric (A 2 com¬ 
mutes with the operation K for complex conjugation), 
while the N x N matrix A 3 , A ±1 and A 0 are Hermitian 
and antisymmetric ( A 3 and A 3 anticommute with the op¬ 
eration K for complex conjugation). The matrix M is a 
solution of the extension problem 


then 


T-L(k) = t 3 (g cr 0 <g> p 0 k 

+ A 1,2,v+°v ®P 2 A l,va) 

i/=0,l,3 

+ T l®( a 2®P„ M l,2,v + Vv®P2 M l,v,2)- 

i/=0,l,3 

(A38c) 

The topological space of normalized Dirac masses ob¬ 
tained by adding to the Dirac kinetic contribution a mass 
matrix squaring to the unit matrix and obeying PHS 
squaring to unity is 

V d Ji, r =8 = { M ' AM 2 = 1} U {N ■ Y\N 2 = 1} , 

M := (M 12 0 ,M 112 ,M 13 2 ) , 

N := {M l Q 2 ,M 12 l ,M 12 3 j , 

X := (iq <g> cr 2 <g> p 0 , t x ® <7 1 ® p 2 , T 3 ^a 3 ® p 2 ) , 

Y := (t x <g) cr 0 ® p 2 , t 3 ® a 2 (g p ± , t x <S> a 2 ® p 3 ) . 

(A39) 

As a topological space, V^JJ I r _ 8 ~ S 2 US 2 is homeomor- 
phic to 0(4)//7(2), for the homeomorphisms 

0(4)/Z7(2) ~Z 2 x SO(4)/[U(l) x SU{2 )] 

~Z 2 x [50(3) x SU{2))/[SO(2) x SU{ 2)] 
~Z 2 x SO(3)/SO(2) 

~S 2 US 2 (A40) 


hold between topological spaces. I n ord er to com pare th e 
Dirac mass matrices entering Eq. (A381 with Eq. (A36b I, 
noti ce fir st that p min = t x and s econd that a v in 
Eq. (A39) plays the role of t v in Eq. (A36b I. 


b. Case n 0 (V) = 0 


There is at least one r • x r 


Hermitian ma¬ 


trix /3 ' min that enters “..in Eq. (A29) such that 
it anticommutes with any one of the (d + 1) matri¬ 
ces a min x ,..., <a mind , /3 min and it squares to the unit 
r min x r min matrix. For any 0 < 9 < 2ir , we may then 
define the normalized mass matrix 


Cl 0}2 = {; K, iK} -> Cl 0t 3 = {; K, iK, M}, (A37) 


which constrains M to be pure imaginary, as can be ver¬ 
ified for Eq. (A36b). The representation (3.15) follows. 


The normalized Dirac mass matrix (A36a) should be 


compared with that in Eq. (A251 for the one-dimensional 


symmetry class Dill with r = 4, and Eq. (A36b) should 


be compared with the normalized Dirac masses for the 
one-dimensional symmetry class Dill with r = 8 that we 
now present. If we impose PHS and TRS through 


H{k) = -n*(-k), (A38a) 

H(k) = +t 2 (g) <t 0 ® p 0 %*(—&) r 2 (g) a 0 g> p 0 , (A38b) 


AninW : = cos 9 f3 min + sin 0 (3' mh] 


(A41) 


that provides a smooth path between /3 min and 
Since /3 min can be chosen arbitrarily in the compact 
topological space V associated with the massive Dirac 
Hamiltonian (A29), we can rule out the existence of 
the disconnected subspace {±/3 •,} in lA . For gen- 
eral values of N = 1,2,3,..., the trivial zeroth ho- 
motopy group 7r 0 (Vb ) = 0 indicates that for any 

’ min ' 


given pair (/?,/?') € 
ily anticommuting), 
muting pairs £ V d 


V(j |r jv x V d ,r minN (not necessar- 
we find the sequence of anticom- 


N X Vd. 




d,r : 


. N X V dr . N , 

min ’ 


(l3 n ,F)eV B 


occurs in Eq. (A26). 


m , n JV> 

d,r . N X V d r N 

’ min ’ min 


(A.) $ 2 ) ^ 

as 
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4. Existence condition for the Dirac masses 


(A29). The question 


Assume that we are given Eq. 
posed and solved in Sec. |IT|was: 

Qi What is the compact topological space spanned by 
the normalized Dirac masses of rank r = r min N in 
the limit N —> oo? 


By definition, the Dirac Hamiltonian (A29) realizes a 


topological insulator if the compact topological space V 
spanned by the normalized Dirac mass /? min consists of 
no more than the set {±/3 min }. It then follows that, on 
any (d— l)-dimensional boundary of d-dimensional space, 
the Dirac Hamiltonian (A29) reduces to a massless Dirac 
Hamiltonian of the form 


H = a 


d 

idx 


(A42) 


where the (d— 1) matrices a ,..., a anticommute pairwise 
and square to the unit, their rank is r = r min /2, and 
accounts for scalar and vector gauge contributions, but 
does not contain Dirac mass terms. Moreover, by 
the definition of the minimal rank r min , no Dirac mass 
is permissible upon lowering the rank of the Dirac ma¬ 
trices entering the Dirac Hamiltonian (A29) holding the 


dimensionality of space and the AZ symmetry class fixed. 
These observations motivate the following question: 

Q2 Given an AZ symmetry class and given the dimen¬ 
sionality d of space, is a Dirac mass of rank r min 
permissible or not? 

Here, r min is the minimal rank for which it is possible 
to write down the kinetic contribution —i<5 • d/dx to the 
Dirac Hamiltonian given the dimensionality d of space 
and the AZ symmetry class. 

We want to know if a massless Dirac Hamiltonian be¬ 
longing to some prescribed d-dimensional AZ symmetry 
class accommodates a Dirac mass. To answer this ques¬ 


tion, we recall that Eqs. (2.15) and (2.17) define for any 


AZ symmetry class and for any dimension d a Clifford al¬ 
gebra that supports at least one normalized Dirac mass, 
provided the rank r of the Dirac matrices is no less than 


the minimal rank r min . If we replace r min by r n 


the 


answer to question Q2 is given by the extension problems 
(the values of p and q are taken from Table |l| 


Cl q -r = {e 


1? ' 


1> ' 


-i} 


q— 1 ? °q 


} = CL 


(A43a) 


for the complex symmetry classes A and AIII; 


Cl p , q - 1 — {e 1; ... 


Op 7 ®p+ 1) ■ • • ! C p+9 

p+17 • * • 7 ^p+q— 17 ^ p+q 


q-l) 


{ e l 7 • • • 7 6 p 7 6 , 


} = Clp, q 

(A43b) 

for the real symmetry classes BDI, D, Dill, CII, C and 
Cl; and 


C^p-1,9 — {f 


1 7 ' 


9-11 


{ei, 


p-li ‘ > °p+g 

e p— 1 1 e p> e p+l ? * * * ? c p+<? 


} = ci p , q 

(A43c) 


for the real symmetry classes AI and AIL In other words, 
for any given AZ symmetry class, we seek to enumer¬ 
ate all the distinct ways there are to construct the cor¬ 
responding Clifford algebra from Eqs. (2.18) and (2.19) 


out of one in the same symmetry class but with one fewer 
generator. Recall here that the Clifford algebras (2.18) 


and (2.19) are obtained after removing the generator as¬ 


sociated with the normalized Dirac mass in each of the 
Cliffo rd al gebras ( 2.15[ ) and ( 2.17| ). The solutions to 
Eqs. (A43) are the classifying spaces (the values of p and 
q are taken from Table [l| 


V = C, 


< 2-1 > 


V = R, 


q-p-li 


V = R, 


'P q~\~ 17 


(A44) 


respectively. On the one hand, if 7r 0 (I?) = Z,Z 2 , then 
the generator e p+q from the corresponding Clifford alge¬ 
bra is unique in that no additional generator exists that 
anticommutes with e p+q (see Sec. A3), i.e., no additional 
generator e p+q+1 that plays the role of a normalized Dirac 
mass represented by a Hermitian matrix of rank r min is 
allowed. On the other hand, if 7r 0 (W) = 0, then the gener¬ 
ator e p+q from the corresponding Clifford algebra is not 
unique in that there exist additional indepen dent gen¬ 
erators that anticommute with e p+q (see Sec. A3), i.e., 
a generator e p+q+1 that plays the role of a normalized 
Dirac mass represented by a Hermitian matrix of rank 
rmin i s allowed. 

For comparison, the answer to question Ql was given 
in Table [I] by the extension problems (the values of p and 
q are taken from Table |T| 


Cl q = { e n • • ■, e q } 

( e i7 ■ ■ • 7 e q , e q+1 } = Cl q+1 
for the complex symmetry classes A and AIII; 


(A45a) 


Clp,q {^17 • • • 7 Op 7 Cp+l7 • • • 7 Op+g} 

^ { e i7 • ■ • 7 Op 7 Op_|_l 7 • • • 7 Op+q 7 Op+q+1 } Cl p q +1 

(A45b) 


for the real symmetry classes BDI, D, Dill, CII, C and 
Cl; and 

Clp, q {oi? • ■ • 7 Op, Cp_j_i, ■ ■ ■, e p + q } 

t {c^, • ■ ■ 7 Op, Cp^_^, Cp_j_2 , • ■ • 7 6p_|_q_(_j} Clp+i q 

(A45c) 


for the real symmetry classes AI and AIL The solutions 
to Eqs. (A45) are the classifying spaces (the values of p 
and q are taken from Table [i]) 

V = C q , V = R q _ p , V = R p _ q+2 , (A46) 


respectively. Thus, the classifying space associated with 
Q2 in d dimensions coincides with the classifying space 
associated with Ql in d + 1 dimensions. The existence 
of nontrivial topological insulators in d + 1 dimensions 
and the gapless Dirac Hamiltonian with no allowed Dirac 
mass in d dimensions are equivalent. 

The final answer to the question Q2 is: 
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• No Dirac mass matrix of rank r min is permissible 
when 7r 0 (V) = Z. If TV channels are added by 
generalizing the Dirac matrices a to a <S> Ijy, then 
no Dirac mass of rank r min N is permissible when 
n 0 (V)=Z. 

• No Dirac mass matrix of rank r min is permissible 
when 7r 0 (V) = Z 2 . If N channels are added by 
generalizing the Dirac matrices a to a® Ijy, then 
no Dirac mass of rank r min N is permissible for odd 
N, while a Dirac mass of rank r min N is permissible 
for even N, when 7r 0 (V) = Z 2 . 

• A Dirac mass matrix of rank r min is permissible 
when 7r 0 (y) = 0. If AT channels are added by gen¬ 
eralizing the Dirac matrices a to a ® l^y, then a 
Dirac mass matrix of rank F min N is permissible 
when 7r 0 (F) = 0. 

Table IVII from Sec. IVIII follows. 


In order to compute the disorder-average over the cor¬ 
relations functions for the Grassmann-valued spinors, it 
is convenient to take the replica limit 

7 n — 1 

In Z = lim -. (B2) 

n-y o n 

If we reserve the index a, b = 1,..., n, we may write 
Z n = j V[tp]D[tp\ exp J d 2 r (L 0 + L int )j , (B3a) 

n 

L o :=i E $“ [ CT * + a v H^y)] ( B3b ) 

a= 1 

n 

A nt :=E|E^^ ( B3c ) 

i a,6=1 

for the disorder average over the nth power of the parti¬ 
tion function. The effective action 

S eS := J d 2 r (L 0 + L int ) (B3d) 


Appendix B: One-loop renormalization group 
analysis in 2D space 

We summarize the one-loop RG flows associated with 
the relevant 2D random Dirac Hamiltonians from Refs. nn 
and 11341 that we have used to deduce the global phase 
diagrams presented in Sec. [V] 

For simplicity, we present RG flows for vanishing mean 
values of all random potentials. A mean value for any 
random mass is a strongly relevant perturbation. This is 
the control parameter for the 2D phase diagrams along 
the horizontal axis. Although we do not present the con¬ 
tributions to the RG flows from non-vanishing mean val¬ 
ues of the random masses, their contributions are essen¬ 
tial to cross phase boundaries. They are thus implicitly 
assumed to be present. 


1. Replica limit 


descri bes an interacting and nonunitary quantum field 
theory!^! Disorder averaged correlation functions for the 
Grassmann fields are computed by first using the effective 
partition function (B3) and then by taking the replica 
limit (B21. 


2. One-loop renormalization group /3 functions 
from the operator product expansion 


Assume that space is the d-dimensional Euclidean 
space. Define the path integral 


Z : = j exp(-5' 0 - S int ), 

S o :=i / d d r^a- 
Sint : =E A « / 

a 


(B4a) 

(B4b) 

(B4c) 


Define the path integral 


Z:= J 

f V[ip}D[ip] exp(-S'), 

(Bla) 

S:= i 

J d 2 r (ipH 0 ip + ^2V L (p L {ip,ip) j , 

(Bib) 

U Q := 

\ L / 

Vx(-' l9 x) +<Jy{-ld y ). 

(Blc) 


Here, </> t (-)/>, ?/;) is a local quadratic form for the pair if) 
and ip of independent Grassmann-valued spinors. For 
any pair i, i', the real-valued functions V L and V L , are 
random with vanishing means and with 

K ( r ) V L , (r') = S L y g L 5(r - r') (Bid) 
as the only non-vanishing moments. 


where a are Dirac matrices that anticommute pairwise 
and square to the identity matrix, a is the short-distance 
cutoff, a (^J,ip ) is a local monomial of the Grassmann 
spinors, and A a is the corresponding dimensionless cou¬ 
pling. Assume the operator product expansion (OPE) 

W) = |r _ r ,|t+x»-* e SctO + ■ • • , (B4d) 

where the summation convention over repeated indices 
is implied, the tensor c abc is real valued, the scaling ex¬ 
ponent x a of 4> a (r) is real valued, and “...” contains 
all functions of |r — r'\ that are subleading relative to 
|r — r'\ x °' JrXb ~ x c in the limit r' —> r. The one-loop RG 
flows for the coupling constants is then 

-jj- = {d — x a ) X a + ~c abc X b X c , (B4e) 
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under the rescaling a —t a(l + dt?) of the short distance 
cutoff. When d = 2, <f> a is quartic in the Grassmann 
spinors, and, after performing the rescaling 

^a ' ^ 9ai 

we find the one-loop /3-functions 
a ^ d 9a 

Pa — 2 ^abc 9b 9c‘ 

3. Renormalization group flows for the AZ 
symmetry classes 

One-loop RG flows for random masse s an d random 
gauge potentials of 2D Dirac Hamiltonians—*!^ are used 
in this paper to decide whether a metallic phase or a crit¬ 
ical boundary separates topologically distinct insulating 
phases in the phase diagrams of Fig. [4] For the symme¬ 
try classes in which metallic conductivities acquire anti- 
weak localization corrections, a metallic phase is stable 
for large bare values of g. Whether a metallic phase per¬ 
sists from large to infinitesimally small bare values of g 
along the boundary between topologically distinct insu¬ 
lating phases is determined by the RG flows in the vicin¬ 
ity of g = 0 . If the characteristic disorder strength g 
is marginally relevant, then it is a metallic phase that 
separates two topologically distinct insulating phases for 
any g > 0 . (Evidently, we are assuming no intermediary 
fixed points between the ones at g = oo and g = 0 .) On 
the other hand, if the disorder strength g is marginally 
irrelevant, then it is a critical boundary that separates 
two topologically distinct insulating phases. (In this case 
there has to be a critical point where the phase boundary 
is terminated by a metallic phase.) 

For these reasons, we revisit one-loop RG flows of ran¬ 
dom masses and random gauge potentials for those sym¬ 
metry classes that are characterized by anti-weak local¬ 
ization corrections, i.e., the symmetry classes D, Dill, 
and AIL® The RG equations are obtained from OPEs of 
operators associated with all random masses and random 
gauge potentials allowed in each symmetry class. The 
OPEs are calculated by taking contr action s in a prod uct 
of two fermion quartic terms in Eqs. (B3c) and (B4c) us¬ 


ing the two-point correlation function of free Dirac fields, 

(B5) 


(f(r)#)) = ^«. 


2t r 


2 u ab~ 


To begin with, we first consider symmetry classes A and 
C and then move on to symmetry classes D, Dill, and 
AIL 


a. Symmetry class A 

Following Ref. [TS] we consider the 2D random Dirac 
Hamiltonian 

ft : = (~\d x + A x ) + G y (-i dy + A y ) +ma z + V 

(B6) 


of rank r = 2 in the symmetry class A. 

We seek the one-loop RG equations obeyed by the cou¬ 
pling constants 


(B4f) 

9 m = in 2 , 

(B7a) 


9 v = V 2 , 

(B7b) 

(B4g) 

9 a = Al + A*. 

(B7c) 


These couplings are the coefficients of the fermion quartic 
terms 


$m = - : V’ a ^ a V’ 6 cr 2 -i/’ { ':, (B8a) 

a,6=1 
n 

<f>y = - Y, (B8b) 

a,6=1 
n 

+ ■V^y'^ a ^ b <Jy^ b :), 


„b—1 


(B8c) 


where the interactions <& v , ( E> A , and 4> M are generated 
by taking the disorder average over the random scalar 
potential, the random vector potentials, and the random 
mass, respectively. Their non-vanishing OPEs ar d 106 f 13 ^ 

®M( r ) ^m(O) = 2 2 ^m( 0)> (B9a) 

7 i r z 

$ v (r)A> v (0) = J^A> v (0), (B9b) 

A> v {r) $ M (0) = —^ [$ v (0) + 2 $ a (0) - $ M (0)], 


7 T r * 


(B9c) 


$mM$4(0) = ^[$m( 0) + <M0)], (B9d) 

<MO <M0) = -!-j t<M0) + $ M (0)], (B9e) 

7 i r z 

where we have only kept singular parts proportional to 
r ~ 2 = r~ 2 on the right hand sides. Then, the RG equa¬ 
tions for symmetry class A read 


d 9 


M 


1 


^ ~ ~ {~ 9 m ~ 9 v 9 m + 9 v 9 a+ 9 a 9 m) i (B10a) 

= “ {. 9 v + 9 v 9 m + 9 v 9 a A 9 a 9 m) > (BlOb) 
d g A _ 2 

"d I~n 9v9M - 


(BlOc) 


As is well known, 15 there is an unstable line of critical 
points g v = g M = 0 and g A > 0 . However, for the 
generic case g v > 0 and g M > 0, the coupling con¬ 
stants always flow to the strong-coupling regime. If so, 
the one-loop RG equations are no longer applicable. The 
(marginally relevant) coupling constants are represented 
by the characteristic coupling g in the phase diagram of 
2D class A in Fig.[4|a). Since a stable metallic phase can¬ 
not exist in 2D systems of class A, topologically distinct 
insulating phases are separated by a phase boundary line 
on which wave functions are critical. 
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b. Symmetry class C 


with itself delivers the RG equatior l — 


Next, we consider the 2D random Dirac Hamiltonian 


H:=- i a x ®r 0 d x - \a y ® t 0 d y + m a z (g) t 0 

+ Ai ’i a i ® T J + Vi a o ® T i 

i=x,y j=x,y,z i=x,y,z 

(Blla) 

of rank r = 4 in the symmetry class C with the operation 
for charge conjugation represented by 

C := a x ® T y K, (Bllb) 

where K denotes complex conjugation. 

There are ten random functions allowed by the sym¬ 
metry constraints that multiply the matrices 


o- z ® r 0 , 


(mass) 

® T x> 

°x® T y, ff *®A 

(vector potentials) 

Vy®T x , 

Oy® T yl °y®T z , 

(vector potentials) (B12) 

0-q® t *. 

°0®T X , CT 0 (g)T x , 

(scalar potentials) 


respectively. Taking disorder average in the replicated 
action yields quartic interaction terms similar to those 
in Eqs. (B 81 . The one-loop RG equations for the quartic 


terms can be found in Refs. 11351 and [1361 Here we discuss 
only a minimal set of the RG equations that allow us to 
deduce the phase diagram in Fig. |4^a). Taking a x (g> r x 
and <7 y ®T x from the vector potentials and <r 0 (g>T a , from the 
scalar potential along with the mass a z ® r 0 among the 
above ten terms, we find that the one-loop RG equations 
for the coupling constants g v , g A , and g M take th e sam e 
form as those in the symmetry class A in Eq. (BlOl. 
Thus, the representative coupling constant g of disorder 
is marginally relevant as in the symmetry class A. This 
indicates that the phase diagram for the symmetry class 
C is qualitatively the same as that for the symmetry class 
A, as shown in Figs.[4ja) and (b). 


d 9m _ 1 2 

d£ 7T 9m ' 


(B15) 


Hence the variance g m of the random mass term is a 
marginally irrelevant coupling and the clean critical point 
g M = m = 0 is stable as long as the mean m = 0. 

The random Dirac Hamiltonian of rank r = 4 in the 
symmetry class D with the PHS symmetry generated by 
conjugation with C = <j x ® r 0 K depends on six random 
functions that multiply the matrices 


a z ® t 0 , & z ®t x , o z ®t z , (masses) 

a x <g> T y , o y ® T y , (vector potentials) (B16) 

<r 0 <g> Ty, (scalar potential) 


respectively. 

Our aim is to explain why, unlike for the case of r = 2 
in Eq. (|B13|), the massless clean critical point described 
by 


:= (“io* d x - i G y d y ) ® t 0 (B17) 

is (i) unstable to disorder, (ii) this instability depends 
sensitively on the choice of the probability distribution 
for the disorder, and (iii) we are led to the choice ( |3. 10 > 
for the probability distribution and to Fig. |4](d) to de¬ 
scribe qualitatively the generic phase diagram in the 2D 
symmetry class D with N = 2. 

We first choose the probability distribution for the six 
random functions such that they are independent and 
Gaussian distributed with the means and variances m € 
R and g 0 > 0 for the mass matrix a. (g> r 0 = /3 0 ; m x € R. 
and g x > 0 for the mass matrix cr z ®T x \ m, € R and g, > 
0 for the mass matrix a z ®t z , 0 and g A > 0 for the vector 
potentials; and 0 and g v > 0 for the scalar potential; 
respectively. With this convention, the parameter space 
for the phase diagram is no less than eight dimensional 
with the point 


c. Symmetry class D 

We consider the 2D random Dirac Hamiltonian 

H := —i a x d x — icr y d y + ma z (B13a) 

of rank r = 2 in the symmetry class D with the operation 
for charge conjugation represented by 

C := a x K, (B13b) 

where K denotes complex conjugation. Disorder only en¬ 
ters through the mass matrix mcr z owing to the sym¬ 
metry constraint imposed by the PHS for the symmetry 
class D. 

The OPE of the fermion quartic term, 

n 

: = - ■'i’ a VzV4’ b ( B14 ) 

a, 6=1 


0 = m = g 0 = m x = g x = m z = g z =g A =g v (B18) 

corresponding to the massless clean Dirac Hamiltonian 
of rank r = 4 in the 2D symmetry class D. 

It is natural to start with the clean limit 

0 = g 0 = g* = g z = 9 a = 9v- (B19a) 

In this limit, we may write [compare with Eq. fl3T0l )] 

U = Ho* d x - icr y d y ) <8> r 0 + V, (B19b) 

where the Dirac mass V is decomposed into 

V = m^ 0 + V 0 (B19c) 

with 

Pq := 0’z G T o, V 0 :=fr Z ®M 0 , M 0 :=rr x T x + m z T z . 

(B19d) 
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The Hermitian 2x2 matrix M 0 has the N = 2 non¬ 


degenerate eigenvalues 


m_ 


:= ±\/mf 


m: 


(B20) 


provided m 2 + m z > 0. The four eigenvalues of V are 
m ± + m and m ± — m. The relevant parameter space 
needed to identify the distinct topological phases is two- 
dimensional. It is spanned by the parameters m and m ± . 
There are three topologically distinct insulating phases 
that are defined by 


(B21a) 


+1) 

m > m + , 

0, 

-H 

V 

T 

-1, 

m < m_, 


where the index v is given by [recall Eq. (3.8)] 


v = ^sgn(m + |m + |) + ^ sgn (m - |m + |) . (B21b) 

The first term on the right-hand side is the dimension¬ 
less (thermal) Hall conductivity of a 2 x 2 massive Dirac 
Hamiltonian with the mass m + |m + |. The second term 
on the right-hand side is the dimensionless (thermal) Hall 
conductivity of a 2 x 2 massive Dirac Hamiltonian with 
the mass m — |m + |. A phase transition between the 
|zx| = 1 and v = 0 insulating phases occurs whenever 

(B21c) 


m = m 


phase boundaries (B21c) in the two-dimensional pa¬ 
ster space, 

{(m, m + ) € M 2 | m + > 0}, (B22) 


The 
rameter space 


meet at the massless Dirac point m = m + = 0. The 
massless Dirac point m = m , = 0 is unstable to any non¬ 
vanishing m or m + . The line parallel to the horizontal 
axis in Fig. [6](a) intersects the phase boundaries (B21c) 
at the two critical points m = ±m + at which the gap of 
one and only one of the two massive Dirac cones in the 

spectrum closes. _ 

The origin of parameter space ( |B22 ) is a multicritical 
point, contrary to the two critical points at fixed value 
of m + depicted by bullets in Fig. Ifita). This difference 
manifests itself in the stability of these three points to 
disorder. Our next task is to assess the effects of disor¬ 
der on the multicritical point at the origin of parameter 
space (|B22h. 


If we choose the probability distribution with 

0 = m = m, = g x = m z = g. = 0, g,£Jv,2A>0, 

(B23) 

then the random Dirac Hamiltonian is reducible, as a 
rotation of r y about t x into r z brings it to the represen¬ 
tation 


"H** = + (A x a x + A y <j y ) <g) t z 

+ M a z ®T 0 + V a 0 ® t z . 


(B24) 


The random Dirac Hamiltonian "H** of rank r = 4 is 
evidently block diagonal with irreducible blocks of rank 
r = 2 that belong to the symmetry class A. Hence the 
variances g 0 , g v , and g A obey the same one-loop RG 
flows as in Eq. (BIO) that we depict in Fig. |6][b). In 
particular, the variance g 0 of the random mass that com¬ 
mutes with all other masses flows to the strong coupling 
as soon as either g v or g A is non-vanishing. 

The reducibility of the rank r = 4 Dirac Hamiltonian 
in the 2D symmetry class D is non-generic. It is broken 
as soon as any one of the two variances g x or g z for 
the anticommuting pair of Dirac mass matrices is non¬ 
vanishing. In the parameter space 


0 = rrn = m_ 


> So 7 


,9ai9v> 0 ^ ( B25 ) 


the variance g is then expected to inherit the marginal 
relevance of the multicritical point (B18) in the neigh¬ 


borhood (B23). If the multidimensional parameter space 


is projected onto a two-dimensional plane parametrized 
by the probability distribution ( |3.9a ), we conjecture by 
analogy to Figs. 3](a) and|4][b) the phase diagram shown 
in Fig. 6j[c) . This phase diagram is not the only one 
dictated by logic. It has the merit of simplicity, how¬ 
ever. For example, its RG flows are more economical 
than conjecturing that the metallic phase precludes any 
direct phase boundary between two topologically distinct 
localized phase for any non-vanishing variance of the ran¬ 
domness. 

The stability analysis of the two critical points inter¬ 
secting the phase boundary at a non-vanishing fixed value 
of m + in Fig. [6j(a) in the presence of a generic disorder 
from the symmetry class D is quite different from that 
of the multicritical point. We assume that the fluctua¬ 
tions are much smaller than the means at the two crit¬ 
ical points defined by |m| = m + with m + > 0 given 
in Fig. [6j[a) . This assumption justifies projecting the 
Dirac Hamiltonian of rank r = 4 onto the eigenspace 
of rn^T x +mZT z with the eigenvalue m + . 1 ^The resulting 
projected Dirac Hamiltonian is of rank r = 2 and has 
the same form as the rank r = 2 Hamiltonian defined 
in Eq. (B13). We can then apply the RG analysis for 
the Hamiltonian of rank r = 2 in Eq. (B13) to find that 


the effective coupling playing the role of g M in Eq. (B13) 
is irrelevant along the phase boundary defined by the 
gap closing condition m = m , . This reasoning explains 
why a non-vanishing mean value m + > 0 is required in 
Eq. (3.10) in order for one of the N clean critical point 
at g = 0 to be the end point of a phase boundary that 
separates a pair among N + 1 topologically distinct lo¬ 
calized phases. This end point is the critical point that 
governs criticality at the phase boundary between two 
topologically distinct localized phases. No metallic phase 
prevents a direct transition between this pair of topologi¬ 
cally distinct localized phases for not too strong disorder. 
The phase diagrams when N = 1 and N = 2 are depicted 
in Figs. [4](c) and |4|]d ), respectively. A phase diagram sim¬ 
ilar to that of Fig.Wc) has been obtained from numerical 
simulations of network models in Refs. m and [S2j. 
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We close the discussion of the 2D symmetry class D 
with a generic r = 4 random Dirac Hamiltonian by ob¬ 
serving that the Dirac mass 


respectively. 

There are two random functions allowed by the sym¬ 
metry constraints that multiply the matrices 


V 0 (r) :=m x {r) a z 3 t x + m z (r) a z 3 t z 

= m + (r) [cos0(r) a z <g> t x + sin 0(r) a z <8> r„], 

(B26a) 


where 


m + (r) := V m 2 x {r) + m 2 z {r ), 9(r) := 

m x {r) 

(B26b) 

is generated by two anticommuting mass matrices, i.e., 
Vq (r) = m 2 + {r ) a 0 3 r 0 . (B26c) 


The Dirac mass ( B26| ) supports point defects defined by 
the two conditions that (i) linv^oo m,[r) = m + and 
(ii) the phase 0(r) winds by an integer multiple of 27 t 
along any closed path T r belonging to any open disk 
centered about the point r vtx that is not too large. The 
point r vtx thus represents a point-like singularity of the 
vortex type. Jackiw and Rossi have shown in Ref. [1141 
that Dirac Hamiltonians of the form 


^vtx : = ^* + KtxM 


(B27) 


support n normalized zero modes if V vtx (r) is given by 
Eq. (B26) with the phase 8(r) winding by 2nn along 


all circles of sufficiently large radii. These zero modes 
are not robust to all remaining random channels com¬ 
patible with the 2D symmetry class D. However, if n 
is an odd integer, at least one zero mode survives any 
perturbation arising from the scalar gauge, vector gauge, 
and the last massive channel, provided this massive chan¬ 
nel associated to the Dirac mass matrix that commutes 
with V vtx (r) is everywhere smaller in magnitude than 
the asymptotic gap m + . This is an illustration of the 
fact that the fundamental homotopy group for the nor¬ 
malized Dirac mass matrices in the 2D symmetry class 
D can be Z 2 . Correspondingly, the localized phase v = 0 
in the phase diagramed) is conjectured to support Z 2 
quasi-zero modes that will contribute to the density of 
states. 


d. Symmetry class Dill 

We consider the 2D random Dirac Hamiltonian 
H := —ia x 3r 0 d x — i(7. y 0T o d y + ma z 3r y + V <j 0 3t x 

(B28a) 

of rank r = 4 in the symmetry class Dill with the oper¬ 
ations for reversal of time and charge conjugation repre¬ 
sented by 

T := \cr y <3 t 0 K, C := a x 0 t z K, (B28b) 


cr z 3 T yi (random mass) 

a 0 3 t x . (random scalar potential) 


(B29) 


These couplings are the coefficients of the fermion quartic 
terms 


$v = - X! : a 0 ® T x V cr 0 3 t x ip b :, (B30a) 

a,6=1 
n 

: ^° a z® T y' , P a ^ b ^ z 3r y ip b :, (B30b) 

a,6=1 


respectively. Their OPEs are 

^m(O) = 2 2 ^m(0)i (B31a) 

7r z r z 

<Mr)<M 0) = -^$ v (0), (B31b) 

7 t z r z 

<M*0 <f> M (0 ) = ^ [<M0) - 4> m ( 0)] . (B31c) 


We note that (i) the OPEs (B31al and (B9a) are identi¬ 
cal, (ii) the O PEs (B31b ) and (B9b| are identical, while 
(iii) the OPE (B31c) differs by an overall sign from the 
OPE (B9c) if one ignores <& A in the latter OPE. The 
following one-loop RG equations follow from this obser¬ 
vation, 


Pm ~ ~ 9m (9v — 9m) > (B32a) 

7T 

Pv — — 9v (flV ~ 9m) ■ (B32b) 

7T 

Addition and subtraction gives 

(B33a) 

d(9 '7/ M) = 1 ( Sv-Suf. (B33b) 

ai 7r 


We have found the marginally relevant linear combina¬ 
tion (B33b) of the couplings associated to the disorder. 
We deduce that a metallic phase must separate the only 
two topologically distinct localized phases in class Dill 
to explain the multicritical nature of the clean critical 
point g = m = 0 as depicted in Fig. |4|e), in agreement 
with numerical simulations of a network model.® 3 


e. Symmetry class All 

The one-loop RG analysis of a Dirac Hamiltonian of 
rank r = 4 in the symmetry class All can be found, 
e.g., in Ref. 103 . Here we give a brief summary of 
the relevant results. The r = 4 Dirac Hamiltonian 
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(a) 

m + 




(c) g D (N=2, m + =0) 



FIG. 6. (a) Phase diagram in the parameter space (B22). There are three gapped phases all meeting at the origin, the clean 

massless Dirac Hamiltonian of rank r = 4 in two-dimensional space. The origin of parameter space is a critical point that is 
unstable to both m or m + as is indicated by the arrows. The gray region with the mass m > m ± > 0 supports the dimensionless 
quantized thermal Hall conductivity v = +1. The gray region with the mass m < — m + < 0 supports the dimensionless 
quantized thermal Hall conductivity v = — 1. The white region with m + > |mj > 0 support the dimensionless quantized 
thermal Hall conductivity v = 0. The critical lines m = ±m + separates two insulating phases whose Chern numbers differ by 
the number 1 in magnitude. The line parallel to the horizontal axis at fixed value of m + intersects the phase boundaries at two 
critical points depicted by two filled circles, (b) Phase diagram in the parameter space <7 m>Sv>5a ^ 0 spanned by the variances 
of a random mass, a random scalar potential, and a random vector potential (all of vanishing means) in the 2D symmetry class 
A for a random Dirac Hamiltonian of rank r = 2p^ The bullet represents the stable flows to the critical point that describes 
the plateau transition in the integer quantum Hall effect. The unstable flow away from this critical point corresponds to a 
non-vanishing mean value of the random mass and is thus not present in this three-dimensional parameter space, (c) Projected 
phase diagram with the vanishing mean masses m + = 0 onto a two-dimensional cut with the same horizontal axis as in (a) 
but with a vertical axis encoding multiples sources of disorder including the variance g associated with the random mass of 
mean m. The direction g can be thought of as a cross section of all the directions associated with the variances of the disorder 
coming out of the plane in (a). 


with the operation of reversal of time represented by 
T := i <7 y g t 0 K has symmetry-allowed random pertur¬ 
bations th at inclu de those allowed in the symmetry class 
Dill [Eq. (B28a)]. Therefore, the set of one-loop RG 
equations for the symmetry class All, which is given in 
Ref, rm-n contains the one for the symmetry class Dill 


as a subset. Thus, we have the same linear combination 
of coupling constants as in the symmetry class Dill that 
is marginally relevant [Eq. (B33)]. This indicates that a 
metallic phase always intervenes between two topologi¬ 
cally distinct localized phases in class All, as is shown in 
Fig. [4|[g) jiniHi05li28l 
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